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ABSTRACT 

Laser radar observations of the nesospherlo sodlia layer often reveal 
wavelike density fluctuations moving through the layer. It Is often assumed 
that these features are a layer density response to gravity waves. The 
density response of atmospheric layers to gravity waves Is developed In two 
forms; an exact solution and a perturbation series solution. The degree of 
nonlinearity In the layer density response is clearly described by the series 
solution whereas the exact solution gives insight Into the nature of the 
response. It Is shown that density perturbations In an atmospheric layer can 
be substantially greater than the atmospheric density perturbations associated 
with the propagation of a gravity wave. Because of the density gradients 
present In atmospheric layers i Interesting effects are observed such as a 
phase reversal in the linear layer response which occurs near the layer peak. 
Once the layer response Is understood y the sodium layer can be used as a 
tracer of atmospheric wave motions. 

In order to Increase the resolution of Ildar data^ a two-dimensional 
digital signal processing technique has been developed. Both spatial and 
temporal filtering are utilized to enhance the resolution by decreasing shot 
noise by more than 10 dB. Many of the features associated with a layer 
density response to gravity waves are observed In high resolution density 
profiles of the mesospheric sodium layer. These include nonlinearities as 
well as the phase reversal In the linear layer response. 
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1. INTRODOCTION 

The availability of laaera haa aade It auoh easier to take advantage of 
resonant scattering In the reaote detection of atmospherlo oonstltuenta. One 
of the Initial reaote sensing applications of tunable dye lasers was In the 
study of the aesospherlo sodlas layer* The layer has been Investigated since 
Its discovery In the late 1930's* A review of the early measurements and 
theory of the aodlas layer ohemistry is contained In Brown [19733* 
NesosiAierio sodlm Is believed to be of meteorio origin [Richter and Seohristt 
1979a; Gadsden, 1970]. In addition to providing a source of atomic sodliai, 
recent work indicates meteoric dust may also contribute to the removal of 
atomic sodlw from the mesos;4iere [Hunten, 1981]* 

Before dye laser based Ildar systems became available as a tool for remote 
sensing, sodlw measurements were largely restricted to studying resonantly 
scattered sunlight* Ground-based measurements of this type [Blamont and 
Donahue, 1961; Gadsden and Purdy, 1970] were able to define seasonal 
variations In colimin abundance, but the sharp layer boundaries were not 
revealed until rocket-bom dayglow measurements were guide [Hunten and Uallaoe, 
1967; Donahue and Meier, 19671* Lldar observations of the vertical structure 
of the sodium layer were first ggade In England [Bowman et al*, 19693* Since 
then, similar measurements have been reported from a variety of locations 
including Frgmce [Blamont et al., 1972], Brazil [Kirohhoff and Clmiesha, 
1973], California [Hake et al., 1972], Illinois [Rowlett et al., 1976] and at 
hl^ latitudes (Franz Joseph Land, USSR) [Juramy et al., 1981]* Because the 
daytime sl^ produces unacceptably hl^ backgroimd photon rates In most sodlin 
Ildar receivers, alggcst all Ildar observations of the sodlw layer have been 
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restricted to the nighttime. In addition, most lidar observations of the 
sodlia layer are zenith neasurements* Steerable measurements have been made 
[Thomas et al., 1976; Clemesha et al., 1980] in an attempt to determine the 
horizontal structure of the layer and the direction of propagation of mavellke 
features in the layer. 

The mesospherio sodium layer is generally confined to an altitude range of 
80-100 km with a sharply defined peak near 90 km where the maximum number 
density is usually about 10^-10^^ m**^. The upper boundary typically has a 
scale hei^t that is less than that of the atmosphere (3-4 km compared to 
about 7 km for the atmosphere at this altitude). The lower layer boundary is 
even more sharply defined; the sodlta density dro;>s an order of magnitude 
wittiin a few km at about 80 km in altitude. Seasonal as well as daily and 
hourly variations in the layer are observed. A winter increase in column 
abundance has been observed in both dayglow [Hunten, 1967; Gault and Bundle, 
1969] and Ildar [Gibson and Sandford, 1971; Megle and Blamont, 1977] 
observations. Current theories of sodiiai layer chemistry [Richter and 
Seohrlst, 1979b] appear to explain the annual abundance variations. Diurnal 
variations in oolunn abundance were observed In ground-based [Gadsden and 
Purdy, 1970] and rooket-bom [Donahue and Meier, 1967] dayglou measurements. 
Daytime lidar measurements [Gibson and Sandford, 1972] have failed to confirm 
these observations. ..jwever, a pre-sunrise increase in column abundance has 
been detected in nighttime observations made at Urbana, Illinois and is 
reported in this thesis. These variations are not fully explained by current 
theories of the sodium layer. 

Short term variations in sodium density are often observed [Kirohhoff and 
Clemesha, 1973; Rowlett et al., 1978; Juramy et al., 1980] and are attributed 
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to tidal and gravity nave osolllatlona. To understand the effeota of 
atmoapherlo dynanioa on the aodlun layer , the response of layers to 
atoospherio wave notions nust be understood. Fortunately, ohenloal effects 
are reduced in the nighttime sodlvn layer so that it oan more readily be used 
as a tracer of atnospherlo notions. The linear density response of minor 
constituents to gravity wave perturbations has been investigated by several 
researchers [Dudis and Reber, 1976; Chiu cuid Chlng, 1978]. These results have 

been applied to 02 ('j;) airgloif measurements [Helnstook, 1978], measurements of 
the green line alrglow at 557*7 nm [Battaner and Mollva, 1980], and the 
nightglou emission of the hydroxyl radical [Frederick, 1979]* It also appears 
that many of the wavelike features observed in the sodlun layer are associated 
with a linear layer response to atmospheric waves [Shelton et al., 1980]. An 
expression for the total layer response to atmospheric waves Is derived in 
this thesis. Many of the features in experimental data are consistent with 
such a layer density response to gravity waves. Features often observed in 
the linear response of the layer Include a phase reversal emd differing 
perturbation amplitudes above and below the layer peak. Nonlinearities in the 
layer response have also been observed. In regions of large density gradients 
or where the linear response is small (i.e. at the point of phase reversal), 
nonlinear effects become important. Sodlun density perturbations associated 
with a double frequency response to gravity waves are consistently observed in 
these regions. Because of the nonlinear nature of the sodlvin layer response, 
density oscillations with short temporal periods and small vertical 
wavelengths can be present. 

A detailed study of the nonlinearities in the layer response to gravity 
waves requires high temporal and spatial resolution. Although shot noise 
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places a practical limit on Ildar resolution, appropriate filtering can reduce 
the shot noise and significantly Improve the resolution In Ildar profiles. 
The spatial filtering technique of Rowlett and Gardner [1979] Is expanded In 
this thesis to Include temporal filtering. As a consequence, spatial and 
temporal resolution have been inoreased and many of the nonlinear features in 
the layer response to atmospheric waves have been observed. 

In Chapter 2, the Urbana Ildar system is described and the principles of 
Its operation are presented. Chapter 3 contains a description of the spectral 
analysis and two>dlmenslonal filtering teohnlques utilized to process the 
received photocount profiles. The sodium layer response to atmospheric wave 
motions Is developed In Chapter 4. In addition to explaining many of the 
wavelike features observed In sodium density profiles, the layer response 
Influences colvmm abundance and must be understood in Interpreting the spatial 
perlodograms used In data processing. Observations of the sodium layer made 
at Urbana, Illinois between October, 1979 and March, 1981 are presented In 
Chapter 5. The wavelike perturbations In density found In many of these data 
exhibit the characteristics of a layer response to gravity waves developed In 
Chapter 4. All of these measurements were restricted to zenith observations 
and, as a result. It Is difficult to deduce the horizontal structure of the 
sodium layer* Since the horizontal layer structure can be Influential In 
determining the layer response to atmospheric waves It la sometimes difficult 
to Interpret Ildar data when only the vertical variations In sodium density 
are known. In Chapter 6, Initial results of steerable observations of the 
sodltai layer are shown, niese observations, made at Goddard Space Pll^t 
Center, Illustrate a technique which may be used to determine the horizontal 
structure of the layer. Chapter 7 describes the stratospheric dust layers 
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observed over Urbans following the voloanlo eruptions of Mt. St. Helens In 
1980. Althou^ these observations are not dlreotly related to observations of 
the aesospherlo sodlia layer 1 they Illustrate the flexibility and the aany 
uses of ground-based Ildar facilities. Conoluslons and reoosmendatlons for 
further work are outlined In Chapter 8 . The appendix contains a description 
of sene of the Ildar systen hardware. 
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2. LIDAR SYSTSI DESCRIPTION AND OPERATION 


2.1 OVERVIEW 

The Urbane sodium Ildar facility consists of a Fresnel lens telescope and 
the associated receiving electronics, a pulsed organic dye laser, and a 
microprocessor system which coordinates the data collection process and stores 
the received data. The output of the dye laser is tuned to the sodium D2 line 
at 5890 % In order to Induce fluorescence In the atomic sodium of the 
mesospheric sodium layer. This Increases the effective backscatterlng cross 
section of atomic sodium substantially. Baokscattered photons are collected 
by the Fresnel lens telesccpe where the detected photons are counted. By 
periodically sampling these counts, a record of photocounts versus height is 
obtained tdilch may be used to determine estimates of the density profile of 
sodium with altitude. The processing required to obtain these estimates Is 
performed off-line on the University of Illinois Control Data CYBER 175 
computer. 

2.2 RECEIVING SYSTEM 

The major elements comprising the receiving systmn are a Fresnel lens 
telescope, photomultiplier tube (PMT) and Its associated electronics, and a 
counter Interface. Each of these Items Is Included In the Ildar system block 
diagram (Figure 2.1). 

2.2.1 RECEIVING TELESCOPE 

The collecting element of the receiving system telescope Is a 1.22 meter 
diameter f/1.56 acrylic Fresnel lens. The lens as well as the remainder of 
the telescope are discussed In some detail by Rowlett and Gardner [1979]* A 





Figure 2.1 


Blook diagraa of the Urbane Ildar ayatea. Double llnea Indicate 
light patha %ihlle alngle llnea repreaent eleotrloal oonneotlona. 
Angled llnea In light path ayabollee beaMplittera and adrrora. 









8 


block diagraa of th« reoeiving tolesoopo la shown In Figure 2.2. The Fresnel 
lens foouses the collected photons <Mito the plane of a field-stop iris. 
Following the iris, a Nikon 35 am f/1.4 lens oolliimtes the beam prior to its 
passage through an interference filter. Because of the 5 nm bandpass of this 
filter, the number of photons due to background sky illumination is reduced 
substantially. A focusing lens contained in the PHT cooled housing foouses 
the photons which are passed by the interference filter onto the oathode of an 
RCA C3 1034 a high quantun efficiency PHT. This PHT is thermoelectrloally 
oooled by a Products for Research Hodel TE-206TS-RF cooled housing. Table 2. 1 
contains a list of parameters for the telescope. All data presented in this 
thesis, with the exception of the measurements presented in Chapter 6, were 
collected with this telescope. Heohanioal considerations have restricted the 
use of the telescope to zenith measurements. 

TABLE 2.1 RECEIVINO TELESCOPE PARAMETERS 


Diameter of objective 1.22 m 
Area of objective 1.17 m^ 
Focal length of objective 1.90 m 
Bandwidth 5 nm 
Field of view 3 nrad 


2.2.2 PHT ELECTRONICS 

The output of the PHT is discriminated by a Princeton Applied Research 
Amplifier-Discriminator (Hodel 1121). The discriminator compares the 
amplitude of pulses generated by the PHT to a preset threshold. If the pulse 
exceeds the threshold, a 12 nanosecond pulse is generated which is counted by 
the counter interface. Discriminator output pulses are NIH standard, l.e. -18 
ma mininim into a 50 ohm line terminated in a 50 ohm load. Another mode of 
operation for the discriminator is the Pulse Height Analyzer (PHA) node. The 



(RCA CSI054A) 


Figure 2.2 SoheMtlo diagraa of telesoope optloa [Rowlett and Gardner 
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prlnolpal utility of the PHA node Is that the aaplltude distribution of pulses 
arriving fron the PMT oan be displayed pemltting the determination of the 
proper threshold level. In addition) this mode oan be used to optimize the 
PHT high voltage setting in order to maximize the output signal to dark count 
ratio. 

There are several possible souroes of error in photon counting systems 
uhic : may introduce significant error at very high oounting rates. One suoh 
source is "pulse pile-up". The PMT output pulse has finite width) usually on 
the order of 10-40 nanoseoonds. As a result) the PMT oan aot as a paralyzable 
detector in that two or more output pulses separated by less than the 
resolving time of the tube will result in one larger and generally wider 
pulse. Another source of error arises from signal induced noise. It hats been 
shown that such noise is proportional to signal strength and inversly 
proportional to t ) where t is the elapsed tine from the signal pulse 
[Pettlfer) 1975]. To minimize these problems and to protect the PMT from 
overloading conditions which are possible due to strong Rayleigh scattering 
fron below 30 Ion altitude) a gain switching or gating capability has been 
added to the receiving system. Prior to its installation) the receiving 
system overloaded at low altitudes and required at least 100 microseconds to 
recover. As a result) Rayleigh scattering fron 30 to 40 km oould not be 
observed. Since the calibration procedure requires these datS) the sodlimi 
data could not be calibrated to yield estimates of the absolute sodium 
density. 

Gating of the PMT is achieved by controlling the voltage on the first PMT 
dynode. When the dynode voltage is reduced fron its normal state (which is 
positive with respect to the cathode) to the cathode potential or slightly 
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negativt with respeot to tho oathodo potential » a aignif leant daoroaao In 
overall PMT gain la aohleved. Measureaenta with thla ayaten indicate that the 
gain la deoreaaed by approxlnately 3 ordera of aagnltude. In addition » the 
recovery tine la abort , on the order of 6 to ^0 nlcroaeoonda. Thua, the PMT 
■ay be electronically awltohed on and off rapidly enough to permit obaervation 
of Rayleigh aoatter froB 30-40 ka without aufferlng from overloading aa a 
reault of atrong aoat taring froa below 30 ka. 

In order to aynohronlze the firing of the laaer and the PMT gating* a 
timing controller waa oonatruoted. Thla element oontalna a triggering circuit 
for the laaer and the blanking olroulta aa well aa a free running oaolllator 
which eatabllahea the laaer pulae repetition rate. The timing la adjuated ao 
that the PMT la blanked 5 mloroaaoonda before a laaer pulae la initiated. 
Thla Inaurea that the PMT will be ooapletely blanked Immediately after the 
laaer flrea tdien the atrongeat return aignal ia preaant. Approximately 200 
mloroaeoonda later (oorreaponding to an altitude of 30 km)* the PMT la 
unblanked and data collection prooeeda. Both the blanking and the timing 
olroulta are daaoribed in greater detail in the appendix. 

Pigurea 2.3. and 2.4 abow idiotooount prof ilea trtiioh illuatrate the 
effeotlveneaa of PMT blanking. In Figure 2.3* overloading haa aeverly 
dlatorted retuma from below 40 km and data calibration* which dependa on 
Rayleigh retuma fi^ 30 to 40 ka* la not poaalble. In Figure 2.4* returns 
from below 30 ka have been attenuated by blanking. Above 30 km* however* 
there la auffloient Rayleigh acatterlng to permit calibration of the apatial 
profile and therefore eatlaatea of the abaoluta aodlia denaity are poaalble. 
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2.2.3 COONTBR INTERFACE 

The disorlalnated pulaes fron the amplifler-dlsorlminator are counted by a 
free running b bit counter in the counter Interface. Thia device , deaoribed 
by Kinter [1977]» containa the logic neceaaary to periodically aample and 
atore the contenta of the counter. Counter interface operation ia 
aynchronized with the laaer firing by a pulae fron a photodiode detector (the 
Laaw Cconand Pulae > LCP). The counter aanpling Interval • or range gate 
width) ia aoftware controlled and nay be varied from 1 to 128 nioroaeoonda in 
1 nioroaecond icorenenta. All the aodiun data preaented in thia theaia were 
collected with a 3 nicroaecond range gate. In addition to the range gate^ a 
delay factor auat be apeoified. Thia deterainea the delay between laaer 
firing and the point at which counter aanpling beglna. Becauae only 256 
f^otooount aaaples are available) auoh a delay may be required to insure that 
the sodium layer is included in the sampled region. A 12 bit word ia used to 
communicate these parameters aa well as commands from the lidar microprocessor 
aystem to the counter interface. Any return data) such as photooount samples) 
are returned in an 8 bit word. Two additional lines provide handshaking 
between the counter interface and Uie mioroprooesaor system. 


2.3 HICROPROCESSOR SYSTEM 

A Digital Group Z-80 microprocessor system forms the heart of the 
nloroprooeasor system. This system %ias initially implemented to control the 
data collection process) provide a control system tdiloh insured that the laser 
tuning remained stable) and to monitor the laser spectral linewidth and output 
energy [Teltlebaum and Seohrist) 1979]. As the lidar system has evolved) so 
has the function of the microprocessor system. Currently) the microprocessor 
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oontrola data oolleotlon and provides a peraanent record of the data on 
oaaaette tape. 

The Z-80 ayatea peraita allgnamt of reoelver and laaer, oontrola data 
oolleotlon and atoragO) and tranaalta the acquired data to the CDC CTBBR 175 
ooaputer for off-line prooeaalng. Rayleigh acatterlng and returna frca the 
aodlun layer are diaplayed to facilitate the allgnaent proceaa. During the 
data collection process » the range bln» delay and nunber of laser pulses per 
profile are specified by the user. This pemlts a great deal of flexibility 
In system operation. After each spatial profile of data la collected , It la 
stored on cassette tape. Thus, a power or equipment outage will never result 
In the loss of more than one profile of data. Data Is transmitted to the 
CYBER 175 computer via a telephone line at 300 baud. Vhen operated In this 
mode, the Z-80 system emulates a remote terminal tied to the CYBER 175. 

2.4 LASER 

Two lasers were used to collect the data presented in this thesis. Both 
lasers are flashlamp pumped, organic dye lasers utilizing Rhodamlne 6G 
perchlorate In an alcohol and water solution. Photographs of the two systems 
are Included In Figures 2.5 and 2.6. Prior to January, 1981, a laser 
utilizing a double elliptloal cavity was used [Richter and Seohrist, 1978]. 
The dye cavity was at the oommcm focus of the two ellipses and a flashlamp was 
located at each of the remaining fooll. During operation, the quartz walls of 
the flashlamps were ablated. Because of a build-up of ablation products, the 
flashlamp lifetime was limited to approxlmatexy 5000 p<ilses. The pulse 
repetition rate for this laser was only 0.5 hertz. This rate could not be 
exceeded without Jeopardizing laser stability. 
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When this laser was raplaoed in January, 1981, the quality of data 
laproved draMtloally. This was partially a result of an Increase In the 
laser pulse repetition rate and the extended life of the new flashlamps* The 
new laser Is a Candela Corporation Model LPDL-5 laser whloh has been Modified 
to Bieet a special narrow llnewidth requlreaent. The narrow llnewldth 
nodlfloatlon has not been used In data oolleotlon activities . It Is 
anticipated that this feature will be used In the future to attempt 
temperature measurements of the sodium. Table 2.2 oontalns parameters 
describing both the old and new lasers. 

TABLE 2.2 LASER PARAMETERS 



Old Laser 

Candela Laser 



Broadband 

Narrowband 

Energy per pulse 

100 mj 

50 mJ 

20 mJ 

Pulse repetition rate 

0.5 Hz 

10 Hz 

10 Hz 

Wavelength 

589.0 nm 

589.0 nm 

589.0 nm 

Llnewldth 

6 pm 

1 pm 

0.2 pm 

Pulsewldth (FWHM) 

4 ps 

2 ps 

2 PS 

Beam divergence 

1 mrad 

1 mrad 

1 mrad 


The operation and characteristics of the Candela laser as well as a means 
of measuring the spectral llnewldth of the laser output are described by 
Nllssen [1981]. The Improved system performance Is explained by the Increased 
pulse repetition rate, extended flashlamp lifetime, and Increased spectral 
energy density at 5890 i. The spectral energy density of the old laser was 
approximately 100 mJ/6 pm or 16.7 mj/pm whereas the Candela laser can produce 
50 mJ/pm In the broadband configuration. In addition, the spatial resolution 
Is Improved somewhat by the shorter output pulse width (2 microseconds as 
opposed to 4 microseconds). The Increased pulse repetition rate available 
with the Candela laser (10 Hz) makes greater temporal resolution possible by 
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decreasing the amount of time required to acquire a spatial photocount 
profile. Actually) the effective laser pulse repetition rate is Halted to 
approximately 4 hertz by the Z»80 system. This system requires roughly 2b0 
milliseconds to collect a laser shot of data and update the video presentation 
summarizing the status of the spatial profile. 

One of the consequences of the extended flashlaap lifetime is that long 
data sets may be collected without dismantling the laser to replace the 
flashl<uiip8 (as was the case with the old laser). The current system has 
consistently run from sunset to sunrise with no major tuning or equipment 
problems. These long data setS) which contain as many as 128 spatial 
profiles I make temporal filtering an attractive addition to the signal 
processing techniques tdilch are applied to the photocount data. These 
processing techniques are described in Chapter 3. 


2.5 DESCRIPTION OF THE PHOTOCOUNT DATA 

The lidar equation relates the received photocounts due to laser induced 

fluorescence of the sodium layer to the sodicn concentration profile. The 

contributions of the various Ildar system parameters ais well as atmospheric 

effects are included. The lidar equation derived by Measures [19771 has been 

modified to Include temporal variations in sodium density. 

00 2 

y(ott) ( 2 . 1 ) 

( 2 . 2 ) 


A(s,t) » 1 j ~ ° 


0 


yU,t) 


6(a.t) 


exp 


-2 I deaia) \ 
0 
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Where X « photoeleotron counting rate (s~^) 

* scattering profile (m~^) 
l - rj^exp(-2s/c3’£), siO (a“^) 

s radiative lifetiae of the soatterer (a) 
p s laser pulse shape 
n = overall system effioienoy 

2 

Aj^ s receiver aperture area (n ) 

h s Planck's constant (J*s) 

V = optical frequency (s“^) 

a s volume extinction coefficient (m*^) 

^ s volume backscatter coefficient (m**^) 

The backscatter coefficient is proportional to the concentration of sodium and 
is therefore the parameter of Interest. Note that this parameter is dependent 
upon both altitude and time. The altitude dependence is a consequence of the 
vertical structure of the sodlm layer and layer dynamics account for the 
temporal dependence. The volume extinction coefficient is negllgable for 
heists greater than 30 km. Hence y we have 

y(z»t) « (2.3) 

for the region of the sodium layer. If i(a) and y(z,t) are assumed to be zero 
for a>0 I the range of integration for both integrals in equation (2.1) may be 
extended to infinity. This results in an expression tor X(ayt) which is a 
double convolution. 

X(a,t) ■ Hz) * p(z) * Y(a.t) 


(2.4) 
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In this case I l(s) and p(s) have the effect of oonvolutionally aaearlng 

The above results apply to the returns from a single laser pulse* 
However I the returns fros 50 to 500 pulses are Integrated to yield a single 
spatial profile of data. Alsot the Ildar systes does not oMasure the 
Instantaneous photon rate. Instead) It Integrates tUe rate over the spatial 
variable) 2 ) by counting the deteoted photons In each range bln to yield a 
photooount) X. The relation between the rate parameter of the counter 
Interface (which Is the counter sampling period) T) and the range) Ah, over 
which the photons are Integrated Is bh^oTll where o is the speed of light. 

If the counting rate Is deterministic) then x Is Poisson and the expected 
value of the photooounts obt 2 dned from one laser pulse oeourlng at time t is 

a+Ah 

<x(a,t)> » I da\(8tt) + <b> (2.5) 

2 

or 


<x( 2 ,t)> ■ <x^( 2 ,t)> + <b> 


( 2 . 6 ) 


where <b> Is the expected photocount due to additive noise and<Xy( 2 ,t)> Is the 
expected signal photocount. The Integration In equation (2*5) nay be 
expressed as the convolution of X( 2 ,t) with a gate function) giz)t where the 
gate function Is defined ast 

(k J 0 < 2 < Ah 


giz) - \ 


(2.7) 


0 ; elsewhere 


When combined with equation 


(2.4)) 


the expression for the expected signal 
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photooounts is ssen to contain a thrsefold oonvolution* 

^ A(s) * p(a) * g(a) * yia,t) (2.8) 


or 


<x^(a,t)> - /(a) * Y(a,t) (2.9) 

r\Ap 

where f(a) • (-r— ) i(a)<*p(a)*^(a) . In the spatial Fourier domain^ this aay be 

hv 

written as 

■ F(wp<r(w^,t)> (2.10) 

where <x^> $ <T> and p are the spatial Fourier transfoms of x^t y and ft 

respectively, u is the spatial frequency. In practice • the returns from a 

3 

number of laser pulses are summed to yield a spatial {diotooount profile. 
These sums are also Poisson. The loss of temporal resolution resulting from 
this integration is discussed in Chapter 3. The photocounts available from 
the counter Interface are the samples of x(s,t) at the points z^kLh* Tbus» 
the photocounts for the range bin and the spatial profile may be 

expressed as x(kLh,m^) tdiere aT ie the time interval separating spatial 
profiles. Since x{kLhtmLT) is Poisson ^ the expected value and variance are 
equal: 

<x{.kLhtmLT)> - VartxOcAft.mAjT)] (2.11) 

The noise photooount contained in equation (2.5) is primarily due to 

background sky illumination and niT dark counts. The expected value of this 
photooount is equal to X f where X is the noise rate and T is the counter 

ft ft 

sampling period. Oenerally« the noise counts due to background sky 
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lllwinatlon aaln the nost slgnif leant oontribution to the total nolae oount. 
The noise souroe nay be reduced by uaing a narrower bandpass Interferenoe 
filter In the reoelvlng telesoope. Beoause a Fresnel lens Is used In the 
telesoopsi It Is not praotloal to reduoe the bandpass of the filter auoh below 
Its current value. Since the telesoope Is not diffraction Halted, 
oolliaatlon errors begin to introduce significant losses and alignaent 
beooaes extreaely difficult. Noise counts due to PMT dark noise are ainialaed 
by cooling the PMT to -20° C. 

Bstiaates of sodiua concentration are obtained by filtering the photooount 
data in such a Banner that the oonvolutlunal effects described by equation 
(2.9) are ainlaized. The filtering process is described in Chapter 3. 
Calibration of the filtered data yields estlaates of absolute sodiua density. 
This process entails noraalizing the data by photooounts arising froa 
Rayleigh scattering and scaling the result to account for the effective 
scattering cross section of sodlvn and the baoksoatterlng coefficient of the 
ataosphere at the altitude of the Rayleigh scattering. The received 
photooounts are related to sodliai density through the following equation 
[Cemy and Seohrist, I960]. 


Pv - 






°e£f 


( 2 . 12 ) 


where 


s soditsB density in the range bin (a~3) 
s ataospherio density at range z (a~^) 

8^ s heli^t of the range bin (a) 

a Signal photooounts froa the range bin 
No a signal photooounts froa the range z 
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« Rayl«lgh baoksoatter ooafflolant (a^) 

o f ■ affaotlva aodlua opoaa saotion (a^) 
exf 

Calibrating the data in this aanner reaulta in aodiua density eatiaates wtaioh 
are independent of ayatea variables suoh as laser pulse energy variations 
and overall systea effioienoy. 


f 
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3. SKMAL PR0CB8SXII0 


3.1 XMTRODOCnON 

Theor«tloally( laaar pula* langtht raoaivar range gat* tildthi data 
integration tiae and saapling interval place the fundaaental liaits on the 
spatial and teaporal resolution of lidar profiles. However f in praotioe the 
resolution of sodiua lidar data is usually United by shot noise. Spatial 
filtering has been shown to be an effeotive teohnique for enhanoing the 
spatial resolution of sodiua lidar profiles by reducing shot noise in the 
photooount data [Hewlett and Gardner « 1979]. Bstinates of sodiua density are 
obtained by spatially filtering the photooount data to renove shot noise and 
to deoonvolve the effects of receiver range gate and laser pulse shape. The 
appropriate spatial filter cutoff frequencies are deterained by coaputing the 
average spatial periodograa (power speotrua) of the profiles oolleoted during 
the ni^ts operation. 

Teaporal filtering of the photooount data is a logical extension of this 
processing technique. Until reoently» howevert the Urbane lidar systea was 
not capable of collecting a sufficient nuaber of profiles in one ni^t to 
peralt teaporal processing. In June, 1979 a new receiving telescope was 
installed and in Januaryt 1981 a new laser was Integrated into the lidar 
systea. A description of these iteas is included in Chapter 2. Systea 
perforaanoe has iaproved draaatioally and« because of the lengthy data sets 
now being oolleoted, teaporal filtering has beooae a viable addition to the 
data processing techniques used to filter sodiia lidar photooount data. 

As in the spatial filtering techniques developed by Rowlett and Gardner, 
out of the aost critical aspects of the teaporal filtering process is the 
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deteralnation of an approprlata outoff fraqutnoy* Aotuallft ainoa tha 
photooount data is now baing two-dlaanalonallj flltarad, tba spatial and 
taaporal f 11 tar outoffs should ba jointly dataralnad. Tha two«dlnanslonal 
flltar daaorlbad in this ohaptar Is a saparabla flltar so that tha tanporal 
and spatial flltar outoffs nay ba salaotad and laplaaantad indapandantly. 
Howevary baoausa of tha dalatarlous affaota of shot nolssy tradaoffs onist ba 
sada In datamlnlng thasa parsaatarsy 1. a.y tha spatial outoff nay ba 
axtended by saoriflolng tanporal raaolutlon and vloa varsa. Thasa outoffs ara 
salaotad by axanining tanporal and spatial parlodograas of tha photooount 
data. Tha pariodograny whioh Is usad ;i» an astlaata of tha powar spaotrun 
the datay is computed using tha Dlaorate Fourier Transfora (DFT). 

3.2 TEMPORAL PERIODOQRAM OP PHOTOOOUNT DATA 

The temporal pariodogram Is usad as an astlaata of tha tanporal powar 
speoti ID of photooount data. Fron this astlaata y tha outoff fraquenoy for a 
temporal filter will be detamlned. Howavary there are several dlffloultl^a 
anoounterad in oonputlng and analyzing this parlodogrsa* Tha tanporal 

variations In the sodium layer density at a fixed altitude ara usually a small 
paroantaga of the average baokground sodliD density. This Is illustrated In 
Figure 3.1 where the signal ^otooounts in tha 80y 65 and 90 kllonatar range 
bins ara plotted as a function of tine for data oollaotad on October 28-29t 
1979. At 85 and 90 kn tha average level Is vary high with rapid shot noise 
fluctuations superimposed on top of snail amplitude y long tarn varlatlms. In 
addition y there are sharp dlsoontinultias In the temporal data at tha 
beginning and end of tha observation period. Tha sharp dlsoontinultias along 
with tha high average level dominate tha tanporal pariodogram and overshadow 
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Fisure 3* 


OCTOBER 28-29. 1979 



TIME (CST) 


T«aporal variations of tlM raw photooount dal& in tha 80, 85* 
and 90 lot range bins COartoer and Shelton, 1981]. 
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spectral ooaponents due to actual teaporal variations in the sodlun layer* 
This problM can be aseliorated by subtracting the average sodiua background 
level f roe each range bin and then aultiplylng the resulting tcaporal data by 
a HlndoH function to eliAinate the discontinuities at the beginning auid end of 
the observation period. A Haamlng window is typically used. The nature of 
the Hanning window is discussed in section 3.6.1. The teaporal periodograa 
can then be conputt.' by taking the Fast Fourier Transfom (FFT) of the 
windowed data and conputlng the square of its aagnitude* 

However f it should be noted that windowing does introduce some spectral 
distortion. To see this let x(.kAh,mtT) denote the photocount in the range 
bin corresponding to the .<n^^ profile. Lh is the altitude separating successive 
range bins and df I 3 the time increment between successive spatial profiles* 
The average photocount level in the range bin due to sodiua and background 
noise is 


“ I x(kAh,mAr) (3. 1) 

where K is the total number of spatial profiles collected and N{kAh) is the 
suo over all spatial profiles of the photooounts in the range bin* If the 
window function is denoted by u(mA7) then the temporal periodogram for the 
range bin is 


HkLhtbi) 


K-l 


N(kLh)^ -imLT 


I w(mtT){x(kth^mt^) - g 

m-0 * 


(3.2) 


To simplify notation, Lh and A? will be dropped and u, x and N will be 
expressed as functions of k and m . 

The expected value of the temporal periodograa may be evaluated by noting 


that X and S are Poisson. This implies that 
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<x(k»m) x(fc,n)> « <x(/c,m)> <x(fc,n)> m 4 n 

<x(k^m) x(JCtn)> ■ <x(JtjW)>^ + <x(,k^m)> m • n 

<N^(k)> - <N(k)>^ + <»(*:)> 

<ff(k) x(fc,m)> ■ <x{ktm)> + <x(fe,m)> <N(k)> 
Taking the expeoted value of equation (3*2) yields the result} 

/f-1 


: I ♦«,«)! 


Jf-l X-1 


X w (m) <x(k,m)> 
n^O 


I I "y u(w) - <x(fc,m)> - <x(fc,n)>] 

m-0 n«0 * * 

Ultimately I this reduces to 

<|*(fc,u)|^> - |<r(fc,«)> * w(bi) - -^^1^ - ^ |fi^(u)r 

«-l 2 

+ J] w (m) <x(fcjW)> 
m»0 

<ra,w)> - J [<xik,m)> ]e~^ 

m»0 ^ 


where 


“turn 


(3.3) 

(3.4) 

(3.5) 

(3.6) 


+ I I Hm) M(n)(<x(fe,m)> - [<x(fe,n)> - (3.7) 

m*0 n»0 * " 


(3.8) 


V(o») ■ J u(m) e 
m*0 

<r(fc,(d)> Is Just the DPT of the expeoted photooount variation In the A:th range 
bln. |<r(X;,<a)>| Is the function we are attempting to estimate by computing 

I ♦(*£,«)! 

As can be seen from equation (3*8), |«|^wlll provide a slightly distorted 
estimate of |<r>| . nie first term In equation (3*8) Is approximately equal 
to |<r>|^. However, the convolution of <r> with W reduces spectral resolution 
and the W/K term Introduces some distortion for m<l/(K6T)» The last term In 
eqiMtlon (3*8) Is due to shot noise and la approximately equal to <N(k)>. At 
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t«aporal fraquenoles this shot noise term Is the aost slgnifloant In 
equation (i.8). To reliably observe teaporal variations in the lidar datai 
|<r>P (i*e.y the first tern in equation (3*7)) nust be large oompared to the 
shot noise level. 'Hie second tern Involving is a distortion tern idiloh 
is significant only at low frequencies where Mote also that the 

periodograa is a biased estiaate of |<r>p. The bias is a constant level 
which depends upon the total number of photons counted in the temporal profile 
and upon the temporal window, u(m)< At wsO, the expectation obtained from 
equation (3*7) is 

<|*(?c,0)|^> « I <x(k,m)> - if ( I w(m))H I w(m) <x(fc,m)>] 

m»0 m»0 m«0 

1^1^ M(m) ^ + Y w(m)(<x(fe,m)> ' ^ (3-9) 

r |m-0 J [m-0 * 

If no window is applied, l.e. if w{m) s 1 for all m, then equation (3*9) 
becomes zero. (x(/c,m) - li(k)/K) is by definition zero mean. However, in the 
isresence of teaporal windowing, b)(m)(x(k,m) - S(k)/K) is no longer zero mean. 

If no teaporal window is applied, the expected value of the periodograa 
becomes 

<h(fe,w)|^> = |<r(fc,w)>|^ + <ff(fe)> 




sinCj) sin(j) 


- 7 Reie 2 ^<r(fc,w)>} 


where XUi^a) is the DFT of it tras assaoed that K teaporal saaples 
were collected at intervals of AT. Except for a> s 0, at frequencies tdiioh are 
integer aultiples of (2ir)/(XAT) the third tens in the above expression is 
zero. Since the FPT computes the values of the DPT only at these points, if 


the PFT is used to compute the periodograa, this last term will contribute 
only for <•> * 0. At this point, its >^lue will be ><iV(fe)>. This effectively 
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oanoals the oonetant shot noise tern peralttlng the antiolpated result that 
(In the absenoe of temporal windowing) <|t(lc 40 )P> sO* 

Unfortunately , the periodogran provides a very noisy estimate of the power 
speotrum [Oppenheim and Schafer t 19771* In addition t It Is often difficult to 
obtain a simple expression for the oovarianoe of a perlodogram. Such Is the 
oase for the oovarianoe of |e(itcja>)p unless it is assumed that the expected 
value of photooounts in the temporal profile, <N> , Is known. This assumption, 
of course, also yields a sllj^tly different expression for the expected value 
of the perlodogram since <N>/K is now subtracted from x^k^m) rather than N/K . 
The new perlodogt'am is 


*TO“0 

and the expected value of this perlodogram becomes 


(3.11) 


<|4-(k,u))|^> - |<r(fe,o))> * W(w)|^ + I w^(m) <x(k,m)> (3.12) 

^ m-0 

As before, the signal spectrum falls into a constemt shot noise level given by 
the second term in the above equation. The distortion at low frequencies Is 
no longer present. The convolutional smearing of the desired power spectrum 
by V(u) remains present, however. 

To obtain the variance of |42(^»‘^)P> the covariance is first evaluated. 

For a photooount perlodogras, the covariance is given by 

/f— 1 K— 1 Jf — 1 1 

cov[|« (fe,u) )|^, |♦_(^:y^i)2)|^] = III I [<Y(fe,j)Y(fe,i)Y(fe,m)Y(fe,n)> 

^ ^ ^ ^ 

- <Y(fejj)Y(fc, ^)><Y(fejn)Y(fe,n)>] e^“l^ ■ ■ ”) (3.13) 


where y(^,j) * w(j) [«(fej J) 

The summation is broken into special cases: kslsn^n, etc. These 
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oases are then evaluated by using the property of Independent Inorenents and 
the required moments of a Poisson prooess [Salehf 1976]. The final result for 
the oovarlanoe is 

^-1 4 2 

« I w^(m) <x(kj/n)> + " “ 2 ^ * " “2^^' 

w*0 

+ 1 ^ 2 ( 0 )^^ + tti^) * <X{k,ui^ + W2)>|^ 

+ 2ffe{[&'(o)^) * <r(fe,Wj^)>ltI/3(«j^) * <XiKuij)>]*} 

+ 2ffe{[.V(o»2) * <r(fe,W2^>ltl/3(a)2) * <X(fe,W2)>l*} 

+ + «2)>Hfi^(«j) * <r(fe,«j^)>]* 

tl/(w 2 ) * <TiKoi^)>]*} + 2Re [ 1 ^ 2 (Uj^ - W 2 > * - “2^^^ 

[Wim^) * <r(fc,u»^)>]* [I/(a. 2 ) * <r(fe,o» 2 )>l> 

K-1 . 

where </(k^o))> » ^ <x(k^m)> e 

WM - 

m-K) 

WAm) = I w^(m)e~^^ 

^ m-0 

The variance is found by setting <i>j^"U 2 *u in equation (3*14): 

, /f -1 . 9 9 

Var(|» 2 (^>“) I ] ■ I w (w) <x(fc,m)> + [ J w (w) <x(Km)>] 

m»0 OT»0 

+ ||/2(2'*>) * <X(.k,2u)>\^ + 4i?e{tl/((tf) * <r(k,u)>] [l/^Cw) * <AT(?C,w)>]*} 

ty 

+ 2i?e{{l/2(2u)) * <f(k,2u»)>Hfc^(«) * <T{k,u)>]* } 

K-l , » 

+ 2 [ I w‘^(m) <«(fe,m)>] ll/(a)) * <r(fe,(u)>r 

m«0 


(3.15) 
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r 2 2 

Aa u becomes largei the varlenoe epproaohea ( ^ w On)<a(k,m)>) * Thus, 

m»0 

In this region the varianoe la equal to the square of the expected value* One 
method of reducing the variance of the perlodogram la to average the 
perlodograms of a number of temporal profiles. This technique (Bartlett's 
prooedure) [Oppenhelm and Sohafert 1977] has been used to reduce the varianoe 
In spatial perlodograms (Rowlett and Gardner, 1979]* If the true power 
speotnmi whloh we are attempting to estimate does not change from temporal 
profile to temporal profile, the resulting average perlodogram will have the 
same expected value, but the variance will be reduced by a factor of *7, where 
J is the number of profiles averaged together* Actually, spectral features 
will vary with altitude as a result of the nature of the layer response to 
wave activity. However, the temporal perlodogram is being computed for the 
purpose of choosing an appropriate temporal filter cutoff frequency* Thus, 
the slight distortion at low frequencies Induced by the temporal window and 
the effects of averaging a nusber of perlodograms are not significant* 

3*3 TEMPORAL PERIOOOGRAM OF SPATIALLY FILTERED PHOTOCOONT DATA 

It Is instructive to detenslne the effect spatial filtering has on the 
temporal resolution of the Ildar photocount data* It Is anticipated that the 
achievable temporal resolution should be Increased In data trtiloh have 
I»*evlously been spatially filtered. The expected value and variance of the 
temporal perlodogram of spatially filtered data may be found by noting the 
fact that, rather than raw photocount data, an estimate of the power spectrum 
of data given by 

{N-D/2 

yikiWi) » \ hik - l)[x{l,m) - b(m)] (3.16) 

7 . 

^ 2 
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Is being sought* N Is the niaber of points in the spatial profile. h{k) is 
the spatial filter iapulse response whioh is oonvolved tiith the photooount 
datSf x(fcjm). b(m) represents the constant background noise present in them 
spatial profile due to sources suoh as PMT dark counts and background sky 
illualnation. This is assuned to be a deteralnistic quantity. In praotioey 
an aoourate estlaate of b(m) can be computed froa the lidar returns obtained 
from above the sodium layer. Thust y(ktm) represents the spatially filtered 
data for the range bin and the m^^ time bin. It is a simple matter to 
determine the expected value and the first four moments of y{k,m) mhioh are 
required to evaluate the covariance of the perlodogram. 

Hk - l)[<xa»m)> - Mm)] (3.17) 



(W-D/2 

<i/(fc,m)> ■ 

’ I 


Z— (i?-l)/2 

<y^{k,m)> 

2 

« <y(ktTn)> 

<t/^(fc,m)> 

- <j/(fc,m)>^ 


(ff-l)/2 


I r(fe - 1) <xU,m)> 
Z»- (ff-D/2 


(3.18) 


(i?-l)/2 

I 

Z— (ff-D/2 


h ik - 1) <x^(Z,m)><x(Z,m)> 


+ I h-’(k - 1) <x(l»m)> 

Z— (ff-l)/2 

(i?-l)/2 2 

<y^(kjm)> - <y(kfm)>^ +6 I nik - Z) <x (Z,m)><x(Z,m)> 
(/?-!) /2 3 

+ 4 \ h (k - Z) <x (Z,m)><x(Z,m)> 

Z— (iV-l)/2 ^ 

(ff-l)/2 2 

+3 ^ n(k - Z) <x(Z,m)>^ 

Z— (ff-l)/2 

(iV-p/2 , 

+ ^ h^(.?c - Z) <x(Z,m)> 

Z— (ff-D/2 

<x^(Z,m)> - <x(l,m)> - bim) 


(3.19) 


(3.20) 


where 
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The teaporel periodospaa of the epetlelly filtered data for the rente 
bln la given hyt 




X-l 


I w(m)tj/(M) -^1 

m-0 


(3.21) 


The expected value of the periodograa can now be oaloulated hy evaluating the 
following expression: 

<1* (k,u)|^> • I f u(m) w(n) <(i/Oc,m) - - 
® m»0 n-0 

. SIM], - «) (3.22) 

As beforsi w(m) is the teaporal window sndS*ik)/K is the average nuaber of 
signal photooounts in the spatially filtered tsaporal profile. After 
aanipulationi an expression which is siallar to equation (3.8) is obtained. 

<!♦ (fe,«)| > * ^ ^ ..2r_\ 1.2 


\ w {m) h (k - r) <x(r,77i)> 
mmo r»-(A?-l)/2 

<r^(fc.«)> * WM - ^1^ - V («»)r 


where <r (k,«)> - I [<y(kttn)> - — e”'*’""' 
® bH) ■ * 


(3.23) 


ir-1 (iV-l)/2 

<H’ik)> • I I h(k - r)t<«(r,m)> - b(m)] 
m-0 r»-(/?-l)/2 

Once again) distortion at low frequenoies is encountered and a constant 
shot noise level is present. Assuaing that spatial filtering does not alter 
the doalnant spectral features in the teaporal periodograa (as should be the 
case if the spatial filter is selected properly), the aain difference between 
the teaporal periodograas of raw and spatially filtered data is in the 
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oonstant abot nolae level. The deoreeae in ahot noiae level will be by a 
factor of G vdiere G la given by: 

I ^ ir{m) h^(k - r) <x(ktm)> 

G - (3.24) 

A“1 2 

w (m) <x(ktin)> 
m»0 

The numerator of the fraction in tbia expreaalon ia aiaply the apatlally 
filtered veraion of the denominator. The iapulae reaponae of thla filter ia 
the aquare of the Impulae reaponae of the apatlal filter. Therefore, aa the 
cutoff frequency of the apatlal filter beoonea lower (l.e., the filtering ia 
more aevere), featurea with longer vertical wavelengtha are removed from the 
filtered data. However, the effect of the filter in equation (3.24) alao 
becomea more aevere and the deoreaae in ahot nolae level apparent In the 
temporal perlodogram beoomea more aignif leant. Since the level la lower, 
featurea which were previoualy maaked may now become apparent. Aa a reault, 
the cutoff for the temporal filter (which la aeleoted by examining the 
temporal perlodogram) can be extended yielding better temporal reaolutlon in 
the filtered data. 

Juat aa the expected value of the temporal perlodogram ia affected by 
apatlal filtering, ao ia the variance. Aa before, it la extremely difficult 
to obtain an expreaalon for the perlodogram variance when the expected value 
of the auB of the apatlally filtered photooounta in the range bin, <y(/c)>, 
ia not known. Aa a reault, it ia aaaumed that thla quantity ia know a priori 
and the covariance and variance expreaaiona are evaluated. Rather than repeat 
much of the development preaented in aeotlon 3.2, only the reaulting varianoe 
ejqpreaaion ia preaented here. The development oloaely parallela the previoua 
work. The varianoe expreaalon ia 
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, X-l (ff-l)/2 . . 

varti*, (fcjw)r] •I I vr(m) h^(k - r) <x(.rtm)> 

mO r— (i7-l)/2 

+ |J?2(2w)[V2(2w) * <2:(fc,2«)>l|^ 

K-1 (/l?-n/2 2 2 2 

+ [ Z I w^(m) - r) <xir,m)>y 

m"0 r"-(iy-l)/2 

+ A/?e{[V(«) * <r^()lc,«)>]tff3(«)(V3(«) * <Jr(fc,u)>)l*} 

+ 2Re{[H2i2o>){W^(2ia) * <X(fc,2u»)>)] (l/(u) * <T^ik,u)>]*^ 

K~l (i?-l)/2 2 2 . .2 

+ 2[ Z Z " w (w) <x(r,m)>] |v(u) * <r (fe,«)>p (3,25) 

m»0 r»-(/l?-l)/2 

iC-1 (i?-l)/2 2 

At high frequenolea, the donlnant tern is (Z Z (k-r)<x(r,m)>) 

mO r— (^^-l)/2 

Thus, as in seotlon 3*2, the variance la the aquare of the expected value for 
large u, Aa in thla prevloua aeotlon, the variance of thia eatiaate la 
reduced by averaging the teaporal perlodograaia fron altitudea between 80 and 


100 ka. 

Figure 3*2 ahowa the average teaporal perlodograa computed froa data 
collected on October 28-29} 1979* Thia data aet oonaiated of 71 apatial 

prof ilea collected at 5 alnute intervale* The 3 dB deoreaae in ahot noiae 
level after apatial filtering ia apparent* In thia oaae} the apatial filter 
waa a deconvolution filter of the type uaed by Rowlett and Gardner with a 
apatial cutoff at a frequency of 0*44 km"^* 

3*4 SPATIAL PBRIODOGRAM OP PHOTOCOOIIT DATA 

Juat aa an average temporal periodogram la uaed to detemlne the 


appropriate temporal filter cutoff} the apatial periodogram ia uaed to aeleot 
a apatial filter cutoff* The nature of the apatial filter ia dlaouaaed in 
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FREQUENCY (min' ') 


Figure 3*2 Temporal periodogram of the sodium Ildar photooount data 

oolleoted during the night of October 28-29 > 1979* The outoff 
frequency for the spatially filtered data Is 0*44 Ian” • [Gardner 
and Shelton I 1981]. 
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section 3.7* The epetiel periodogrMi of pbotooount date hes been enelysed In 
detail by Hewlett and Gardner , [1979]* Their results are presented here for 
comparison with the spatial periodogram of temporally filtered data. 

The spatial periodogram of photooount data for the ^ spatial profile is 
given by 

I S ^ 1 • 2 1 ^ 2 

I - I [*(k,m) - (3.26) 

k -0 ^ ') fe -0 * 


Notioe that the background counts are subtracted from the raw lAiotooounta. It 
is assumed that this background count is known a priori. Using the fact that 
xikttn) is Poisson, the expected value of is readily found. 

2 ff-1 2 

<le(u»,m)| > - I <a(k,m)> + |<e(u),m)>| (3.27) 

k-0 


The first term is a oonstant shot noise level and is equal to the expected 
total number of ikiotons counted in the spatial profile. The seoond term is 
the power speotrma of the expected signal photooounts. This latter term is 
the quantity we are attempting to estimate by computing the periodogram. At 
<i» s 0, the expectation is <P(m)> <P^(m)>^, where P(m) is the total number of 

photooounts and P^(m) is the total number of signal photooounts in the 
spatial profile. As w grows large i the oonstant shot noise term dominates the 
periodogram. The shot noise level is down from the value of the periodogram 
at w 8 0 by a factor of 



<P (m)>^ 
, I - 
<P(m)> 



(3.28) 


As mentimed in lurevlous seotions, if this shot noise level is too high* it 
will obsoure features in the power opeotrum. From equation (3*26) it is clear 
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that inor«Mlng th« mabar of aifoal photom will doopoMO tht shot noise 
level* 

The vsrlsnoe of the spstlsl perlodogpsa hes also been derived by Rowlett 
and Oardner. 

vart|e((D,ffi) 1^3 ■ 2<P(m)> 1<6(«,*:)>|^ *f 2i?e{<X(2w,m)><e*<Wjm)>^} 

+ 4/?e{<2r(u,m)><e*(w,m)>^} + |<2r(2«,m)>|^ + <P(m)>^ + <P(m)> (3.29) 

where <X(u,m)> ■ I <x(kjm)>e 

k*0 

At <D s 0, the first two terns donlnate yielding a varianoe of U<p(m)><P^(in)>^» 
For large u, the varianoe is approxiaately <P0»)>^ * As in previous oaseai it 
is seen that for large o > , the varianoe of the perlodograa is the square of its 
expected value. Perlodograa averaging is used to reduce the varianoe of this 
estlaate of the power speotrua by averaging the periodograas of individual 
spatial profiles together* nils procedure was disoussed in earlier sections 
on teaporal periodograas. 

3*5 SPATIAL PBRIODOQRAM OP TEMPORALLY FILTERED DATA 

Just as spatially filtering the data prior to ooaputing the teaporal 
perlodograa lowers the shot noise flooPt teaporal filtering ovi be used to 
lower the shot noise floor in the spatial perlodograa. This peraits the 

spatial cutoff frequency to be extended by extending the frequency at which 
the true power speotria beooaes obscured by the shot noise level. The 

teaporally filtered data is given by: 

(M-D/2 

yAKm) ■ \ g(,m - i) (3.30) 

® t-(i»^-l)/2 

Vhe perlodograa of the teaporally filtered spatial pro..'ile is 
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^ k-0 * 


(3.31) 


b*(m) is ths bsokground nolss count ppsssnt in the spstlsl profile of ttoe 


teaporelXy filtered date. M is the length of the tenporsl filter iapulse 


response* The ohsreoteristios of the filter ere diaouased in section 3*6.1. 


The expected value end verisnoe of \ ^ wmy be found in e aenner 


siailsr to that presented in seotion 3*3 ehere the teaporel periodograa of 


spetislly filtered date ves exaained* This problsa is analogous to the 


present problsa of evaluating the expeoted value and variance of the spatial 


periodograa of teaporally filtered data. In the present problsa y however* no 


window is applied to the data before the periodograa is o<»puted* The 


e^qMoted value is 


<|0j,(u,m)|' 


N-l (f1-l)/2 , » 

I I g^im ~ i) <x{M.,i)> + I<e„(u*m)>| (3,32) 

fe -0 i — (/<-!) /2 ^ 


where |<e^<Dym)>|^ is the true power speotrua we are attsapting to find* The 


constant shot noise floor Is due to the presenoe of the first tera. At u « 0 


the expeoted value of the periodograa is 


5 («-l)/2 , 

>■ \ I g (m - i) <x(fc,i)> + 

fc.O i— (W-D/2 


1?“1 (A^l) 

I I 

kmQ 


^l)/2 


g(m - t) <x{k,i)> 


The shot noise floor is below this value by a factor of 


/? - 1 + 


p?-l (M-l)^ 

I I 

Uc »0 t »-( Af»l 


( W - p /2 

I <x(kti)> - b'(m) 

nn (^i)/2 ; " 

I I g^im^) <xikti)> 

[fc-0 t— (W-D/2 


(3.34) 


i 
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The value of toiporal filtering In deoreaslng the shot noise level is also 

illustrated by this equation* As the teaporal filter cutoff is extended t the 

value of equation (3*34) generally beooses smaller until it beoones equivalent 

to equation (3*26) when no temporal filter is used* As the temporal filter 

cutoff becomes low (l*e* severe filtering) « the value of this expression 

beoones small* nils principle may be demonstrated by asavning that the 

teaporal filter impulse response is a rectangular function. In order to 

insure that the magnitude of the filter frequency response la unity at u s 0| 

the value of g(jn) must be \/Mt trtiere M la the length of the FIR Impulse 

response* As the filter cutoff becomes smaller, the length of the filter 

impulse response becomes larger* It is also assumed that xOCftn) - x(k) does 

not vary withm. This, of course, is an unrealistic £issunption. However, it 

simplifies equation (3*34) and serves to clarify the effect that varying the 

temporal filter cutoff has on this equation* Under these conditions, equation 

(3*34) simplifies to 

^-1 -1 

R ~ (1 + Z M<x(k)>) (3.35) 

As the filter cutoff is reduced, M grows larger and R decreases. This 
indicates that the separation between the spectral peak at o> = 0 and the shot 
noise floor is increased* Figure 3*3 contains plots of average spatial 
periodograms for raw and temporally filtered data. The separation between the 
spectral peak and the shot noise level is 11 dB greater after a temporal 
filter was applied. In this case the temporal filter cutoff Is 0*03 min~^* In 
addition, the spatial periodograms were averaged for the entire data set to 
reduce the variance* 

As in previous oases, the perlodogram is a noisy estimate of the power 
spectrum* Using the same methodology reported in earlier sections, the 
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variance expresalon was found to be 


, ff-1 (j«M)/2 . 

Varl|e_(w,m)| 1 •I I g {m ~ i) <x(.k,i)> + |<i?,(2u),w)>| ' 

^ fe-0 ^ 


N-l (M-U/2 , 

I I 0 (m - i) <xikti)> 

k-0 i— (W-D/2 


+ 4i?e{<0y(a),m)><i?2 (u,m)>> 


* 2 

+ 2i?e(<i?2(2wjm)><0y (a),w)> } 


+ 2 


ff-1 (W-D/2 , 

Y Y g^im - i) <xikti)> 
k~0 i=(W-l)/2 


|<0j,(u,m)>l' 


(3.36) 


N-l («-l)/2 ^ 

where <f? (oj,m)> * ^ g (m - i) <ac(k,t)> 

” i=^(W-l)/2 


The form of this expression la quite similar to the variance expression for 
the spatial periodogram of raw photocount data (equation (3*29) )• As before, 
at hi^ frequencies the variance la the square of the mean. To reduce the 
variance, the temporally filtered spatial perlodograms are averaged. 


3.6 TWO-DIMENSIONAL FILTERING 
3.6.1 TEMPORAL FILTER 

The perlodograms discussed In previous sections are used to determine the 
cutoff frequencies for t«nporal and spatial filters. First the average 
temporal periodogram Is computed and examined. A teaporal filter with the 
desired cutoff frequency Is designed and applied to the data. At this point, 
the spatial filtering technique developed by Rowlett and Gardner Is used to 
complete the two-dimensional filtering process. The implementation of the 
temporal and spatial filters are quite different. Aspects of the temporal 
filter are discussed in this part of the section and the spatial filtering 
technique Is reviewed at the end of this section. 
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technique is reviewed at the end of thla section. 

Because the sodium layer has infinite temporal extent, as opposed to the 
limited spatial extent of the layer, temporal windowing wais utilized in 
computing the tonporal periodogram. The limited temporal extent of the data 
set presents problems in the filter implementation as well. To avoid 

distortion Induced by the sharp discontinuities at the beginning and end of 
the observation period, a Finite Impulse Response (FIR) filter is used. This 
has an added advantage in that it is a simple matter to insure that the filter 
has a linear phase response. 

Since the returns from a number of pulses are integrated to yield a 
spatial profile of data, it is possible to apply a temporal deconvolution 
filter to reduce the distortion introduced by this effect. This distortion is 
much the same as the distortion Introduced by the range bln width in the 
altitude profile. If the returns are integrated for minutes (by summing 
the returns from all the laser pulses in this interval), the returns are 
effectively smeared in the temporal domain. This may be illustrated by noting 
that the photocounts for each laser pulse, ^p(^.'i), are actually simimed over a 

period T to obtain the integratec photocount x(/c,m). 

P 

M -I 

x(feA)i,mAr) ^ I X {kLhtmtT + nAx) (3.37) 

n«0 P 

where is the number of laser pulses per spatial profile. More precisely 
T s /d^Ai Where Ax is the period separating laser pulses. This is actually the 
convolution of x(7c,n) with a rectangular filter Impulse response. The 
convolution is equivalent to a multiplication in the temporal frequency domain 
by ain(ir/7’p)/(7r/5’p) where f is the temporal frequency. The first zero occurs 
in this function at / s ^/T^t where is usually on the order of 2 minutes. 
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Examination of temporal perlodograma Indloates that there Is no disoernible 
structure whloh oan be attributed to temporal variations In sodium density 
with periods less than 20 to 40 minutes. Thusy the distortion introduced by 
this convolution Is negligible in current data and there is little need for a 
temporal deconvolution filter. 

A linear phase low pass filter will be utilized to reduce the shot noise 
oontrlbutlons in the temporal data. A filter with short Impulse response is 
desirable in order to insure that the length of the temporally filtered data 
set is comparable to the length of the raw data set. The nunber of points in 
the temporally filtered data set (A^) is less than the number of points in the 
original data set (K), The relation between and K is 

^ K - M I (3.38) 

where M is the duration of the filter impulse response. To maintain linear 
phase y the Impulse response must satisfy the relation 

g(m) ■ g(M - 1 - m) (3.39) 

This amounts to a requirement that the Impulse response be symmetric about its 
midpoint. If this condition is satisfied , the filter is linear phase and the 
various spectral components in the temporal response of the layer will all 
suffer equal delays. This can be an important consideration when trying to 
Interpret filtered data in terms of a layer density response to wave activity. 

In addition to the requirement for a short Impulse response, we desire to 
have small sidelobes and a sharp transition from the pasaband to the stopband 
in the filter frequency response. Unfortunately, these requirements are at 
odds with each other since a short impulse response usually means large 
sidelobes or a gradual pass band to stopband transition. A compromise is made 
by selecting a Hamming function as the Impulse response. The Hamming 
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fwotlon, which la a meaber of the family of ralaad ooalne ourveat'la 
aa 


gW) 


'0.54 - 0.46 coa(^); 0 <m <M - 1 

k 0 ; elsewhere 


defined 

(3.40) 


The magnitude aquared of the DPT of thla impulae reaponaey |{7(2i;)| « haa a 
atopband which la In exceaa of 40 dB below the paasband. The width of the 
paaaband la 2/(Af47), where AT la the interval between apatlal prof Ilea. The 
paaaband la defined aa the range of frequenolea between u s 0 and the first 
null In the DPT of the filter frequency response. Figure 3.4 contains a plot 
of the magnitude of the filter frequency response for a 13 point Haaoilng 
Impulse r<38ponse when AT = 5 minutes. In thla case, the cutoff frequency la 
0.0308 mln~^ (corresponding to a period of 32.5 minutes). It la evident that 
the response In the passband Is not flat. Thus, the higher frequency 
components In the passband will suffer substantial attenuation as the cutoff 
frequency Is approached. By selecting cutoff frequencies which are somewhat 
greater than the frequency at which the shot noise level dominates the signal 
portion of the temporal periodogram, this problem Is circumvented to some 
extent. Unfortunately, the requirements for a flat passband, linear phase and 
short Impulse response cannot all be meet simultaneously. The Hamming Impulse 
response filter la a reaisonable compromise. 

The effect of temporal filtering on the expected value and variance of the 
photocount data Is easily determined. Since x(7c,m) Is Poisson, the 
photocounta collected In non-overlapping time or range Intervals are 
statistically Independent. The expected value and variance of filtered data 
corresponding to the range bln and the spatial profile, yiktm), are 
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(«-l)/2 

<y>r^Km)> - I g(m - i) <x(kti)> 
i— (W-D/2 

(W-D/2 , 

Var(j/(fc,m)] - I g (m - i) <x{kti)> 

i— (M-D/2 

The signal to noise ratio (SNR) at the filter output | defined as 

2 


SNR 


is seen to be 


Var[t/j,(fc,w)] 


L 


(W-D/2 

I gOn - i) <x(kji)> 
SNR . 


-i2 


(W-D/2 , 

I g^(m - i) <xik,i)> 
t— (Af-l)/2 

The next step in the processing teohnlque is spatial filtering* 


(3.41) 


(3.42) 


(3.43) 


(3.44) 


3*6.2 SPATIAL FILTER 

The spatial deconvolution filter which is applied after temporal filtering 
was developed by Rowlett and Gardner [1979). In equation (2.7) the Ildar 
equation was used to express a spatial profile of data as the three-fold 
convolution of the radiative lifetime of the sodlimiy l(s)| the receiver range 
gatSf 9(^)f the laser pulse shape, p(s), and the sodium scattering profile, 
Y(s). The application of temporal filtering does not affect the nature of the 
spatial deconvolution filter. It C€m be seen by examining equation (2.8) that 
the proper deconvolution filter is simply 


B(u) 


0 |u| > (1) 


(3.45) 


where is the cutoff frequency determined by examination of the spatial 
perlodogram of temporally filtered data and H(u) is the DPT of the spatial 
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filter Impulse response t h(k)» The baokground count is subtracted from eaoh 
spatial profile and the resulting profile Is transformed via FFT teohnlques. 
At this point) the spatial deoonvolutlon filter) A(u)t is applied to the 
temporally filtered data by multiplication In the spatial frequency domain. 
After an Inverse Fast Fourier Transform) the spatially and temporally filtered 
data are range squared corrected (l.e.) the filtered data are multiplied by 
the square of their range). By applying the calibration technique discussed 
In Chapter 2) the final estimates cf sodium density can now be calculated. 

The expected value and variance of the two-dlmenslonally filtered data are 
easily found. 

iN-\)/2 

<z(k,m)> • I Hk - j)[<yAd,m)> - b(m)] (3.46) 

j— (ff-l)/2 ' 

(N-D/2 (M-D/2 - - 

Var[a(k,m)l •‘I I h‘‘{k - j) g^(m - t)(<x(j)i)>-fc(m)) (3.47) 

(/7-D/2 i— (M-D/2 

where (m ) Is the baokground count present In the spatial profile and 
ziktm) Is the two»dlmenslonally filtered data. It Is assumed that b(m) is 
known a priori. The signal to noise ratio may be expressed) as In equation 
(3.**3), 


SNR » 

" (A?-l)/2 (M-D/2 (N-D/2 

I I h(k-J)g(m-i)<x(j,i)^~ I h(k-c)b(m) 

J— (H-D/2 i— (M-D/2 .i—(H-D/2 

(N-D/2 (M-D/2 , , 

I I hUk-j) g^(m-i) <x(j,i)> 

Jm.(n-D/ 2 i— (M-D/2 


(3.48) 


As mentioned previously) this procedure ultimately yields a 
two-dlmenslonally filtered data set for which the two-dimensional filter Is 
separable. The frequency response of this filter is given by (7 (w^) A (u^), where 
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Is the tMsporal frequency and Is the spatial frequency. The paraaeters 
and iapleaentatlcn cf each filter are independent* The order in which the 
filters are applied has nc effect on the filtered data. Temporal filtering 
has been applied first simply because a more dramatic decrease in shot noise 
level is apparent in the spatial periodogram than is evident in the temporal 
perlodogram after spatial filtering* The temporal and spatial cutoff 
frequencies may also be chosen by examining the two-dimensional periodocram of 
the i^otcoount data. The charaoteristios of this periodogram are discussed in 
the next section. 

3.7 THO-DIHENSXONAL PERIODOGRAM 

The two-dimensional periodogram discussed in this section can be used to 
choose the optimal set of cutoff frequencies for the spatial and temporal 
filters discussed previously. The constant background count, is subtracted 
from each spa tied photocount profile. It is assuaed that this quantity is 
known and, to simplify the following analysis, is constant from spatial 
profile to spatial profile. Once again, the nature of the temporal variations 
in sodiiD density present a problem. These variations are a small percentage 
of the background sodium density. Unlike the temporal periodogram, however, 
the average photooounts at each altitude cannot be subtracted since the 
spatial variations in the data set are also of interest. A temporal window is 
applied to prevent the large discontinuities at beginning and end of the data 
set from contributing significantly to the periodogram. The periodogram is 


expressed as: 
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|r(w^,wj,) r 

N-l iV-1 K~l K-l . . tj, , /_ 

■ I I I I tfW »(n) X (M) 

kmo Z-0 m-0 n-0 ’ ’ 

Ug and (d^ ara apatlal and tanporal frequenolaa, raapaotivalyf and z^(m,n) s 
xim^n) - b » w(m) Is the temporal window. The expected value of this 


expression Is evaluated yielding the result! 


, N-l K-l , 2 

<|3’(u-,a)„)r> - I I W (m) <a;(/c,m)> + |»/(wj * <r(uj-,(«)„)>| (3.50) 

* fc»0 m-0 ^ ^ ^ 

2 

where W(u ) Is the DFT of u(m) and |<r(b)eM>>m)>| two-dimensional power 

I ST 

spectrum we are attempting to estimate. As in previous perlodogramS} a 
constant shot noise level is present. 

The variance expression for the periodogram is now evaluated. This is 

accomplished by evaluating the covariance t cov[|T'((i)^^,a)^^ 

for s (0^^ and covariance expression for a two-dimensional 

photooount periodogram is given by: 

, , It-l IV-l ff-l If-l k-l K-l K-l K-l 

cov(|r(u„pU„.)|% - llllllll w{n^w{n)wi8)w(t) 

lc-0 1-0 p-0 r-0 m-0 n-0 s-O t-0 

• [<®^ (fe, m) ( It n) (p, a) x^ (r, t)> - m) x^ ( 1, n><x^ (p, a) x^ (r, t> >] 


iui^^(p-r) iuj,j(m-n) (3.51) 

The eightfold summation is broken into special oases: kslansntpsrssstt 

k^ttn^dpsTsahtf etc. in order to apply the property of independent 

increments and the first four moments of a Poisson process. 

<x^(fc,m) x^(l,n)> - <x^(fc,m)><x^(l,n)> k 4 I m 4 n (3.52) 


<x^{ktm)> - <x^(fc,m)>^ + <x^(fe,m)> 


(3.53) 


3 3 

<x^ (k,m)> - <x^(lc,m)> + 3<x^(l:,m)><x(k,m)> + <x(fe,m)> 


(3.54) 


(3.55) 


<x^ (Km)> - <«^(fcjm)> + 6<x^(k,m)> <x(k,m)> + i*<x^ik,m)><xik,m)> 

2 

-¥ i<x(k,m)> + <x(k,m)> 

After tedious algebra! o manipulations the oovarianoe expression is found to be 

2 2 4 

cov(|r(w^j,«-,j)| , |r(u^2»“3'2^l 1 “ I I w (m) <x(k,m)> 

k^O imO 

2 

2 

+ I<5<«5i"“ 52’ “7’r“3’2^’' * ^2Wr“i'2^^ 

+ 2i?e{ [<3’(u)^j, a)yj)> * Ii^(wj,j)H<5(w^j, Wj,j)> * } 

+ 2/?e{[<J'(aj^2* “r2^^ * ^"^'^^“52* “1*2^^ * 

+ 2^e{ [<r(u^j, W2,j)> * J/(wyj)][<3’(w^y Wj,2>> * 

[<S(u^j+w^2»“j'i'*^j- 2^^ * *^2Wl"*™3’2^^^ + 2i?e{ [<r(u)^j, Uj,j)> * fc^(o)2,j)l 
[<r((»)^2* <*>2*2^^ * l^^^***5l"’***52* ***n~***2'2^^ * J^2^***2'l~**’j>2^^ ^ (3.56) 

—ibiffk -'Z-W2^ 

idiere <5(a)_ u_)> ■ / A w(m) <x(k,m)>e e 

^ kmO ffl.0 

<>^((«>m) ■ I w(m)e 

m*0 

fc'’ (uj ■ I w’*(nj)e 

” mO 

Setting u„, s s and w„, * * w_ yields the varlanoe expression: 
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5 It-l A-1 , , 

Var(|2'(«^,Wj,)ri • 2 I I w^m) <x(k»m)> |<2’(u>^,Uj,)> * J/(wj,)r 

k^O w*0 

+ 2i?e{[<T(w^,«j,)> * J/(aij,)]^[<5(2a)^ 2wj,)> * ^ 2 ( 20 )^)]*) 

+ 4i?e{t<!T(fc)^.<.)j,)> * V(Wj)][<5(u^,u) 2,)> * f/3(a.2,)]*> 

, ff-1 X-1 . 

+ |<r(2u<,,2u„,)> * V,(2u_)| + I I W (w) <x{k,m)> 

^ ^ ^ fc-0 wO 

lV-1 /f-1 , 1, 

*11 <x(fe,n?)> ^ (3.57) 

fc“0 m»0 



When u- and u_ equal zero, the first two terms dominate and the variance 

U 1 

is approximatelyt 
Var[|3’(0,0)l^) 

7v-i x-1 , "irff-iif-i 1- 

»A J I w^im) <xik,m)> I I w(m)i<x(k,m)>-b) (3.58) 

fe-O w-0 _ jc-0 m-0 

Aa and increase, the second term falls off rapidly because of its 2u^, 
dependence. For large a>, the last term dcwninates. Since the variance 
does not approach zero as and grow large, this is not a consistent 
estimate of the two-dimensional power spectrum. As was typical of the 
periodograms analyzed previously, the variance is the square of the expected 
value for large frequencies. Following Oppenhelm and Schafer [1977], a 
possible method of reducing the variance of the periodogram is to apply a 
spatial window, u(fc), before computing the periodogram. This has the effect 
of smoothing the periodogram and thus reduolng its variance. Note that the 
same effect was observed when the temporally windowed temporal periodogram was 
analyzed. Applying the spatial window has the effect of convolving the 
frequency response of the window, Uiu„)$ with the spectral content of the 
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signal. The varlanoe expression beooaes sodlfled by this windowing so that 
the term whloh dominates at large frequencies Is 

|2 


N-l K-l , , 

I I i/im) uik) <x(k,m)> 
fc»0 m"0 


and at and « 0, the varlanoe Is 

Var[|r(0,0)|^] 

Ta^-I K-l 


(3.59) 


“ 4 


I J w^(m) u(k) <x(k»m)> 
k^O m“0 


N-l K-l 

I I w(m) u(k)i<x(,k,m)> - b) 
k“0 m*0 


(3.60) 


The expected value of the perlodogram becomes 

, X-l , , , 

<|3’((D^«^)r> • J I w^(m) u^(k) <x(ktm)>+ *•> U(Uf,) W(u^)\ (3.61) 


'5*"r 


m»0 

where the second term is the two-dimensional oonvolutlon of the spectrum of 
signal photooounts and the two-dimensional frequency response of the window, 

i/(w^) V ( Wy ). 

We see that when u and w- « 0 the ratio of the variance to the expected 

ST 

value is 


N-l K-l 


Var ^ U-0 fffO 

• 8 — ^ 


, , "IpV-1 AT-1 

I I w^(m) uik) <xikttn)> \ \ 

J Jc -0 m »0 

pt-l K-l 

*>”)> '*'11 
B1«0 


w(m) uik)i<x(ktm)> - b) 


< > N-l K-l , , 

I I w^(m) u ik) <x(k 
k^O m*0 


w(m) uik)i<xiktm)> - b) 


(3.62) 


and for large frequencies this ratio Is 

n?-i K-l 


Var 
< > 


I I w^(m) u(k) <x(k 
k^O OT-0 


,w)>j 


(3.63) 


Oppenhelm and Schafer [1977] have estimated the varlanoe reduction ratio, 
whloh Is the ratio of the variance after windowing to the varlanoe before 
windowing. For raised cosine curves which are eiq>ressed as 
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u(k) - e - aco*(^) (3.64) 
they estlMte the reduotlon in vmrlenoe to be a faoccr c.' 

R « 2a^ + (3.65) 
Nhere N « length of the window and 0 * 0.54 and a s 0.46 for a Hanning window. 
Under theae oonditionsy R « 0.3974. However) it ia clear fron the varianoe 
expression derived here that the varianoe at high frequenoies is still the 
square of the expected value. Unfortunately) this is the region trttere the 
power apeotruB ia obscured by the shot noise level. Further) the reduction in 
varianoe is fairly aeager and is obtained by saorifioing speotral resolution. 
As a result) spatial windowing does not ueen to offer suoh inproveaent. 

The two-dlnensional power spectra should yield sore insight into the 
nature of sodlia layer dynaslos as well as provide a acre preoise neans of 
deteralning appropriate filter cutoffs* Unfortimately the varianoe of the 
perlodograa oannot be reduced substantially without saorifioing the speotral 
resolution of the perlodograa. Further) the teaporal variations in 
photooounts beooae aasked by the large average levels upon whioh they are 
superiaposed. This (Toblea was solved by subtraoting the average counts froa 
each altitude in the case of the teaporal perlodograa. This teohnique oannot 
be applied here without destroying the spatial variations in photooounts whioh 
are also of interest* ThuS) little oan be done to oorreot these sbortooalngs. 
The draaatio decreases in varianoe possible by averaging the perlodograas is 
not possible because we have only one saaple of the two«diaensional 
perlodograa. 
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4. THE LAYER RESPONSE TO GRAVITY HAVES 


4.1 INTRODUCTION 

Lldar measurenents of the nighttime aodlm layer often reveal wavelike 
density perturbations moving through the layer. It appears that the motion of 
these features Is assoolated with the propagation of Internal gravity waves 
[Rowlett et al., 1978; Clemesha et al., 1978; Negie and Blamont, 1977]. Thus, 
In order to Interpret sodlm layer dynamics the response of an atmospheric 
layer to wave activity must be considered. 

In this chapter two mechanisms by trtilch atmos;dierlc waves may perturb the 
atomic sodlun density are considered. The first mechsuilsm, the corkscrmr 
mechanism^ Is considered only briefly because It is most effective only In the 
upper half of the sodium layer. Further y It Is difficult to account for the 
large density variations observed in the sodlun layer, particularly In the 
lower lutif of the layer, by this process alone, ^e layer density response to 
atmospheric waves, which is derived In section 4.4, appears to provide a more 
satisfactory explanation for the observed sodium layer dynamics. This 
mechanism Is discussed in some detail. It Is found that the layer response 
can be highly nonlinear in certain areas of the sodiuD layer. 

Although it appears that gravity waves are responsible for many of the 
observed sodlun layer dynamics, the results presented In this chapter are more 
general in nature and may be applied to other types of atmospheric waves. 
Gravity waves, however, are an Important type of wave motion at mesospheric 
helots, probably dominating atmospheric motions [Houghton, 1977]. Since some 
of the results la this chapter are discussed In terms of gravity wave 
parameters, a brief re lew of Internal gravity waves Is Included In the next 


section. 
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4.2 INTERNAL GRAVITY HAVES 

Internal gravity waves are a type of ataospherlo wave trtiioh can exist in a 
stably stratified fluid. Since the aedlua is stably stratified > a snail fluid 
parcel which is displaced vertically from its equilibrium level %illl be 
subjected to a force which tends to return it to its equilibrium position. 
Once an atmospheric parcel is displaced and released, a wave may be generated. 
Internal gravity waves are a very important subset of the general class of 
gravity waves. Internal gravity waves exhibit phase variations with height 
tdiile the remaining type of gravity wave, sometimes called an external gravity 
wave, does not. Thus, only internal gravity waves propagate vertically. Many 
of the gravity waves observed in the region of the sodim layer are probably 
generated much lower in the atmosphere and propagate upward to mesospheric 
heights . 

The linear theory of gravity waves is well developed and has been 
considered in detail by Hines [I960]. Rather than repeat the mathematical 
development describing these waves, results which are applied in following 
sections are summarized. The assumptions made in obtaining these results are 
as follows; 

1. Perturbations are sufficiently small to permit nonlinear effects to 
be ignored. 

2. The atmosphere la stationary in the absence of wave activity (no mean 
winds are considered). 
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3* Only forces due to pressure » gravity » and inertia are oonsldered. It 
is assumed that the gravitational field is constant in magnitude and 
direction. 


The linearized equations of motion (4.1), adiabatic state (4.2), and 
continuous mass conservation (4.3), as given below, are solved to obtain the 
perturbations in pressure and density and the horizontal and vertical 
velocities associated with the propagation of Internal gravity waves. 

®o(H) “ 

|| + V . VPq - [II + I- VPq] 

Sp (^-3) 

f + Z-^Pq + pQ^*i^ • 0 

From the solution to these equations it is found that the wave numbers and 
wave frequency are related by the dispersion equation; 


0)^ - + (Y - 1) + iygi^K^ = 0 

X Z X z 

where u = wave frequency (s~^) 

<? s speed of sound (m/s) 

K s k s horizontal wave number (m**^) 

K s + vl2Hs vertical wave number (m“^) 

z z 

Y s ratio of specific heats 

2 

g 3 acceleration of gravity (m/s ) 

2 

p 3 pressure (N/m ) 

3 

p 3 density (g/m ) 

V_ = velocity (m/s) 

H - atmospheric scale height (m) 
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The solution has the fora} 


1 

o 

^ 


- 

1^0^ J 




V 

z 


4exp[i(wt - K^x - 


(4.5) 


The unperturbed density end pressure are given by and respeotively. 

Beoause of the ooaplex nature of K t the aaplitudes of the density and 

B 

pressure perturbations and the vertical and horizontal winds are seen to grow 
exponentially with altitude. The polarization factors » Bt P$ X$ and z» 
deteraine the relative aagnltudes and phases of these quantities. It aay be 
shown that they have the fora 


P - ymh - iygui^/C^ 

B 


(4.6) 


B - + i(Y “ - iygi^^K^ 


(4.7) 


i 


X - oiKKC^ - iguK (4.8) 

Z A 

2 - - (akJc^ (4.9) 

X 

hn exaaination of equation (4.4) reveals that for any pair of real wave 
nuabers 0^ and k )* two distinct » real values of u exist. One of these 

X z 

values aust be greater than - ygllC (the natural frequency) and the other 

l/2 

aust be saaller thanu. • (y-l) gtC (the Brunt-Vaisala frequency), u is 

D a 

always greater than If u > the wave is an acoustic wave. Likewise » 
if <i) < the wave is an internal gravity wave. No internal waves will 
propagate in the frequency gap < u < ot^. 

In this chapter t only Internal gravity waves are considered. These 
ataospherlo oscillations will always have a period greater than the 
Brunt-Vaisala period) At 90 kSf Is approxiaately 4 alnutes. In 


addition) the processes of reflection and dissipation serve to restrict the 
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possible values of k » and w. Beoause of the struoture of the atnosphere 
below 90 km, regions where ^ ^ minutes exist. Speoifioally, near 5M km the 
Brunt-Vaisalla period is probably closer to 8 minutes. This will have the 
effect of preventing some of the higher frequency waves from reaching the 
sodivB layer. In contrast to this effect is viscous dissipation which beoomes 
stronger with increasing altitude. This effectively places an upper limit on 
the propagation of these waves which is a function of the frequency and wave 
mnber of the wave. Figure M.1 [Hines, I960] shows the allowed propagation 
modes at meteor radar heights. 

Although this figure indicates that waves with periods less ttum 10 
minutes c 2 ui be present, only longer period perturbations in sodium density 
have been observed. If a low frequency approximation is made, the 
polarization relations become much simpler. On the basis of experimental 
data, this simplification appears to be Justified. If u « a>, , the relation 

D 

between the horizontal and vertical winds and the density perturbations 
become; 
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X 

- »o] 


- a {y - 1) 
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^0 1 

(A) 

iJ- + g(y - 1)K^ 


(4.10) 
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P - Pq 
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-^u)(l - iyHK) 


. 2 


P - P/ 
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(4.11) 


The winds are expressed in terms of density perturbations to simplify the 
formulation of the layer density response to these waves in a subsequent 
section. 

It is Interesting to note that under the conditions 


k^ » i//c^ 
z a 


(4.12) 


and 


W « U), 


( 4 . 13 ) 



Figure 4.1 Gravity wave propagation modes at oesoapheric heights. Gravity 
wave periods t measured in minutes » are shown in boxes on the 
corresponding constant period contours (solid lines). The 
limits of the permitted spectrum as determined by viscous 
damping are shown for heights of 60, 70, 80, 90, 100, and 110 tan 
(dashed lines). Modes lying above and to the right of these 
curves are excluded. The modes subject to reflection at heights 
of 54 and 79 tai are also shown (dotted lines). Nodes lying 
below these curves cannot proceed from the lower atmosphere to 
the upper CHines, I960]. 
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the dispersion relation slaplifies to 


2^2 . 2. 2 
(rife - <ri, fe_ 

2 D X 

and the polarization relations oan be reduced to 


Z 

X 


-fe 


R . i(Y - 1) 
X~ o 


1/2 


(4. 14) 


(4.15) 

(4.16) 


In this oassf energy flow is almost horizontal trtille the phase progression is 
nearly vertioal* As in all intez*nal gravity waves » the vertical oomponent of 
phase progression will be opposite to the vertioal oomponent of energy flow. 
Thus, an upward transfer of energy aooompanies a downward phase progression. 
Since many of the gravity waves at mesospheric heights probably originate 
lower in the atmosphere and transfer energy upward, they should be associated 
with a downward phase progression. 

4.3 CONVERGENCE OF IONIZATION 

The motion of long>llved ions, such as Na*** in the nighttime mesosphere, 
can be very sensitive to neutral winds. A magnetoionio process exists through 
which horizontal winds blowing across the geomagnetic lines of force Impart a 
vertioal velocity to ionization. As ions are swept along by the neutral wind, 
a Lorentz force will be established due to interactions with the magnetic 
field. Ions will aoquiru a velocity perpendicular to their previous velocity 
and the direction of the magnetic field. Above 145 km, the ion-neutral 
collision frequency is quite low (much less than the ion-gyrofrequenoy) . In 
this case, the ions are effectively trapped on a line of magnetic force and 
will exhibit little vertioal motion due to horizontal neutral winds. At lower 
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heights t below 90 kmt the lon>neutral oollision frequenoy is much greater than 
the lon-gyrofrequenoy* In this oasei the Ions will be swept along by the 
neutral wind. However » In the region between approxiaately 90 and 145 kOf a 
iMt vertical motion will be imparted to the Ionization due to the Lorentz 
force. 

A wind blowing towards the east across the geomagnetic field will produce 
an upward welling of positive ionization in this region. Similarly , a wind 
blowing towards the west will cause a downward movement of lonlzaton. If the 
winds are produced by a gravity wavey for example » a region where the wind 
blows towards the esist will be bounded above and below by winds blowing 
towards the west. At the interface where westward winds are above eastward 
winds, ionization will tend to converge. Where an eastward wind is above a 
westward wind, ionization will be depleted. 

This mechanism has been used to explain sporadic E [Chimonas and Axford, 
1968; Macleod, 1966; Axford, 19633. Gravity waves which propagate upward will 
exhibit a downward phase progression. Thus, layers of ionization trapped at 
the wind nodes will move downward and be dumped around 90 km irtjere this 
meehanjam becomes ineffective. Figure 4.2 [Chimonas, 19693 contains the 
result of a computer simulation indicating how ionization can be collected and 
transported downward by a gravity wave. Richter and Sechrlst [19783 have 
calculated the Na* ion drift velocity (Figure 4.3). Below 90 km these drift 
velocities are quite small. However, this mechanism can cause layers of 
sodium ions to be swept downwards and dumped in this region. When these ions 
are converted to atomic sodium, corresponding enhancements in Ha density would 
be apparent in Ildar observations of the sodium layer. Since the 
[Na3/[Na'^3 ratio is high (probably 30 to 90 [Richter and Seohrist, 19783), it 
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Computer simulation of the domward motion of Ionization under 
the Influence of gravity wave motions. Initially the Ions were 
equally spaced In altitude (at T s 0). The downward motion and 
accvnulatlon of the Ions Is Indicated by the solid lines. 

Dashed lines Indicate the motion of the wind nodes. A maximum 
horizontal wind of 50 m/s and a vertical phase velocity of 1 m/s 
were assuned. [Chimonas. 1979]* 
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appears that this aeohanisa is not responsible for the large sodlua density 
variations observed In the nlghttlae layer. The aeohanlsa Is not very 
effeotlve below 90 ka where the largest fluotuatlons In density are found, 
niese large density variations are aore readily explained In teras of 
perturbations In the neutral layer associated with the passage of ataospherlo 
waves. 

4. it LkIBR DENSITY RESPONSE 

4.4.1 RESPONSE TO PERTURBATIONS; THE CONTINUITY EQUATION 

Althou^ we are prlaarily oonoerned with the aesospherlc sodium layer, the 
following development Is less restricted in nature and applies. In general, to 
atmospherlo layers composed of a neutral alnor constituent. The density 
response of a neutral atmospherlo layer to wave activity may be found by 
examining the equation of continuity; 

V(n£) - P - e (4.17) 

where n • density of the alnor constituent 
V s velocity of the alnor constituent 
P s source teras 
Q 8 loss teras 

The velocity field, V , associated with an ataos;^erlo 
general, be written as follows; 

V ■ exp(i(ii»t - X»r)) 

where A s rave amplitude 

e 8 fi JC B J6 

— X B 


wave can. In 
(4.18) 


<i> 


8 rave frequency 
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K»Kt*K&* wave v«otor 

— « B 

r • xt * at • position vsotor 

V is sxprsssed as a ooaplsx sxponsntial in ordsp to simplify ths problsa. 
Beoause equation (t.17) is nonllneart only tbs real part of V should be used* 
However I for the first order approxiaatlon aade in section 4.4.2 little error 
is introduced. In the perturbation series solution presented in section 4.4.3 
the baok-ooupling of high frequency perturbations to lower frequency 
perturbations is neglected beoause of this approximation. The osoillatory 
density variations associated with the second order terms are accurately 
represented. Higher order terms suffer larger errors. However < the 
oontributions of the higher order terms to the total layer response are not 
significant for low amplitude waves. 

In general, B and g are expressed in terms of other quantities perturbed 

X z 

by wave activity. These relat 'hips nay be deduced by utilizing the 

polarization relations of the typw > . wave activity being considered. He 
desire to express the layer density response in terns of the wave induced 
atnosi^erie density perturbations. Letting i4exp(i(i)t - iX*r) represent the 
atmospheric density perturbaticns , the polarization t'elatlons are used to find 
the appropriate B,. In this case, B relates the magnitude and phase of the 
atmospheric density perturbations to the horizontal wind, V , and & relates 

X 2 

the magnitude and phase of the density perturbations to the vertical wind, V * 
Several other assumptions are made in obtaining the solutions presented in 
this section. It is assumed that the diffusion time of the minor constituent 
forming the atmospheric layer is much greater than the period of the 
atmospheric waves inducing the observed density perturbations. This implies 
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that tha valooity field of the ainor oonatltuent equals the ataoapberio 
ireloolty field. Only wave induced dynaaios are oonsidered) no oheaioal 
effects associated with atoaio soditai or sources of new atoaio sodiua are 
inoluded in these solutions. This has the effect of reducing the source and 
loss terast P and Q in equation (4. 17)» to sero. Pinallyt it is assuaed that 
the sodiua layer baa a density profile, n (r), in the absenoe of wave activity 
that is a function of positim veotor only. Tiae variations in this 
baokground or steady state profile are not oonaidered. Under these 
oonditions, two expressions for the solution of the continuity equation are 
found. 

4.4.2 THE EXACT SOLUTION 

It is assuaed that the solution of equation (4.17) has the fora 

n(r,t) -s^n-Cr-f e) (A. 19) 

The factor, represents an **aaplifioation" tera which will induce 

perturbations in the aagnitude of n_ • In addition, the arguaent of n is 
perturbed by the factor, 6. e has a horizontal ooaponent, e , as well as a 
vertical coaponent, 6 , since n is a function of both horizontal and vertical 
variables. This tera accounts for vertical and horizontal displaoeaents in 
the layer as it is swept along by the ataospherio waves. Hhen equation (4.19) 
is substituted into (4.17), three differential equations desoribing * 

X 

and 6^ are obtained. 

+ ^*7^) Aexp(i(u)t - X*r)) (4.20) 

dd 

• -(8^ + Aexp(i((iit - X»r)) (4.21) 

36 

■ -(6jj + Acxp(i(a>t - K’r)) 


(4.22) 
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In Moh oas«i m / appro«oh«s stroi « » e i and e mist alae approach aaro. 

X z 

Tha aolutlona to those aquations ara not laMdiataly obvious. Houavary 
powar sarlas taohnlquas aay ba used to solva them. Nota that all thraa 


equations hava tha ganaral fora 


• “(a + ^*Vy) M 


(4.23) 


a and ^ ara ooaplax constants. AS la tha souroa function. In this oasa« It 
Is saan that AS • A^xp^iut - Tha solution of aquation (4.23) Is 

axprassed In powar sarlas fora; 

«D 00 

r- I Y, - I CAAS)^ (4.24) 

l-l £-1 * 

By substituting (4.24) Into (4.23) and oollaotlng teras whloh contain slallar 


powers of AS, the coefficients, C^, aay be datamlned. It Is saan that 

3y, , 

^-£i<^,(M)^ 

and , 

X*Vy^ - -li(\»K) C^(ASr 

Coablnlng these results, It Is concluded that 

-o( 

C, 

. ‘ £(iu) 

*“** « £-1 

r o(tX* K ) ...,x£ 


£(iw)* 


(4.25) 

(4.26) 


(4.27) 


(4.28) 


This result aay be slapllflad soaewhat by utilizing the following 


relationship} 


£nO - a) • - I 


As a result, It is seen that 


~ £n 1 AN 

X’K (It 


(4.29) 


(4.30) 


Is a solution of equation (4.23). 
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This result is now used to find the solutions for end 6 , To 

X is 

obtain the solution of equation (4.20) » tre let !.*§. and To solve 

equations (4.21) and (4.22) i we let 1 ^ and a > 6 and a « B , respeotively. 

X A 


The results of these manipulations are 

4> * -in 1 - 

-iV 


K>V 


in 


K-V 

(j) 

K'V 


1 - 


u 


( 4 . 31 ) 

( 4 . 32 ) 


Thusy the exact solution for the layer density response is given by; 


n(r,t) = 


i V 

K-V 


r - -r;-— £n 

1 . 


— K-V_ 

0) 



u 


K-V 


1 - 


( 4 . 33 ) 


At this pointy the gravity wave polau*ization relations have not been utilized 
so that the solution is valid for a variety of atmospheric waves. The 
numerator of the solution is simply a nonlinear mapping of the steady state 
layer proflley n . The shape of the layer is distorted further by the 
amplification factor; 


K-V 


1 - 


ui 


If wave induced perturbations are small y equation (4.33) may be simplified. 
Specifically, if \K‘^I*h\ « 1 y the layer density response may be approximated 
by 


n(r,t) - 



K-V 


' iV 


1 + 

n 

r + — 



s 

— 0) 


( 4 . 34 ) 


The argument of n is perturbed by a sinusoid with the ssm temporal and 
spatial dependence as the atmospheric wave responsible for the density 
perturbations. If the steady state density profile Is a function of altitude 
(xilyy the argument is perturbed only by the vertical component of the wave 
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Induced wind. Under this condition) a wave with a period of 180 ninutea and a 
■axlauD vei'tical velocity of 1 neter/aecond (auoh aa sone low frequency 
gravity waves) would perturb the argument of n by aa auoh aa 1.72 km. Since 
the aodluB layer is typically 15 to 20 kn wide and contains large gradients) 
such a shift in the layer's features would be quite obvious. Lower frequency 
tiaves and larger vertical velocities would create larger displacements in the 
layer's features. It should also be noted that the maximum upward 
displacement in features occurs one quarter of a wave period after the maximus 
vertical wind. This is a reasonable result. Features are swept upward by the 
vertical wind) attaining their maximum height when the upward wind becomes 
zero. As soon as the vertical wind changes direction (1/4 of a period after 
the maxifflun upweu'd wind) the layer features will begin to be swept downward. 
The undulatory behavior of the sodium layer's topside, bottomside and peaks 
are consistent with the layer behavior indicated by these arguments. Such 
periodic oscillation in the structure of the sodium layer has been observed in 
many Ildar data sets collected at the University of Illinois [Richter et al., 
1981 ). 


In addition to the perturbation in the argument of n , equation (4*34) 

o 

indicates that the layer density is Influenced by the amplification factor, 
(1 + K'V/di ) . For low frequency gravity waves, the phase of this amplification 
factor with respect to the vertical perturbations in layer features gives rise 
to an interesting effect. It is still assuoed that the steady state density 
profile is a function of altitude only. Using the polarization relations for 
low frequency gravity waves, the amplification factor may be related to 
atmospheric density perturbations by noting that 


K»V 


- - j Aexp[£(u)t - ^*r)] 


u 


( 4 . 35 ) 
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where y is the ratio of speoifio heats* Further » the ateospherio density 
perturbations will lag the vertioal wind by 1/4 period for low frequency 
gravity waves* Thusi the amplification factor will lead the vertioal wind by 
1/4 of a period, i.e*, the amplification factor is at a maximum 1/4 of a 
period after the vertical wind is at a minimm and at a minimum 1/4 of a 
period after the vertical wind is at a maximum* One quarter of a period after 
the maximum vertioal wind, the largest upward displaconent of layer features 
occurs as does the smallest amplification factor* Conversely, when the layer 
features are perturbed to their lowest hei^t, the amplification factor 
reaches a maximum* Obviously, the magnitude of the resultant layer density 
perturbations at a fixed altitude are very dependent on the local density 
gradients of the layer* However, the density perturbations below the layer 
peak will generally be larger than those above the peak as a resul^ of the 
relative phases of the amplification factor and the perturbation in the 
argument of . 

In addition to the Imbalance in the size of density perturbations above 
and below the layer peak, the layer density perturbations below the peak are 
180 degrees out of phase with atmospheric density perturbations* Above the 
steady state layer peak, these perturbations are in phase* Of course, these 
are general observations and the specific layer response can only be 
determined when an exact description of the layer is known. The phase 
reversal in the layer density response mentioned above is easily explained if 
a steady state density profile which is symmetric about a single peak is 
assumed. The atmospheric density perturbation is at a maximum 1/4 of a period 
after the maximum vertical wind. At the same time, the layer peak reaches its 
maximum upward dlsplacMient . Above the steady state location of the layer 
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peaki this haa the effect of enhancing the density of the minor constituent* 
Belou the layer peak} the layer density is depleted as the low density 
portions of the layer are swept upward by the wind. The net effect is that 
when the atmospheric density perturbation is at a maximum, the minor 
constituent density reaches a maxlmua above and a minimum below the layer 
peak. k similar argument may be used to show that when the atmospherio 
density perturbation is at a minimum, the layer density perturbation will be 
at a minimum above the layer peak and at a maximum below the layer peak* As a 
result, the density perturbations of the minor constituent 2 u*e generally 180 
degrees out of phase with atmospheric density perturbations below the layer 
peak and in phase above the layer peak. 

Although equation (4*33) is a concise expression for the exact solution of 
the continuity equation and yields insight into the layer response, it is not 
necessaurlly the simplest formulation of the solution to work with* Another 
form of this solution may be obtained which separates the linear and nonlinear 
portions of the layer response. This is accomplished by expanding the 
solution given in equation (4*33) in a Maclaui*in series. It comes as no 
surprise that the results obtained in this manner are identical to the 
complete perturbation series solution presented in the next section. 

4.4.3 THE PERTURBATION SERIES SOLUTION 

In order to separate the linear and nonlinear portions of the layer 
response, a series solution of equation (4.17) is sought which has the 
following form; 

00 00 

» I “ I P»(AiV)*' 
t-0 i-0 


(4.36) 
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As In section 4.4*2) M * Aexp(t(i>t • iK»p)» This solution iS) in fact) the 
perturbation series solution of equation (4.17). It is apparent that will 
be the linear response and *>■2 will represent the nonlinear portion of the 
response associated with a double frequency oscillation in density. In 
general) the term will describe the portion of the layer response 
associated with density oscillations whose frequency is the harmonic of 
the atmospheric wave. 

To determine P^) equation (4.36) is substituted into (4.17) and terms of 
equal power in LN are collected. This is the same technique used in obtaining 
the power series solution for the exact layer response (section 4.4.2). 
Unlike this example) however) the coefficients) P , are functions of the 

X 

position vector) r. By applying the boundary condition 


ilim - n (r) 


(4.37) 


it is determined that P = n (r) • The ooef f icients of higher order terms must 

u o 


satisfy the recursion relation; 
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and the series solution is given by; 
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(4.39) 


In order bo simplify the result) it is asstmed that n (r) - n ( 2 ); the 

S 8 

density profile in the absence of wave activity is a function of altitude 
only. This is not an unreasonable restriction. Without this restriction on 
the horizontal variability of n ) the expression for the series solution 

o 

becomes quite cmbersome. The simplified result is; 
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vrtiere 
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Therefore, the first three terms In the series are; 
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(4.41) 


(4.42) 


expt2i(ut - X*r)l 


(4.43) 


Note that derivatives of n (a) up to the order are contained In the 

8 

term. 

The polarization relations for gravity naves may now be used to determine 
the nature of the layer density response to gravity waves. Since low 
frequency gravity waves are apparently responsible for mmiy of the dynamic 
features observed In the sodium layer, the polarization relations for low 

frequency gravity waves [Hines, I960] are used to relate b and g to wave 

X z 

Induced perturbations In atmosiNierlc density. Using these relations, which 
are expressed In equations (4.10) and (4.11), It oan be shown that; 
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(4.44) 
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The linear layer responae obtained by utilizing theae relation la 
Identloal to the linear layer reaponae reported In prevloua reaulta [Chiu and 
Chlng) 1978; Welnatoek, 1978]* 


-A 


n. 


Y - 1 


n dz 
8 


exp[i(wfc -^*r)] 


(4.46) 


Experimental evidence Indicating auoh a reaponae In ataoapherlo layera haa 
been reported for the meaoapherlo aodlum layer [Shelton et al*» 1980]) In 

) alrglow meaaurementa [Welnatoekf 1978]* and In nl^tglow emlaalona of 
the hydroxyl radical [Prederlok, 1979]* 

Several of the oharacterlatlc featurea of the linear layer reaponae are 
aummarlzed below* The denalty perturbatlona In the atmoapherlc layer 
aaaoclated with the linear reaponae can be aubatantlally larger than the 
atmoa;^erlc denalty perturbatlona aaaoclated with the gravity wave. Thla la a 
reault of the advectlon of denalty gradlenta* In general, large denalty 
gradlenta will give rlae to large "ampllfloatlcn" faotora, l*e*, large 
perturbatlona In the denalty of the minor conatltuent. Aa a reault of the aum 
In the rlghthand aide of equation (4*46), denalty perturbatlona below the 
layer peak (where dn /da la generally poaltlve) will typically be larger than 

3 

thoae above the peak (where dn / da la negative)* In addition, a phaae 

® yH dn 

reveraal In the layer reaponae oocura at the point irtiere — = -1 * Below 

8 

thla point, atmoapherlo and aodlun denalty perturbatlona are 180 degreea out 
of phaae tdille above thla point theae denalty perturbatlona are In phaae* 
Recall that thla type of behavior waa deduced by examining the exact aolutlon* 


Since the nonlinear part of the layer reaponae doea not generally approach 
zero at the point of phaae reveraal, the layer reaponae will uaually be 
dominated by the nonlinear terma In thla region* Evidence of auoh 
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nonllnearltlas in the sodium layer response has been reported [Shelton et al.» 

1980 ]. 


4.4.4 SPATIAL POWER SPECTRUM OP THE LINEAR UYBR RESPONSE 

The vertical structure associated with the linear layer response 
complicates the Interpretation of Ildar data. One of the consequences of the 
phase reversal In this response is that peaks in the sodl\SB layer will not 
necessarily be separated by distances equal to the gravity wave's vertical 
wavelength. The spatial power spectna will reflect this fact. If it Is to 
be used as an Interpretive ald| the effects of the vertical structure of the 
layer response on the power spectrum must be understood. If only the linear 
portion of the layer response to gravity waves la Included, the sodium layer 
density may be expressed as; 

r ii 1 r 

n(r,t) - n (a) ^ n (a) + yH — t •exp[t(o)t - K x - K z)\ (4.47) 

8 A j ^ u2 J X Z 


The spatial Fourier transform of this equation Is readily found. The power 
spectrum due to the linear layer response Is simply the square of the 
magnitude of this transform. The transform will be e.xprossed as a function of 


the vertical spatial frequency, f , where K 

Z Z 


2vf^. The following equations 


are written as functions of spatial frequency rather than wave muber because 
estimates of the power spectrum calculated from experimental data are 
traditionally plotted as functions of spatial frequency. The Fourier 

transform relating the spatial variable, a, to /_ is evaluated yielding the 

z 


result; 




exp[t(u)t - x)l 
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(4.48) 
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or 

- K^*n 

•II + - /■,) r*l *,(r, - /,) (4.W 

It and N are the Fourier tranaforms of n and n reapeotively and /' = i/x 
uhere la the vertioal wavelength of the gravity wave. The tranafora of the 

A 

linear portion of the layer reaponae ia in the ri^thand aide of each of theae 

equationa* The linear reaponae, aa expreaaed in equation (4.46), la the 

product of a complex exponential, exp(iwt - iK»r)f irtiich deacribea the time 

and apace dependence of the ainuaoidal variationa in denalty and an envelope, 
-1 **”8 

(y _ , Which determlnea the vertioal atructure of the magnitude 

of the denalty oaolllationa. The tranaform of the linear reaponae la the 
convolution of the tranaform of thia envelope and the tranaform of the complex 
eiqwnentlal . 

The power apectrum ia eaally found. 

2 

+ [74^) 111 

- (t^t) s.u>n>Ci(»t - vmi * “•W'. - !'*> ".W’P W - ( 4 . 50 ) 


In practice, the average apatlal periodogram of experimental data la 
evaluated. In order to make direct compariaona convenient, the average 
apat.lal power apectr>jm (ASPS) la found by averaging equation (4.50) over one 
period of the gravity mve. The reaulting equation ia; 


ASPS - k (/ )r + 


8 


AJi‘ 


Y - Ij 


1(1 + i2v(/^ - Q yH) 




(4.51) 


Thia expreaaion la very dependent upon the tranaform of the ateady atate layer 
profile, iV.(/„) . The nature of thia dependence may be examined by evaluating 

C ** 
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a apeoial oaae. 

The steady atate sodium layer density profile may often be approximated by 
a gauasian shape. Thus, it is assumed that 


« ( 2 ) « Pexp 

8 


“U - Sq) 




2o" 


(4.52) 


tfhere P is the peak density, Zq is the altitude of the layer peak and a is 
related to the layer width. The ASPS for such a layer becomes; 

ASPS - 2noV[exp(-4itV/^) + (l+4Tr^(/; - yV) 


•exp[-4w^a^(/' - f )^J 

The second term in the above equation peaks at 




(4.53) 


This behavior may be understood by examining equation (4.48). For a gaussian 
layer, the peaks of (1 + i2vf yH)N(f ) occur at /' - ±Af . These peaks are 

2 2 Z a 

shifted by the convolution in equation (4.48). A plot of the calculated ASPSs 
for several values of A and several vertical wavelengths are shown in Figure 
4.4. Aside from the lack of a constant noise level, these ASPSs are very 
similar to the periodograms used to estimate the power spectra of experimental 
data. One of the unusual oonsequenoes of the vertical structure of the layer 
response is that a local minimus oocurs in the ASPS at the spatial frequency 
corresponding to the vertical wavelength of the gravity wave. At this 
frequency, the contribution to the ASPS from the linear layer respc-nse la 
approximately (A/(y - i))2. 
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Plots of the simulated average spatial power spectra resulting 
from the linear layer density response to gravity wave 
perturbations. The wave amplitude Is given by A, The vertical 
wavelength of the gravity wave is (a) 10 km* (b) 5 km, and 
(c) 3 km. 
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4.4.5 COLOHN ABUNDANCE VARIATIONS 

The Fourier transform of the layer response as given in equation (4.49) 

has another interesting use. By evaluating this expression at f - Of the 

z 

column abundance of the layer may be tooui. However t high order terms (i.e. 
nonlinear effects) also contribute to column abundance variations. If the 
perturbation series solution (equation 4.40} is integrated over the spatial 

variable! Sf an expression for the total oolwn abundance la found. 

00 00 

C(«,t) ■ Re|j dsn(r,t)| « Re|| dz I (4.55) 

Because n(r.t} is zero for s < 0, the lower limit of integration may be 


1 r r r 


f • ReM 1 
^ 'm»0 £*0 

u 


“6 K 

B B 


-th 




(4.56) 

The m'’“ order term is responsible for variations in column abundance with a 
temporal frequency cf mu . The linear term (ms 1 ) dominates the variations in 
column abundance. This becomes obvious when the magnitudes of the first 


several terms in this series are examined. Recall that K 


3 


2,f^. 


<^0 ■ W 


0 ) 


(4.57) 


W ■ (4.58) 

(1 - 2.V;\V . jr^Cf, - 2/;) 


<^2 ■ 

The gravity wave polarization relations have been used to obtain the above 
expressions. The factors responsible for the rapid decrease in the magnitude 

of C with increasing m are the dependence of C on (^/ (r - D)'” and the value 

m w 


of N at f - mf 

8 B Z 


The magnitude of the gravity wavef A t is quite small so 
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that it'" rapidly daoraaaaa with Incraaslngm* Typioally, W falls off rapidly 
with inoraaaing For axasplOt if ia gaussian (as given in aquation 

(4.52)), then 

^^(/g) • ^ oP expl-2irV/^] exp(-t2v/^SQl (4.60) 

As a result, the magnitude of deoreases with increasing m at the 

approximate rate of 

exp(-2vV(m/^)2l (4,61) 

due to these factors. This decrease is opposed, to some extent, by an 

increase in the size of the coefficient ot N (f “ nf) in the expression for 

s z z 

(equation 4.56)). 

ffl 

Because the magnitude of the term, C^, deoreases rapidly as m 

increases, the column abundance is usually described adequately by the first 
few terms in equation (4.56). As a result, the approximate column abundance 
is obtained by utilizing only the first three terms from this equation. 


C(x,t) s + + C^} 


(4.62) 


It is evident that C(x,t) varies sinusoidally about Ke{N (0)}. For a 

s 

gaussian layer given by equation (4.52) the column abundance becomes} 

C(x.t) - ^ op[l - exp(-2vV/;^j 

•(1 + 4vVlSV)^^^ co8(«t - Kje + tan"^(2ii/'Yi?) - 2irra.) 

A X Z Z V 

+ (tTt) - 2vV;W)^ + 

(4.63) 

This expression oontains the sun of three terms; a constant (the zero order 


cos 


-1 1 
2o.t - 2K x + tan ^ 

ll - 2«V\V ^ ® 
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tera), a tara with a taaporal fraquanoy of » (tha first ordar tara)» and a 
tara with a taaporal fraquanoy of 2u (tha aaoond ordar tara)* Tariationa in 
ooliam abundanoa dua to tha linaar layar raaponaa ara asaooiatad with tha 
first ordar tara. Tha aagnituda of tha linaar variation in ooluan abundanoa 
noraallsad by tha avaraga ooluan abunduioa ia givan by} 

AC- - «xp(-2wV/'^J (1 + (A. 64) 

AC^ is tha aagnituda of tha noraallzad variation in C(x,t) dua to tha linaar 

layar raaponaa about Ita aaan value. Pigura b.5 shows how tha ratio of AC^ to 

A varias with and o. For a gaussian layar, tha largast oolunn abundanoa 

variations ooour whan /' ■ — ^ ^ — ]*! At this point, tha ratio of 

* 2ir [o2 

AC to A is approxiaataly 0.6065 'yH/iaiy - 1)). If o*3» this oorrasponds 
to a gravity wava with a vartloal wavalangth of 19*6 kn and ooluom abundanoa 
variations whose aagnituda is roughly 5.2 tiaas tha aacni.tuda of tha 
ataospherlo dansity variations. It is clear froa Figure 4.5 that large 
variations in ooluan abundance will be evident whan tha layar is narrow (i.a. 
a is saall). 

Tha magnitude of tha ooluan abundanoa variations assooiatud with the 
saoond ordar t«>ra in tha layar response (aquation 4*59) for a gaussian layer 
ara 

AC • exp(-8v Vr^Hd - 2 vV;^yV)^ + 9vV'^yV)^^2 (4.65) 

Figure 4.6 shows how tha ratio of AC, to A^ varias with o and X . Although 

£ z 

AC2 is usually auoh saallar than AC^ , for large wava aaplltudas and vartloal 

wavelengths, nonlinaaritias beoona aora iajwrtant* Figure 4.7 shows tha 

■axinua wava aaplitude, as a function of X and o for which AC, *iAC.. It is 

Z I z 
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Figure 4.6 Plots of the ratio of the normalized perturbation in column 
abundance due to the second order layer response (&C 2 ) to the 
Have amplitude squared A gaussian layer with width a has 

been assumed. 
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clear that can aafely be nagleotad for ka and ^ 3 ^ 15 ka* 

Beoauae of the preaenoe of high order teraa In equation (4*56) » oolian 
abundanoe variatiraa tilth teaporal perloda leaa than the period of the gratrity 
cave nay add aignifloantly to Aa la apparent froa evtuatlon (4.63) t 
the relative aagnitudes and phaaea of theae nonlinear teraa are very dependent 
on the layer ahape and the vertloal wavelength and aaplitude of the . gravity 
wave. To illuatratOi the peroent ooluan abundanoe varlationa due to and ^2 
for a gauaalan layer are plotted in Figure 4.8. Aa tho layer beooaea narrow 
and aa the wave aaplitude Inoreaaeai the relative oontributlon of the firat 
nonlinear tern, Inoreaaea. 

4.4.6 NONLINEAR NATURE OP THE LATER RESPONSE 

The degree to whioh the aodlun layer reaponae ia nonlinear depends upon 
aeveral factora. Since the layer denalty reaponae ia highly dependent on the 
density gradients ooourlig in the layer, the steady state layer density 
profile beooaes very laportant. Large density gradients encourage 
nonlinearities in the layer response. Two other laportant factors are the 
aaplitude and the vertloal wavelength of the gravity wave inducing the density 
perturbations. The degree to whioh these two factors influence the nonlinear 
nature of the layer response has been investigated* Reoall that the 
polarization relations for low frequenoy gravity waves are being used. Since 
the sodlua layer shape can often be approxiaated by a gaussian function, this 
shape was used as the steady state sodlua density profile for slaulations in 
whioh the gravity wave aaplitude and vertloal wavelength were varied. In each 
case, a gaussian density profile centered at 90 ka with c • 3 ka was assused. 
The dependence of the total RHS error on vertloal wavelength and wave 


A«0.05 

A-0.03 

A-0.01 





Figure 4.8 The percent varletlons In sodium column abundance due to linear 
and second order density perturbations. The variations are 
plotted for 2 periods of the gravity wave, k gausslan layer 
with width o has been assumed. The dependence of A ^2 on 

A and a Is. Illustrated. 




89 


•■plltude Is given In Figure 4.9* Wave eaplltudes ere given by the peak 
percent ataospherlo density variations associated with the nave propagation. 
The total RMS error was calculated by taking the square root of the ratio of 
the Integrated) squared nonlinear layer density response and the Integrated) 
sipiared total layer density response as Indicated In the Figure. Equation 
(4.33) was utilized to compute the total layer response and the nonlinear 
layer response was obtained bf subtracting the linear response from this. The 
time Integrals were performed over one period of the gravity wave) and the 
limits of the spatial Integral were 70 and 110 km. Outside of this altitude 
range the density perturbations are negligibly small. If the layer response 
were completely linear) the total RMS error would be 0%» Similarly) If the 
layer response were completely nonlinear) the total RMS error would be 1009. 
Figure 4.4 Indicates that the degree of nonlinearity Is only slightly 
sensitive to vertical wavelength) but very sensitive to wave msplltude for 
gausslan steady state layers. Even for gravity waves which create only 
peak atmospheric density variations) at least 27% total RMS error will be 
Incurred by assualng a strictly linear response. 

As shown before) the nonlinear portion of the layer response will tend to 
dominate In restricted reglmis of the layer. Fl^e 4.10 shows an IMS error 
analogous to that of Figure 4.9. In this case the altitude Integral Is not 
performed and the resulting error expression Is a function of altitude) wave 
amplitude) and vertical wavelength. Since the errors are relatively 
Insensitive to changes In vertical wavelength) a vertical wavelength of 10 km 
Is assuaed and the simulation results are displayed as r. funotlon of altitude 
and wave amplitude only. 
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Pigur* 4.9 








Plgur* 4.10 shows thst ths RHS srror pssks to 100S st 90*92 tos. This 
oorrssponds to ths point of phsss rsvsrssl in ths linssr Isysr rssponssy i.s. 
ths point whsrs ths linssr Isysr rssponss is ssro. In gsnsrslt ths RMS srror 
is snallsr bslow thsn sbovs ths Isysr pssk (whiOh ooouni st 90 loi) bsosuss ths 
linssr rssponss is Isrgsr bslow ths Isysr pssk. Ths RMS srror spprosohss 1001 
ss ths distsnos fron ths Isysr pssk inorsssss* nis nonlinssr tsm in ths 
psrturbstion ssriss solution doninsts ths linssr tsm in thsss rsgions. Hots« 
howsvsr» thst ths sodiun drasity psrturbstions boooas vsry sasll st thsss 
sxtrsBSS whsrs ths Isysr is quits tsnuouo. Pigurs 4.11 illustrstss ths sizs 
of dsnslty psrturbstions sssooistsd with snd . Ths snsslops of 
psrturbstions for ssoh tsm hss boon plottsd. Ths gravity wsvs saplituds wss 
ohossn to bs oonsistsnt with 51 pssk to pssk vsristions in staosidisrio 
dsnsity* It is quits svidont thst ths ssoond ordsr tsm» doainstss tho 
Isysr rssponss in ths 90.3 to 91*2 ks rsgion. Hots» howovor* thst ths third 
ordsr torst is Isrgsr thm ths ssoond ordsr tsm st tlaos. 
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5. SODIUM LXDAB OBSBRVATIONS AT UBBANA, ILLINOIS 

5.1 BACKOROOND UFOIIHATION 

Ltridmoe o/ mv« AOtltrity n oft«n found In lldnr observations of the 
sodluB layer aade at Orbana, Illinois. These observations are aade with the 
Ildar systea pointing at zenith. Because low frequency ataos^erlo waves my 
be responsible for aany of the observed layer dynaalosi data sets spanning 
long observation periods are essential. In soae data, density osolllatlons 
with periods as long as 8 and 12 hours are present. Extended observation 
periods are necessary to detect these low frequency oscillations. The length 
of the observation periods varied with each data set. The period of 
observation as well as other pertinent paraaeters are auamrlzed In Table 5.1. 
The **0 Indicate that data was collected with the original laser. If no * Is 
present, the Candela laser was utilized. Not all data were oalilwated. In 
some oases, overloading distorted all Rayleigh scattering mklng calibration 
Impossible. In each case where calibration could be achieved, the time 
variations In column abundance are presented. After March 11-12, 1981 the 
overloading problem was permanently solved by gating the PMT. This was 
accomplished by electronically switching the first fUT dynode voltage as 
discussed In Chapter 2. As a result, all data collected after Naroh 11-12, 
1981 have been calibrated. 

5.2 DATA COLLECTED WITH THE ORIOINAL LASER 
5.3.1 OCTOBER 13-lA, 1979 

On tiM evening of Ootober 13-14, 1979 the Ildar system operated for almost 
3 hours oolleotlng 49 spatial profiles at 3 minute intervals. During this 
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period 5 profllM mt* lost du« to oqulpMnt fdiluro. Thoao profilM iioro 
roploood with MtlMtoo obthlMd by oppi/lng o linoor intorpolotion proooduro 
to tho riMlnlng proflloo. Tho Avtrofo Toaporal Porlodograa (ATP) la ahown in 
Pifura 5.1. It waa diaoovarad that tho ATP was only alifhtly affootad by 
laaving tha 5 intarpolatad profilaa out antiraly. Thia indioatad that no 
aisnifioant arror waa introduoad by uaing intarpolatad data to raplaoa tha 
■iaaing profilaa. 

Tha ahot noiaa floor In tha ATP ia 10 dB btlow tha paak at 0.0063 
■in'^. Thia paak oorraaponda to a taaporal pariod of roughly 120 ainutaa. 
Raoall that tha low fraquanoy portion of tha ATP ia diatortad aoaawhat by tha 
affaota of windowing. It la olaary howavary that tha doninant taaporal 
varlationa in aodlua danalty hava long parioda. Baoauaa thia data aat ia 
fairly abort y tha apaotral raaolution of tha ATP ia not vary fO*aat. Tha firat 
null in tha DPT of tho window oooura at 0.012 ain**^. Aa a r^aulty thara la 
aignifioant diatorti«n at fraquanoiaa lowar than thia. Tha natura of tha 
diatortion ia daaorlbad by aquation (3.6). In addltiony tha oonvolutlonal 
aaaaring of tha ATPy aa daaorlbad by aquation (3.6) y la fairly aavara. 
Although fina datail ia obaourad by thia anaaringy a Mttoff for tha taaporal 
filtar nay atill ba ohoaan froa tha ATP. In thia oaaa tha ahot noiaa laval 
doainataa tha ATP at fraquanoiaa graatar than 0.02A ain*^. Baoauaa tha 
paaaband of tho ta^^al filtar la not naty a outoff at O.OAAt ain*^ waa 
ohoaan. Thia inaia^ that tha gradual rolloff of tha filtar paaaband did not 
aignifioantly altar tha algnal apaotriai bolow 0.02A aln~^. 

Tha Avwaga Spatial Pariodograa (kSf) of tha taaporally filtarad data ia 
ahown in Pigura 5.2. Thara la a aignifioant nount of anargy at tha apatial 
fraquanoiaa 0.166 and 0.255 Iob*^. At thaaa pointay tha va?.ua of tha ASP ia 
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Figure 5* 


Figure 5. 
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TEMPORAL FREQUENCY (min') 


1 Plot of the normalized average temporal periodograa of data 
collected on October 13-14, 1979* The temper^ frequency 
associated with the gravity wave (0.0083 aln**^) is indicated 
at a. 



Plot of the normalized average apatlal periodogram of temporally 
filtered data collected on October 13-14, 1979* The spatial 
frequency associated with the gravity wave (0.21 km~') la 
indicated at b. 
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about 22 dB less than the peak value. It Is olesTf therefore, that features 
with vertical wavelengths of 3*9 and 6*0 km are present. This ASP is <iuite 
similar to the simulated power speotm plotted in Figure 4.4(b). Taking into 
account the nature of the linear layer response as developed in Chapter 4, it 
appears that a gr*avity wave with a vertical wavelength of X is responsible 

B 

for these spectral peaks, irtiere 1/x s.5(0.25St>0.l66) ka*^. The amplitude of 

s 

this wave may be estimated by noting that the magnitude of the periodogram at 
f « 1/X^ is roughly 3*1 X 10**^ and, therefore, the amplitude of the gravity 
wave is approximately (y-D/ 3*1 X 10”3, Thus, a gravity wave with a vertical 
wavelength of roughly 4.7 km associated with 4f peak to peak variations in 
atmosphci'lc density may be present. Since the shot noise level dominates the 
ASP for frequencies greater than 0.5 km~^ , a spatial filter cutoff at this 
frequency was chosen. This filter was applied to the temporally filtered data 
to yield the final two«dlmenslonally filtered data. 

A plot of the spatial variations in estimates of the sodium density is 
contained in Figure 5.3* The layer is fairly narrow; it does not extend below 
83 km nor above 102 km* Further, the main layer peak gradually moves from 90 
km at 2330 to 88 km at 0112 CST. The peak InltUIly located at 87 km moves 
steadily downward with a velocity of 0.64 m/s. Althou^ the relative 
amplitudes of the layer features above 90 km change, it is more dlffioult to 
discern the same consistent motions above the layer peak. The diagonal lines 
drawn in this figure indicate the apparent phase progression of the wavelike 
features. These lines are separated vertically by the previously deduced 
value of 4.7 km. Below 90 km peaks in sodium density lie on the lines. Above 
90 las, these peaks become valleys. The valleys appear to share the peaks' 
downward velocity of 0.64 m/s* Note that the peaks below 90 km are more 
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Figure 5*3 Tlae history of th*. estlasted altitude profiles of sodiiaa 
density observed on October 13>14, 1979* The spatial and 
tesporal filter outoffs we.*e 0.5 )(■'' and 0.044 sln~'t 
respeotively. The diagonal lines Indloate the apparent phase 
progression. 
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pronouno«d tban the valleys above 90 km. It appears that density 

perturbations beloir the layer peak are greater than those above the peak. 
This ohange In the magnitude of wavellke features and the observed phase 
reversal near the layer peak are oharaoterlstlo of a layer response to gravity 
waves. 

Another way of displaying sodium density variations Is shown In Figure 
5.4. The t«nporal variations In ooHm density about their means at speolflo 
altitudes have been plotted. Baoh of these temporal i»*oflles Is then 

normalized so that the peak to peak variations are equal at eaoh altitude. A 
downward phase progression In the density variations below 90 km Is evident. 
Associated with this wave motion Is a 120 minute period and a 4.8 km vertloal 
wavelength. The wavelength Is consistent with the ASP and the downward jrtiase 
velocity matohes the velocity of the layer features observed In the spatial 
profiles of sodium density. A gravity wave with this temporal period and 
vertloal wavelength would have a horizontal wavelength of approximately 120 
km. At an altitude of 89.1 km# a double frequraoy oscillation tflth a period 
of 60 minutes Is notloeable. nils Is consistent with the prediction that 
nonlinearities should be present In the layer response » particularly near the 
layer peak. Above 90 icSf these periodic fluctuations In sodlua density are no 
longer as apparent. Instead, a gradual decrease In density Is observed. This 
may be caused by t«sporal changes In the steady state sodlia profile or it may 
be a layer response to atnosjdierlo waves of mudi lower frequencies than can be 
observed In this short data set. 
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Figure 5*4 



2330 ' 2357 ' 0027 ’ 0057 
2342 0012 0042 0112 


TIME (COST) 


Plot of the temporel verlatlone In eodiiai density observed on 
October 13-14, 1979> The SMtlal and taaporal filter outoffa 
tiere 0.5 ha"' and 0.044 ain”', reapeotlvely. The diagonal lines 
indioate the apparent phase progression. 
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5.2.2 OCTOBER 28-29» 1979 

The llder system operated for approximately 6 hours on the evening of 
October 28-29 t 1979. 71 spatial profiles were oolleoted at intervals of 5 
minutes. The ATP for these data is shown in Figure 5.5. The effects of the 
temporal window on this ATP have been disoussed in Chapter 3. Two distinct 
peaks near 0.0059 min*^ and 0.0118 min~^ are obvious. These fre<]uenoies 
correspond to temporal periods of 170 and 85 minutes » respeotively. The faot 
that the seoond frequency is twice the first suggests that nonlinearities in 
the layer response to wave activity nay be observable. A temporal filter with 
cutoff at 0.033 min~^ was applied to the photooount data. The ASP of the 
temporally filtered data is shown in Figure 5.6. A single large sldelobe is 
iresent. If this sldelobe is attributed to the layer response to a gravity 
wave* it would appear that the vertical wavelength of this wave is greater 
than or equal to the wavelength associated with the local minlnvn in the ASP 
at 0.102 km**^. Thus, it is probable that the vertical wavelength of the wave 
Is greater than 9.8 km. The portion of the ASP due to the steady state layer 
profile dominates the portion of the ASP resulting from the layer response to 
waves since the latter is confined to small spatial frequencies. As a result » 
the exact vertical wavelength cannot be determined from an examination of the 
ASP. This is apparent in the simulated power spectrum shown in Flguie 4.A(a). 
However » the ASP is used to select a spatial filter outoff at 0.44 km'^. 

TWo-dimenslonally filtered spatial profiles are presented in Figure 5.7. 
In order to conserve space) estimates of the spatial variations in sodium 
density are plotted at 15 minute Intervals. In addition) a three dimensimial 
plot of these data is presented in Figure 5.8. Below the layer peak) which 
occurs near 90 loS) a great deal of aotivity la present. The peak) which is 
initially located near 88 km) moves down to 84 km over a period of 4.25 li?urs. 
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TEMPORAL FREQUENCY (min') 


Pigur« 5<5 Plot of the noraellzed average teaporal periodograa of data 
oolleoted on Ootober 28-29t 1979* 



Flffire 5*6 Plot of the noraellzed average apatial periodograa of teaporallf 
filtered data oolleoted on Ootober 28*29 f 1979* 


ALTITUDE (km) 


OCTOBER 28-29, 1979 



TIME (CST) 


Figure 5«7 Tine history of the estiaated altitude profiles of sodlua 
density observed on October 28-29» 1979* The spatial and 
tanporal filter cutoffs were 0.44 kB~' and 0.033 ■In'* f 
respectively. Profiles are plotted at 15 minute Intervals. 
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MMr 0100 CST» this stoondary paak haa alaoat dlsappaarad. Aa a raaoltt tha 
layer beoosea narrower and than broadens as the secondary peak reappeara. 
Tseporal variations In the two-dieenslonally filtered data are shown in Pigure 
5«9> A downward phase progression of wavelike features is apparent. 
Assoolated with this phase progression is a vertioal wavelength of 12 Ion and a 
tsnporal period of 170 ninutes. Both the temporal and spatial perlodograns 
are oonslstent with these parameters. Near 91 km density oaolllations with a 
period of 85 ninutes appear. These are probably assoolated with the second 
order (double frequenoy) portion of the layer response and are related to the 
peak in the ATP at 0.0118 nln"^. They appear at this altitude beoause the 
linear response becomes small Just above the layer peak near the point of 
phase reversal. At the point of phase reversal, the linear layer reapcmse 
vanishes. This idiase reversal, as discussed in Chapter 4, is apparent in 
Figure 5.9. 

Figure 5.10(a) contains a plot of the maximum and minimum sodiimi densities 
observed during the evening's data collection (dashed lines). The average 
altitude profile of sodliai density is shown by the solid line and is seen to 
be roughly gaussian. A gaussian layer with os3 km centered at 90 km was 
utilized to predict the maxlmvD and minimun densities associated with the 
linear layer response (Figure 5.10(b)). There is good correlation between 
Figures 5.10(a) and (b) in several respects. The amplitude of the density 
variations below the layer peak are greater than those above the peak. For 
the gaussian layer, the largest dsnslty perturbations occur at 87.6 km and are 
9.1 times largsr than the atmospheric density variations. The magnitude of 
the density variations in experimentrl data are approxiaately tha same. This 
indicates that ths amplitude of the atmospheric density variations associated 
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Plot of tbo toaporol vmriotlons In sodliai donnitjr oboortrod on 
Ootobor 28-29 1 1979* Tbo s^tlnl and twpornl flltor outoffa 
iiara 0.44 loa*^ and 0.033 nln’^ raapaotivaly. Diafonal llnaa 
Indioata tha apparant pbasa prograaaion. 
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Pifur* 5*10 Plot of tbo oovolopo of sodlta donslty vurlotiOM (daobod lines) 
•bout the averafe daaalty (solid lloa) obaareed on Ootober 
28-29* 1979 (a) and tba aiaulatad oneolopa of danaity 
perturbatlooa (daahad llaaa) for a fauaaiaa ataady state layer 
(aclld line) (b)« A gravity uave aaplltude oeoalateat with 3f 
peak ataoapbarie density variations (A a 0*03) ana aaauaad. 
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with th« fPAtrity wavt ar« protebly n«hr 31 «0.03)* TIm oloMat 
approach of tha ounraa doaorlblnf aaxirai and ainlaiai danalty profllaa in 
Plfura 3.10(a) oeoura at 91.5 la (1 )■ ahova tha layar paak at 90.9 la), nila 
point probably oorraaponda to tha point of phaaa ravaraal. Tha two ounraa do 
not aaat at thia point baoauaa hichar ordar taraa (aapaoially tha aaoond ordar 
tara) in tha layar raaponaa ara not aaro at thia altituda. Tha altituda of 
tha point of phaaa ravaraal daduead abova oorraaponda to tha point of phaaa 
ravaraal obaarvad in tha plot of tMporal variationa In danaity (Pifura 5.9). 
For tha gauaaian layar » tha point of phaaa ravaraal oooura 0.92 km abova tha 
layar paak. Thia ia approxljHktaly tha aaaa aa tha aaparation of tha layar 
paak and tha point of phaaa ravaraal daduoad froa axpariaantal data. 

5.2.3 DBCBffiBR 13-1t, 1979 

On Daoeabar 13>1t» 1979 a aariaa of 58 apatial profilaa wara oollaotad at 
tiaa intarvals of 5 winutaa. Tha ATP of thaaa data la ahoun in Pifura 5.11. 
Baoauaa tha data apannad laaa than 5 hourat tha low fraquanoy diatortlon 
introduoad by windowing tha ATP auat ba oonaldarad. Tha firat null in tha OPT 
of tha window oooura at 0.0069 ain*^. Tha paak of tha ATP» whioh oaoiva at 
0.0031 ain'^i ia wall within thia ragion. Aa a raault* tha taaporal fraquanoy 
of danaity oaoillationa cannot ba praoiaaly datarainad by axaainlng tha ATP. 
It ia olaaTf howavar* that tha o^?lllationa *hava a ^ry low fraquanoy. Tha 
pariod of tha oaoillationa ia probauly in axoaaa of 300 ainutaa. Thara ia 
algnifioant anargy at taaporal fraquanoiaa balow 0.022 ain'^. At hlghar 
fraquanoiaCf tha ahot noiaa floor doainataa tha ATP. In ordar to paaa 
apaotral ooaponanta with fraquanoiaa laaa than 0.022 ain*^ without algnifioant 
diatortionf a taaporal filter outoff at 0.033 ain**^ waa aalaotad and applied 
to tha photooount data. 
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Fleur* S.lt Plot of th* nomllMd *v«r*e* t«ipor«l ptrlodotrM of dat* 
ooll*ot«d on D*o«ib*r 1979* 



Pifur* 5*12 Plot of tb* norMllMd *v*r*a* apatial porlodogris of taaporallp 
filtar*d data oollaotad on D*o«ibar 13-14, 1979* 
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The ASP of the teaporelly filtered data is ahoim in Figure 5* 12. This 
plot la interesting beoause it indioates that a series of peaks in the ASP 
exist uhioh groir aaaller trith Inoreasing spatial frequency. It is difficult 
to explain this in terns of a single gravity wave unless the aagnitude of the 
nonlinearity in the layer response is large. The aagnitude of the higher 
order ooaponents in the layer response are quite dependent upon both the nave 
anplltude and the steady state layer profile. As the wave aaplitude and the 
density gradients in the steady state layer beoone large , so do the higher 
order ooaponents of the layer response. Other possibilities are that several 
waves are present or waves are "breaking" into structures with saaller soale 
due to nonlinear! ties in the propagation of the waves. These last 
possibilities have not been treated in this thesis. On the basis of the ASP* 
a spatial filter cutoff at 0.3^5 ka~^ was selected and applied to the 
teaporally filtered data. 

Plots of two-diaensionally filtered spatial profiles are contained in 
Figure 5.13. Estinates of the profiles of sodliai density are plotted every 10 
ainutes. The layer Is observed to be fairly broad. It generally extends frota 
77 to 103 ka. The Initial bifurcated layer grows into a single peaked layer 
over a period of 2.5 hours. Near the end of the observation period, the layer 
is once again beo<»lng broader as the wavelike features below the layer peak 
grow in anplltude. The changes in layer width can be explained in terns of 
the layer response to gravity waves. Beoause of the phase reversal in the 
layer response, waves with vertloal wavelengths of the sane order as the layer 
thickness can have the effect of causing the layer thickness to oscillate. 
The layer appears to beoMe narrow when a nininun in atnospherio density 
occurs above the layer peak and a naxiaun occurs below. Conversely, when a 
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■axlauB ooours abov« the layer peak and a alnuaua ooours below the peakf the 
layer aeena to broaden. This last effeot nay partially explain the 
bifuroatlone often seen In eodiun Ildar data. The hl^eat peak aay oorreapond 
to a peak In ataosidterlo density foroed by a gravity wave. The sodiUB density 
■axlaui found below the layer aldpoint nay oorrespond to a slnlaua In 
ataospherlo density. Since the sodiun density perturbations due to the linear 
layer response and the ataospherlo density perturbations are 180 out of jrtwse 
at this pointy a peak In sodlua density would result. 

Teaporal variations In sodlua density are shown In Figure 5.14. Most of 
the density variations observed are consistent with very low frequenoy 
osolllatlons. The observation period y In this oassy was not long enough to 
cover one full period of the oscillation. The ATP was used to predict these 
low frequency osolllatlons. Near 92 and 97 losy higher frequenoy density 
oscillations are observed. These aay be associated with higher order teras In 
the layer density response. 

5.3 DATA COLLECTED WITH THE CANDELA LASER 

5.3.1 FEBRDARI 24-25y 1981 

A series of 97 spatial profiles were collected during this observation 
period. Five ainutes separated the starting tlaes of successive profiles. 
Because of the long observation period (roughly 8 hours) y the severe low 
frequenoy distortion in the ATP is confined to teaporal frequencies below 
0.004 aln~^. An exaalnatlon of the ATP shown In Figure 5.15 Indloates that the 
aaln peak at 0*0033 aln**^ Is on the edge of the region of distortion. The 
shot noise floor Is approxlaately 14 dB below the spectral peak and appears to 
doalnate the signal spectrua at frequencies greater than 0*022 aln*^. As a 
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Figure 5*15 Plot of the normalized average temporal periodogram of data 
collected on February 24-25 » 1961. 



Figure 5*16 Plot of the normalized average spatial periodogram of temporally 
filtered data collected on February 24-15t 1961. 
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r«sult» a tawporal f 11 tar outoff at 0.0267 was aalaotad. Tha ASP of tha 
taaporally filtarad data la ahown in Pigura 5.16. Tha ahot nolaa floor la 
38.3 dB balow tha apaotral paak and appaara to donlnata tha pariodograa at 
apatlal fraquanoiaa graatar than 0.4 ka**^. "nta apatial flltar outoff waa 
ohoaan to oolnolda with thla fraquanoy. 

TNO-diaanalonally filtarad apatlal profllaa ara oontalnad In Pigura 5.17. 
Thara la auoh activity below tha aaln layar paak. Thla paak la looatad naar 
94 ka at 2230 CST and by tha and of tha obaarvatlon period haa aovad down to 
90 ka. A ladga In tha aodlua danaity batwaan 82 and 86 ka la praaant at 2300. 
Thla ladga dlaappaara by 0030 and raappaara by 0130. Aftar 0300 tha ladga 
dlaappaara peraanantly. Ona of tha Intaraatlng faaturaa that later data alao 
ahara la a high altltuda Inoraaaa In aodlua danaity ooouring In tha 
pre-aunrlae houra. Beglnlng at 0200, tha aodlua danaity above 100 ka appaara 
to ateadlly Ineraaaa. 

The taaporal varlatlona in aodlua danaity ahown in Pigura 5.18 aaka It 
clear that tha Inoraaaa la not confined to hl^ altitudaa. Thla plot waa 
obtained by examining tha temporal varlatlona In calibrated data. An InoracMa 
In column abundance appaara aa a general Inoreaae In aodlua danaity at all 
altitudaa above 82.2 ka. Prevloua data could not be calibrated and, aa a 
reault, the temporal varlatlona In aodium danaity were obtained by aaaumlng 
the column abundance did not vary. In addition, theae data did not typically 
apan long perloda and did not extend Into the pre-aunrlae houra when column 
abundance Inoreaaea are obaerved. The enhanced aodium denaltlea aaaoolated 
with the early morning houra tend to maak wavelike featurea In the temporal 
plota. . To alleviate thla problem, the temporal varlatlona In aodium denalty 
have been plotted aaaunlng no varlatlona In column abundance oooured (figure 
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Figure 5*18 Plot of tbe temporal variationa in sodliaa density observed on 
February 24-25, 1981. The spatial and temporal filter cutoffs 
! Here 0.4 tali'*' and 0.0267 min'*', respeotively. 
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FEBRUARY 24-25, 1981 
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Figure 5<19 Plot of the teaporal varletione In sodlva density observed on 
February 24-2$, 1981 • The spatial and teeporal filters used In 
Figure 5.16 free used. These data have been noraallted to 
reaove the effeots of large ooluen abundance variation. The 
diagonal lines Indloate the apparent i^se progression. 
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5* 19)* This rmofs th» •ffsota of tho drautio InoFoaot la ooluan abundanoo 
and aakaa iiavallka faaturaa olaarar* The doalnant faaturaa ara apparantly 
relatad to a gravity nava with a vartloal wavalangth and taaporal period of 
11.5 In and 300 ■lnutes» respeotlvely. A horizontal wavelength of roughly 500 
km would be associated with auoh a wave. Reoall that the ATP contained a peak 
at 0.0033 aln~^. This corresponds to a period of 300 minutes. At a spatial 
frequency of 0.087 kn**^) corresponding to the wave's vartloal wavelmgthf a 
valley occurs In the ASP contained In Figure 5.16. This Is consistent with 
the predicted power spectrum presented In Chapter A. Hear 88 km a double 
frequency component with a temporal period of 150 minutes appears. In 
addition, a jAtase reversal In the features exhibiting a 300 minute period Is 
evident near 88 to 90 km. This probably corresponds to the phase reversal 
associated with the linear layer response. The double frequency oscillations 
become evident In this region because of the small size of the linear term. 

Figure 5.20(a) shows the average density profile (solid line) and the 
maxlaw and minimum densities observed during the evening (dashed lines). The 
layer peak at 92 km Is evident. Additionally, It Is seen that the density 
deviations below the layer peak are greater than those above the peak. Figure 
5.20(b) contains curves representing the maximal and minimum densities (dashed 
lines) predicted by a linear simulation of the layer response assunlng the 
average density profile of Figure 5.20(a) and gravity wave Induced 
perturbation In atmosi^erlo density of 65 peak to peak (4«0.03). The point of 
phase reversal Is observed at approximately 93 km. Figure 5.20(o) shows the 
result of a similar simulation In which the first nonllnesr term, k 2 from 
equation (4.36), Is Included. This second order simulation agrees more 
olosely with the exp«*lnental observations. It Is olsar that the seoond order 
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ENVELOPE 

OF SODIUM DENSITY 
PERTURBATIONS 
FEBRUARY 24-25. 1961 


ALTITUDE (km) 

Plsur« 5*20 (•) Bovelop* of sodliai dooaltjr voriatloM (dasbod llnaa) about 

tba atraraga danaltjr (aolld llna) obaanrad on Fabruarj 211-25. 
1981. Lifwar (b) and aaoond ordar (o) alaulatlona of tba 
anTalopa of danalty parturbatlons aasiaing tba aaa* avaraga 
danslty proflla and a grarlty uava anplltuda ooaalstany witb 3S 
paak varlationa In atnoapbarlo danslty ( d « 0.03) ara also 
plot tad. 



OP doublo fpoquonoy tors nkM slgnifloant oontrlbutlono noor tho layor poak 
•ad to • loaaor oxtont oox^r 82 ka. TIimo alBulatlon rosult* aro oonalatant 
with the tMporal profllaa oontainad in Pifupa 5>19* Tha doublo fpo^noy 
ooapoaonta obaorvod in this figura appaar at altitudaa uhara tha siaulation 
indioatas tha saoond ordar tara should aaka substantial ooatributiona* 

Tha tiaa variations in ooluan abundanoa ara shown in Pigurs 5*21* Baoausa 
of lasar alignaant problaas, astiaatas of tha absolute ooluan abundanoa oould 
not be aada* Tha ralativa ohangas in abundanoa oan ba datarainadf howavar* 
It is olaar that tha ooluui abundanoa aora than doubled bstwaan 2230 and 0500 
CST* Tha anhanoaaant prograssas rapidly after 0230* This was roughly 3*5 
hours bafora ground sunrise* Tha ragions in whioh tha largest abundanoa 
inoraasas ooour ara also dasoribad by Pigura 5*19* Piva kiloaatar slioas of 
tha layer ware axaainad and tha taaporal variations in tha ooluan oontant for 
asoh slioa ara plotted* It is olaar that tha aain inoraasas in ooliam oontant 
ara duo to anhanoad sodivat densities above 85 ka* Taaporal variations in 
ooluan oontant due to tha layer density response to wavs aotivity should ba 
•pparant. Howavar» those variations ara obsourad soaawhat by the aoming 
anhanoaaant* Thraa paaks in tha total ooluan oontant ara apparant at 2310» 
0150 and 0355 CST* Tha paaks at 2310 and 0355 ara saparatad by 285 ainutas* 
Tha faot that this period is olosa to tha 300 ainuta period astiaatad for the 
gravity wave indioatas that those paaks aay result froa wave aotivity (sao 
Chapter 4}* 

5*3*2 HMCH 11-12, 1981 

Oooasionally, lidar observations, suoh as those aada on Mar oh 11-12, 1981, 
Indioata a groat doal of aotivity in tbs sodiua layor that oannct olaarly ba 
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Figur* 5*21 Taaporftl variations in tha ralativa ooltan abundanoa obaarved on 
Fabniary 2b-25» 1981. Curvaa (a), (b), (o), (d) and (a) 
rapraaant tha total oolunn abundanoa and tha ooluan oontmt in 
tha 80-85» 85*90| 90-95# and 95-100 taa rancaa# raapaotlvaljr. 
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attribiittd to tiM lajror rooponot to o olaglo otaotphopio novo. On this 
ovonlag 60 opotlol proflloa wort ooUootod at 5 alauto intorvala. Tho kft, 
ahoNn Plfurt S>22» Indioatoa a groat doal of aotlvlty. ?oaka ooour at 
taoporal froquoaoioa of O.OO270f 0*011 » 0.0175 and 0*0262 oia*^. Booauao of 
tho longth of tbo oba«*vatlon porlod and tho offoota of taoporal windowinCf 
foaturoa at froquonoloa holoo 0.0059 aln”^ art dlstortod. A taoporal f 11 tor 
Mltb a outoff at O.031 oln*^ waa appllod to tho data. Tho ASP of tho 
taoporallj filtorod data la ahown In Plguro 5.23* Slnoo tho slfud apootruo 
aoooa to dlsappoar Into tho ahot nolao Itvol noar a apatlal froquonoj of 0*356 
ko*^t thla froquonoy waa aolootod aa tho outoff of tho apatlal flltor* 

TNo-dloonalonally filtorod apatlal prlflloa art oontalnod In Plguro 5*24. 
Ibo layor poiJc oooora noar 91*5 Is* It la dlffloult to dlaoom any auatalnod 
oovooant In tho layor'a foaturoa. Is. gonoralt howovoPt tho layor appoars to 
broadon aa tho ovonlag progroaaos. Thla growth la not oonflnod to tho uppor 
roaohoa of tho layor« howovor. Tho bottoooldo of tho layor bovoo froo 83 toi 
at 0025 to 79 Is at 0505 C8T* Tho largo donolty gradlont at 95 Is boooMO 
■ttoh oaallor by tho ond of tho oboorvatlon porlod aa a roault of tho upward 
oxpanalon of tho lajror* 

Taoporal varlatlone In aodiuo donolty aro plottod In Plguro 5.25* T^iOao 
ounroa havo boon nomallaod to rooovo tho offoota of Inoroaaoo In oolsm 
abundanoo and aako tho wavollko varlatlona In donolty oloaror. Thoro aro oany 
wavollko foaturoa proaont In thooo proflloa that nay roault froo tho prooonoo 
of oovoral wovoa. Aa a rooultt It i* ^mrf difficult to dotorolno tho 
paraoetoro of any olnglo wavo. Taoporal varlatlona la ooluon abundanoo aro 
plottod In Plguro 5.26. Aa laoroaoo la tho ooluon abundanoo by a factor of 
ooro than 2.5 la evident botwooa 0025 and 0505 C8T. Thla draoatlo Inoroaao 
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Figure 5*22 Plot of the normalized average temporal perlodogram of data 
collected on March 11-12| 1981. 
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Figure 5.23 Plot of the normalized average spatial perlodogram of temporally 
filtered data collected on March 11-12. 1981. 
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Figure 5*25 Plot of the teaporal variations in sodivn density observed on 

Haroh 11-12» 1981* The SMtlal and temporal filter outoffs were 
0*356 km~' and 0.031 nln'*^> respeotlvely. These data have been 
normalized to remove the effects of large variations in oolumn 
abundance . 
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Figure 5>26 Plote of the teaporal variations in ooliam abundanoe observed on 
March 11>12( 1981. Curves (a), (b), (o), (d), and (e) represent 
the total oolunn abundanoe and the ooluwi content in the 80-65 , 
85-90. 90-95t and 95-100 km ranges, respectively. 
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ooours priaarily In the 90-100 ka region. Approxiaetely 759 of the total 
inoreaae ia aooounted for by the increased ooluan content in this altitude 
range. 929 of the Inoreaae ooours between 90 and 95 ka and 339 ooours between 
95 and 100 ka. 

5.3.3 HARCH 12-13t 1981 

On the evening of March 12-13t ^981 103 spatial profiles spanning 8.58 
hours were oolleoted. The ATP for these data is presented in Figure 5.27* 
This ATP differs silently froa previous ATP's in that the spectral peak is 
broad and the shot noise floor is quite low (-22 dB). It appears that two 
overlapping spectral peaks are present. Even though the first speotral peak, 
at 0.0022 ain'^t is in the region of distortion, it is sli^tly larger than 
the peak at 0.0040 aln**^. Because of windowing, the speotral features below 
0.0039 min*'^ are distorted. This distortion, in general, deoreases the 
magnitude of the periodograa below 0.0039 aln**^. These peaks indloate the 
presence of long period waves (450 and 250 minutes, respectively). Althou^ 
the frequencies at which these two peaks appear suggest that the second peak 
nay be associated with the second harmonic of the 0.0023 ain~^ oscillations, 
an exaalnation of the filtered data indicates the presence of a separate wave 
with a frequency of 0.0040 iBln"^. A teaporal filter with a cutoff at 0.0308 
ain*^ was applied to the data. The ASP of the teaporally filtered data is 
shown in Figure 5.28. The low shot noise level (at -47.0 dB) results froa the 
large number of signal photons in each spatial photooount irofile and the 
effects of teaporal filtering. No distinct peaks appear in the ASP and the 
aagnltude of the ASP falls off gradually indicating a narrow sodiua layer. 
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Figure 5*27 Plot of the noraalized average tenporal perlodogrui of data 
oolleoted on Marob 12-13i 1981. 
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Flgur« S<28 Plot of the noraalized average spatial perlodograa of taaporally 
filtered data oolleoted on Maroh. 12-13* 1981. The t«aporal 
filter outoff Has at 0.0308 aln** « 
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A spatial filter with a outoff at O.A76 was applied to the data. Ihe 
resulting spatial profiles of sodlta density are shown In figure 5*29* As 
antlolpated, the sodlta layer remains narrow during most of the evening. 
Before 0300 CST Its e*^ width varies between 7 and 8 kilometers. After 0300 » 
the layer rapidly grows wider attaining a width of over 13 km after 0515. 
Peaks separated by 5.8 km are present below the layer peak. There Is a 
striking and oonslstent downward movement of these features at 0.4 m/s. The 
temporal period assoolated with this motion (240 minutes) corresponds to the 
spectral peak at 0.004 mln*^ In the ATP. The bottomslde of the layer also 
appears to osolllate up and down as these peaks sweep doimward. The 5.8 km 
spacing of the peaks Is lost near the main layer peak, nils probably results 
from the phase reversal in the layer response. The diagonal lines In Figure 
5.29 Indicate the apparent phase progression of the wave. The phase reversal 
appecurs to occur slightly above 90 km. Below this altitude, peaks In sodium 
density are aligned with the lines. Above this altitude, valleys In the 
sodlisB layer appear on the lines. These valleys are not as distinct as are 
the peaks below the point of phase reversal. This probably results from the 
fact that the density perturbations below the layer peak due to the linear 
layer response are. In general, larger than those above the layer peak. 

Temporal variations In sodlta density are plotted In Figure 5.30. The 
phase progression so evident In the plots of spatial profiles Is not clear 
beoause of the oonfused j table of wavellke features In this plot. This may 
result from the superposition of 2 dominant waves. A peak In the ATP also 
oooured at a very low frequenoy (0.0022 mln'b. Low frequency density 
osolllatlons may be contributing to this confusion. In addition, the narrow 
nature of the layer serves to Increase the magnitude of the nonllnearltles In 
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*29 TIm history of tbs sstlastsd sltltuds profllss of sodiiai 

density observed on Maroh 12«Ht 1981. The spetlsl and teaporal 
filter outoffs were 0.t76 UaT' and 0.0308 ftin~^» respeotlvely. 
Bstlaated profiles are plotted at 15 ainute intervals. The 
diagonal lines indloate the apparent phase progression. 
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Figure 5*30 Plot of the teaporal verlations in aodiue denalty obaerTed on 

Mar oh 12>p, 1961. The spatial and teaporal filter cutoffs were 
0.476 las**' and 0.0308 aln~', respeotlvely. Those data have been 
noreallsed to remove the effects of large variations In ooluan 
abundance. 



134 


the layer denalty response. A double frequency ooaponent with a temporal 
period of rou^ly 125 minutes Is definitely present near 87 end 94 la. The 
ooablnatlon of layer nonllnearltlea and the presence of two waves may explain 
the aany wavellke features apparent in this figure. 

Temporal variations In column abundanoe are shown In Figure 5.31* The 
peaks In abundanoe ooourlng at 2325 and 0345 CST may be Induced by the 250 
minute wave discussed previously. In addition to these peaks , a general 
Increase In abundanoe Is observed. The column abundanoe for the first 100 
minutes of the observation period was roughly 1.65 X 10^^ atoms/m^. This 
Inoreased by 431 to 2.65 X 10^^ atoms/m^ by the last 100 minutes of the 
observation period. As in previous evenings f most of this Inorease oooured 
above 90 km. 

5.3.4 MARCH 13-14, 1981 

The ATP for data oolleoted on this evening Is shown In Figure 5.32. 113 
spatial profiles oolleoted at 5 minute Intervals were oolleoted during the 9*5 
hours of this observation period. As a result of the length of this data set, 
the low frequency distortion In the ATP resulting from windowing Is confined 
to temporal frequencies less than 0.0034 min*^. A speotral peak near 0.002 
ain“^ Is evident Indicating the presence of long period osolllatlons In sodlrns 
density (with periods near 500 minutes) . A temporal filter with a cutoff at 
0.0308 aln"^ was applied to the photooount data. The ASP of the temporally 
filtered data Is shown In Figure 5.33. Beoause the shot noise level dominates 
the ASP at spat*<^l frequenoies greater than 0.4 km’S a spatial filter cutoff 
was selected at this point. 
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Figure 5*31 Plot of the teaporal verietlons in oolunn abundenoe observed on 
Maroh 12-13i 1981. Curves (a)» (b), (o), (d), and (e) represent 
the total ooltaui abundance and the ooliaan content in the 80>85t 
85-90t 90-95| and 95-100 ks rangeSf respeotively. 




Plsure 5*32 Plot of th« nonMllMd •vorag* twporal periodograa of data 
oollaotad on March 13-14, 1981* 
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Pigura 5n33 Plot of tba DorBaliaad‘‘avaraga spatial psriodograa of taaporallp 
flltorod data oollaotad on March. 13-14, 1961* Tba taaporal 
filtar cutoff Naa at 0.0308 aln~\ 
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Tht fllt«r«d spatial profilaa art plettad In Plgura S*34* TIm layar la 
not tanaralljr as narrow as It was on tha pravlows avanlnt* At 2025 CST tlia 
a**^ width of tha layar is 10 ka. It narrows until attaining a width of 6 la 
at 0030. Aftar 0100, it broadana, finally raaohing a width of 11.5 la at 
O5A0. A ragion in whioh tha dansity gradiant is larga novas downward fron its 
initial location naar 95 la batwaan 2200 and 0100 CST. Howavar, tha obvious 
downward novaamt of faaturas apparant on Maroh 12-13f 1961 is not prasant. 

Plots of tha tanporal variations in density ara shown in Pigura 5*35. Tha 
diagonal linos indioata tha apparant phase progression of a wave with s 
tanporal period of approxinataly 490 minutes and a vartioal wavelength of 
roughly 25 kn. A phase reversal ooours naar 95 ka (above the layar peak at 
92.5 kn). In addition, higher fraquanoy dansity osoillations ara prasant naar 
89 and 95 ka. Tha latter osoillations ooour naar tha point of phase reversal. 
Larga dansity gradients ooour at both of these altitudes. Sinoa larga 
gradients proaota nonlinaaritias, these high fraquanoy osoillations appear to 
be avidanea of tha nonlinear layar response. Tha 490 ainuta period of tha 
wave suggests ataoaphario tidal wave aotivity. 

The variations in oolunn abundanoa ara plotted in Plgura 5.36. A 505 
inoraasa in ooluan abundanoa begins to ooour two hours before ground sunrise. 
As in previous obeervatlons, the large inoreases in ooluan abundanoa ara 
assooiatad with a broadening of the sodlua layar and an inoraasa in ooltan 
oontant at high altitudes in the pra-sunrise hours. 

5.3*5 MARCH 16-17, 1981 

On the evening of Maroh 16-17| 1981, 126 spatial photooount profiles ware 
oollaoted over a period of 10.66 hours. Tha ATP of this data is shown in 
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Flfurt 5*34 TIm history of tbs ostlastod sltltudo profllos of sodlua 

dsnolty obsorvod on Hsrob 1961. Tbo SMtinl and tsnporal 

flltsr outoffo voro 0.60 loi~^ and 0.0308 aln* , raspaotivaly. 
Bstlaatsd profilaa ara i»Iottad ovary 15 ■Inutaa. 
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Figure 5*35 Plot of the temporal variations in sodium density observed on 

March 1981* The sMtial and temporal filter outoffs were 

0.40 10B~' and 0.0308 min'^i respeotively. The diagonal lines 
indicate the apparent phase progression. These data have been 
normalized to remove the effeots of large variations in column 
abundanoe. 
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Plot of the tenporal variations in ooluan abundanoe observed on 
Harob 13-lb, 1981. Curves (a), (b), (o), (d), and (e) represent 
the total oolvsm abundanoe and the ooluen content in the 80-85* 
85-90, 90-95* and 95-100 kB ranges, respeotlvely. 




Figure 5* 37* This ATP Is unusual beoause the shot noise floor is so far beloM 
the spectral peak (-17*5 dB). This is a oonsequenoe of the large nuaber of 
signal i^otons in eaoh temporal profile and the size of the tMporal 
variations in density. In addition » the ATP falls off into the shot noise 
level fairly rapidly. The temporal window used in computing the ATP results 
In distortion of spectral features below 0.003 min**^. The spectral peak occurs 
below this point at a frequency of 0.002 min*^. Beoause of the rapid fall off 
of the ATP with increasing frequency « a temporal filter cutoff at 0.025 
mln*^ was chosen. The ASP of temporally filtered data is shown in Figure 
5.38* Since the signal spectrum falls into the shot noise floor near a 
spatial frequency of 0.476 km**S a spatial filter outoff was selected at this 
frequency. 

Spatial profiles of two-dlmenslonally filtered data are plotted in Figure 
5.39. The layer is seen to evolve from a broad layer early in the evening 
into a narrow, sharply peaked layer at 0330 CST. By 0500, the layer is once 
more becoming broad. One of the major factors contributing to the formation 
of the narrow layer at 0330 is the decrease in amplitude of features below 90 
km at this time. The undulatory behavior of atmospheric layers discussed in 
section 4.4 is evident in this figure. A peak in sodium density near 95*3 km 
is observed at 2005 CST. This peak steadily moves downward reaching a minimum 
altitude of 91*2 km at 0150. By 0450, the peak has moved upward to 93*5 km. 
The period of this oscillation is approximately 11.5 hours. The exact 
movement of the layer peak la plotted in Figure 5.40. In this figure, the 
peaks and valleys for each spatial profile of sodium density collected during 
the evening are plotted. The movement of the layer i^k is clearly evident 
here. Since the period of this undulation is near 12 hours, this behavior la 
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Figure 5*37 Plot of the nomalized average temporal perlodogram of data 
oollected on March 16-I7t 1981. 



Figure 5.38 Plot of t^.j nonullzed average spatial perlodogram of tmaporally 
filtered data oolleoted on March 16-I7f 1981. The tsaporal 
filter outoff tias at O.0308 min**'. 
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Figure 5.39 Tiae history of the estiaated altitude profiles of sodiua 

density observed on Maroh 16-17» 1981. The spatial and teaporal 
filter cutoffs were 0.476 ka“’ and 0.0308 ain"\ respectively. 
Estiaated profiles are plotted every 15 ainutes. 
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suggaatlvs of ataoaphorio tidal aotlvlty* If this is the oassi the BOveMiit 
of the layer peak indioates a aaxiaum upward wind oooured near 1700 CST. The 
amplitude of this wind may be estimated by observing the maximiB displaoement 
of the layer peak. In this oaset lb is estimated that the maximum vertical 
wind was approximately 0.3 m/a. In addition to the undulatory motion of the 
layer peak, there are large density variations between 85 and 90 km. Many of 
the features in this region have peak to peak separations as small as 2 or 3 
km. The period of osolllation of these features is muoh less than the 690 
minute period assooiated with the movement of the main peak. This fine 
structure may be a layer response to gravity waves. Thus, the total layer 
response appears to be due to the presence of more than one atmospherlo wave. 
Since the layer response is, in part, nonlinear, there can be an Interaotlon 
between the responses to these various waves. This can cause density 
osolllatlons irtiose periods are not a hamorlo of any of the atmospheric waves. 
Such nonlinear interactions have not been treated in this thesis. 

Temporal variations in sodlm density are plotted in Figure 5.41. Density 
osolllatims with a temporal frequency of 290 minutes appear to be present. 
The vertioal and horizontal wavelengths assooiated with the wave apparently 
responsible for these osolllatlons are 18 km and 800 km, respeotlvely. These 
observations are tentative, however, because of the difficulty in clearly 
establishing the phase progression of the mve. It does appear that some form 
of wave activity is present in addition to the possible tidal motion disoussed 
previously. 

An increase in column abundance in the early morning hours was observed in 
these data as in previous observations. This is illustrated by Figure 5.42. 
The column abundanoe doubled from about 1.7 X 10^^ atoms/m^ to over 
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Figure S.41 Plot of the temporal variations in soditan density observed on 

March 1981. The spatial and temporal filter cutoffs were 

0.476 km~^ and 0.0308 mln*^, respectively. These data have been 
normalized to remove the effects of large variations in column 
abundance. 
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3*5 X 10^3 atOM/a^ during the 3 hours b«for« ground aunrino* As In prsvlous 
observations, aost of this Inorssss resulted froa an Inorsass In the ooluan 
content shove 90 ka* 

5.4 SOtMARY OP OBSERVATIONS 

It has been denonstrsted that aany of the features In the sodlua layer 
behave In a Banner oonslstent with the layer density response to ataos^erlo 
imves. The features assoolated with the linear layer response are often 
observed. These Include phase reversals In density fluotuatlons and larger 
density variations below the layer peak than above It. The periodic ohanges 
In thickness and the appearance and disappearance of bifurcations In the layer 
are probably assoolated with the layer response to ataospherlo waves. In 
addition, the aaplltude of ataospherlo waves required to Induce the observed 
sodlua density perturbations are reasonable* Evidence of nonllnearltles in 
the sodlua layer response Is evident and appears to cluster around the point 
of phase reversal in the linear response and around regions wnere large sodlua 
density gradients are located. Spectral features In both teaporal and spatial 
perlodograas are consistent with a layer response to gravity waves* 

An unexpected feature In these observations Is the early aornlng sodlia 
enhancement that was observed In every calibrated data set. This enhanceaent 
appears to be assoolated with the Influx of atoalo sodlua in the upper half of 
the sodlua layer. Yet, it is difficult to explain this enhanceaent with 
current theories of sodlua layer oheaistry. It aay be the result of aeteorlo 
deposition of sodlua. The morning enhanceaent which was also observed during 
our Ildar aeasureaents at the Goddard Optical Facility (Chapter 7) Is 
discussed further In Chapter 8. 
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6. 8TBDUBLB SODZUH LZDAB OBSBBVATZONS 

f 

t 

t 6.1 MBASORBMBNT TBCHNIQUB 

t 

Steerable lider meaeureaenta have been need by some groups to study 
horizontal density variations In the sodlua layer [Thoaaa et al.y 1976} 
Cleaesha et al.| I960]. These aeasureaents are laportant for several reasons. 
The exlstenoe of significant horizontal structure In the steady state sodlua 
layer (l.o. In the absence of ataosiRierlo wave aotlvlty) will be reflected In 
the nature of the layer response to wave aotlvlty. As a result* the 
horizontal structure of the layer aust be considered to aoourately predict the 
layer response. If su<^ horizontal structure exists* estlaates of the 
aagnltude and direction of the background wind aay be aade by observing the 
structure as It Is carried along by the wind. Additionally* the horizontal 
wavelength of waves which are Inducing density perturbations In the layer aay 
be deduoed froa steerable observations. Without steerable aeasureaents* this 
quantity must be Inferred froa estlaates of the vertical wavelength and 
frequency of the gravity wave obtained by examining data oolleoted at a fixed 
pointing angle. 

In most oases* these aeasureaents are aade with a steerable Ildar 



operating at a fixed el«‘vatlon angle and pointed at several azlauth angles. 
Typically* three azlauth angles separated by 120° are ohosen and spatial 
profiles of data are collected at each azlauth angle In a sequential aannar. 
As a result* every third spatial }>roflle la oolleoted at the saae azlauth 
angle. Large horizontal separations can be achieved by decreasing the 
elevation angle. Decreasing the elevation angle also Introduces larger losses 
due to ataospherlo soatterlng and absorption. Thus* acre laser pulses are 
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required to foneroto dat« proflloa of tqual quality as tho olovatloa anglo 
STOWS SBSller* This offsotlssly doorsasss tbs taaporsl rssolutlon of tbs data 
slnos tbs rats at wblob prof lies oan bs oollsotsd Is rsduosd* 

Tbs rsduotloQ In slsnal strsngth as a function of slsvatlon angls Is shown 
In Plfurs 6.1. for assoaphsrlo hoists » tbs vsrtloal ataospharlo tranaalsslon 
faotori T( la approxlaatsly 0.6 [Bhattaoharyya st al.» 1978]. Thus, at an 
slsvatlon ansis of 60°, It la antlolpatsd that only 64| of tbs snsrsy rsoelvsd 
wben tho syatsa la pointed at zenith will be available. At elevation anslea 
below 40*^, the effeots of ataospherlo attenuation beoom prohibitive. When 
three azlauth angles are soanned at an elevation angle of 60**, the tsaporal 
resolution la reduced to approxlaatsly 219 of the resolution which could bs 
aohleved with zenith asasuresents . This factor takes Into account both 
ataospherlo attenuation and the reduction In the profile oolleotlon rate 
resulting froa the scanning of the Ildar systea. Therefore, Insight Into 
horizontal variations In sodlia density oan be obtained by ssorlflolng 
tsaporal rssolutlon. 

6.2 PAST OBSERVATIONS 

The first steerable Ildar aeaaursaents of the sodlia layer were aade 
between Septsaber, 1975 and January, 1976 at Vlnkfleld, U. K. (51.4^N, 0.7^W) 
[Tboaas et al., 1976]. Soae of these observations Indloatsd obanges In the 
height distribution of sodlia ooouring over saall horizontal separations 
(approxlaatsly 15 ka) and short tlas periods (approxlaatsly 2 to 3 alnutes). 
These density fluotuatlons were attributed to perturbations Induoed by gravity 
wave activity. The raaalnlng observations, however, failed to Indloate any 
algnlfloant dlfferenoes over horizontal separations as great as 100 ka. 



Figure 6.1 Variation in received signal strength as a function of the 

elevation angle (a) and the vertical atnospherlc tranraiasion 
(7) to the altitude of interest. For aesoupherle heights and at 
a wavelength of 589.0 nm, 7 « 0.6 [Bhattacharyya et al.» 1978]. 
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Subs«qumit ■•Muraasnts by th« group (ooll«ot«d in August t 1976) l«d to 
tiM oonoluslon that th« differonoo* in denaity observod at u<r«w points 
'jsp^atsd by 80 loi wnrs duo tc horlaontal variations in ths steady stats 
sodiun layer ratber than tine variations in the density of a horiaontally 
unifom layer (Thonas et al«i 1977]* It should be noted, houever, that in the 
majority of the data it was difficult to detect any systeaatio density 
variations and therefore no oonoluslons were reached. In addition, long tern 
variations in sodiia density oould not be observed sinoe the data sets did not 
span sore than 2 hours. 

Here extended observations have been made In the southern henlsphere by 
Clenesha et al. [1961] at Sao Jose dos Canpos, Brasil (23^3, A5°W). On the 
evening of July 21, 1979» ten hours of scanned data were oolleoted at an 
elevation angle of 69.1^ and at three azinuth angles. The horizontal 
separation of points achieved in this manner was roughly 60 Iob. A 
oros8>oorrelatlon analysis was perfomed for data oolleoted at heights between 
82 and 99 kn. Velocities were calculated fron tias«shlfts derived In this 
nanner. These velocities were found to vary in an osoillatory nanner with 
amplitudes increasing with height. The east-west component of wind showed a 
vertical wavelength of approximately 10 km and the north-south ooaponent 
showed a vertical wavelength of 5 km. Wind velocities were estimated to vary 
between 6 and 146 m/s. 

On August 25-26, 1979» using the ssme steerable system, this group 
observed s dramatic enhancement in sodius density [Clamesha et al., 1980]. 
This enhancement, associated with a tenfold increase in sodium density at 95 
kn, was restricted to s narrow time interval and a height interval of little 
more than 2 km. It was concluded that this enhaneaftsnt was fw^bably due to 
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the Influx of Mtoorlo aaterial. The tiae lage assoolated with the passage of 

this patoh of enhanoed aodlta suggested a veloolty of 205 n/a at a heading of 

0 

201. In addition* It vraa estlaated that the extent of the patch In the 
direction of motion waa on the order of 100 ka. The horizontal extent of the 
patch perpendicular to the direction of motion waa oonalderably 0 *eater. 
Deducing the motion and shape of this patoh illustrates one of the problems 
inherent in estimating velocities from this type of data. Both the shape and 
velocity of the patch affect the density variations observed at widely spaced 
points. A simple time lag approach is not sufficient to deduce wind 
velocities since the shape of the patch also determines the nature of the 
observed density variations. 

6.3 EXPERIMENTAL RESULTS AT GSPC 

In May and June, 1981, a 48 inch telescope located at the Goddard Space 
Plight Center Optical Test Site was usea In conjunction with components of the 
University of Illinois sodium Ildar systms to make steerable lidar 
observations of the mesospheric sodiiai layer. The telescope facility, located 
at 39.1*^N, 76.8Plf, consists of a 48 inch Cassegrain astronomical telescope 
with Coudd focus and supporting equipment which enables this instrument to be 
used in a variety of laser ranging experiments. The telescope aperture area 
is approximately the same as the objective area of the receiving telescope 
used in Illinois. However, as well as making steerable measurements possible 
because of the sophisticated telescope mount, the high quality optics permit a 
small receiver field of view. This serves to reduce the noise counts due to 
background sky illumlnaticn. Photographs of the telescope are contained In 
Figures 6.2 and 6.3. 
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A photograph of the telescope facility Is shorn in Figure 6.2* The 
telescope rests under the dome at the top of this photograph and aooess to the 
CoudI focus is on the ground floor. mie Inteferenoe filter and PMT are 
positioned near the fooal point of the telescope. No oolllaatlng lens is 
required due to the long fooal length of the primary/seoondary mirror 
combination in the telescope. The aperture of the focusing lens in the PMT 
cooled housing limits the telescope field of view to less than 1 mrad. This 
field of view can be reduced further by inserting a field stop iris in front 
of the PMT. Figure 6.3 contains a photograph of the telescope body showing 
the location of the laser head. The head was mowited on the telescope to 
insure that the laser beam tracked the telescope motions accurately* 
Coincidence of the transmitted laser beam and the telescope field of view was 
insured by adjusting the folding mirror (seen to the right of the laser head 
in Figure 6.3). Because the laser beam divergence angle was larger than the 
telescope field of view, a beam expander was inserted in the optical path 
between the laser and the folding mirror. This served to reduce the beam 
divergence to roughly 0.1 mrad. Laser tuning was monitored by examining the 
scattering of the laser output in a grating spectrometer (for rou^ tuning) 
and in a sodium cell (for fine tuning to the sodium D2 line). A fiber optic 
cable was used to couple the laser output to these instruments* 

The position of the telescope was computer controlled and therefore the 
telescope was easily and rapidly positioned at the desired azimuth and 
elevation angles, nie range of permissible angles was rather restricted) 
however) by mechanical considerations related to the laser head. The cables 
and hoses connecting the laser head to the dye reservoir and power supply were 
only 15 feet long* Since these last two items remained stationary) the 
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azimuth angle oould only be varied by 70 or 80 . In addition) the elevation 
angle oould only be varied from 0° to 90^ beoauae of the teleaoope mount 
design* It was not possible to move the telescope past zenith. Steerable 
measurements were Initially attoipted on the evening of June 8-9 » 1981. nie 
telescope %ns scanned between three points; point A (azlmuths277'' ) 
elevatlons90° } ) point B (azlmuths277'^ t elevations70*^ ) t and point C 
(azimuths 307° ) elevatlons70° ) . A greater horizontal separation was achieved 
In later observations of the sodlm layer. In order to examine 3 points which 

were equally spaced horizontally) the observations alternated between zenith 

0 

measurements and measurements at an elevation angle of 60 idille the azimuth 
angle was varied by 60° . Observations of sodium density at the Increased 
horizontal spacing were made on June 17»18) 1981. The locations of the points 
at which sodium observations were made on both evenings are shown in Figure 
6.4. Within the mechanical constraints mentioned above) an attempt was made 
to select points over sparsely populated regions In order to avoid areas idiere 
sodium vapor llj^ts were used in large numbers. It should be noted that 
because data was gathered at depressed elevation angles ) the horizontal 
separation of the observed points was a function of altitude. For the data 
collected on June 17-18) 1981) the observation points were separated by 50 km 
at an altitude of 100 km whereas at an altitude of 80 km the separation was 
only 40 km. 

Nine spatial profiles of ^otooount data were collected on June 8-9) 1981. 
Each profile consisted of the Integrated returns from 500 laser pulses. 
Because only three photocount profiles were collected at each pointing angle) 
It was not possible to temporally filter the data. Receiver range gate and 
delay times of 3 pa and 60 pS) respectively) were utilized at all three 


*CC>Q 



Figure 6.4 Horizontal spacing of observation points corresponding to an 
altitude of 90 km for measurements made on June 8-9, 1981 
(points A, B, and C) and June 17-18, 1981 (points 1, 2, and 3). 
The elevation angle (a) and azimuth angle (B) are given for each 
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pointing angles. The ASP of the three profiles oolleoted at zenith Is plotted 
In Figure 6*5. The large peaks at 0.12 km*^t 0.17 and 0.28 km'^ 
Indicate that wavellke features are present. These data were spatially 
prooeased In the manner described In Chapter 3* A spatial filter with a 
cutoff at 0.4 km**^ was applied to the zenith data. An effective cutoff of 0.4 
km~^ for vertical variations In sodium density can be achieved for the 
remaining profiles by taking Into account the elevation angle of off>zenlth 
observations. As a result » a spatial filter with a cutoff at 0.376 was 
applied to the data collected at a zenith angle of 30°. 

Plots of the estimated sodium profiles for each pointing angle are shown 
In Figure 6.6| 6.7 and 6.8. A narrow^ sharply defined peak In sodium density 
at 95 km Is apparent In all three figures. This peak is approximately 2.5>3>0 
km thick. The main layer peak appears near 90 km. Additionally, a peak In 
density near 102 km evolves during the 90 minutes covered by these 
observations. It is unusual to observe distinct features at such high 
altitudes. Above the main peak, layer features are consistent at each 
pointing angle* Below the layer peak, however, features are not so well 
behaved. Features below 90 km In the profiles oolleoted at points B and C are 
similar whereas the profiles collected at point A differ more markedly from 
these* To Illustrate, at 0346 EDST (point B) and 0356 EDST (point C), peaks 
In sodium density at 84 and 87 km and a valley at 87*5 km are observed. At 
0406 EDST (point A), a peak at 84 km Is apparent but no valley Is observed at 
87*5 km. The similarity In profiles at B and C is not surprising since B and 
C are closer to each other than to point A. The variability of the layer 
below 90 km Indicates that variations In sodlia density with small horizontal 
scales (perhaps as «Dall as 30 km) are present. Layer features also appear to 
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.6 Altitude profiles of the estimated sodium density observed at 
point A (a s 90^) on June 8-9) 1981. The spatial filter cutoff 
was at 0*4 km' * 
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JUNE 8-9, 1981 



Figure 6.7 Altitude profiles of the estimated sodium density otMerved at 
point B (a s 70^, gs 277°) on June 8«9» 1961. The spatial 
filter cutoff was chosen to give an effective vertical cutoff at 
0.4 km*', as discussed In the text. 
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JUNE 8-9« 1981 



.8 Altitude profiles of the estiaated sodita density observed at 
point C (o « 70°, B* 307°) on Juiia 8-9, 1981. The spatial 
filter cutoff was ohosen to give an effective vertical outoff at 
0.4 km” . as discussed in the text. 
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ohang* algnlfioantly tiithin a 60 mlnuta pariod» Dnfortunatalyi thaaa 
obaarvatlona apan only 1*5 hours and it la not posslbla to daduoe details 
oonoernlng tha horizontal atruoture of the layer or the direotion of 
propagation of ataospherlo waves whloh may be present. 

Beoauae of poor weather, additional data could not be oolleoted until June 
17*18, 1961. On this evening, a series of 8l spatial profiles were oolleoted* 
Reoall that three points, separated horizontally by 45 km at an altitude of 90 
km, were observed (Figure 6.4). Each spatial profile consisted of the 
integrated returns from 500 laser pulses. A 3 ks range gate and a 196 PS 
delay factor were used in obtaining this data. The delay factor was chosen to 
insure that Rayleigh scattering from 30 km as well as sodium returns lould be 
received as the elevation angle was varied frcm 90° to 60°. The integration 
time for each profile was approximately 2 minutes and a new ^*ofile of data 
was started every 4 minutes* Thus, 12 minutes separated data profiles 
oolleoted at the same pointing angle* Because of vser malfunction, data 
was not oclleoted between 0120 and 0200 EDST. For the purpose of temporal 
filtering, an interpolation procedure was utilized to provide estimates of the 
missing profiles. The temporal and spatial processing techniques described in 
Chapter 3 were adapted to process the off-zenith data* 

The average temporal periodogram for zenith measurements is shown in 
Figure 6.9* To insure that no temporal variations in sodloa density were 
inadvertantly removed, a temporal filter cutoff at 0.033 mln*^ was chosen and 
applied to the data oolleoted at all three pointing angles* The average 
spatial periodograms of the temporally filtered data are shown in Figures 
6.10, 6*11 and 6.12. For off-zenith measurements , spatial frequencies have 
been scaled to reflect spatial photooount variations in the vertical rather 
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9 Plot of tb« noraallsod avortgo tesporal parlodograa of sonlth 
data oollaotad on June 17*16| 1961* 
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Figure 6.12 Plot of the normalized average spatial periodogram of temporally 
filtered data collected on June 17-18, 1S8l at point 3. The 
temporal filter cutoff was at 0*033 mln**^. 
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than the slant range* As a result > the off-zenith perlodograms may be 
compeu*ed directly with the perlodogram of zenith photocount data. A spatial 
filter with a cutoff at 0.4 km**^ was applied to the temporally filtered zenith 
data and a spatial filter with a cutoff at 0.345 km'"^ was applied to the 
off- zenith data. Taking Into account the elevation angle of the off zenith 
data, thi*-' laist filter gives an effective cutoff at 0.4 km~^ In the vertical 
direction. Thus, wavelike features with vertical wavelengths of 2.5 km or 
less will be removed from all three sets of data while features with longer 
wavelengths will be retained. 

The two-dlmenslonally filtered data Is presented In Figures 6.13y 6.14 and 
6.15. It is evident that the data collected at each pointing angle are very 
similar. The most significant difference Is the presence of a pronounced peak 
near 95 km at the start of the observation period In data sets 1 and 2 and Its 
absence In data set 3* A great deal of activity Is evident In all three data 
sets. In each case, the peak located near 98 km Is observed to move downward, 
reaching 95 km by 0430* In addition, the sodium density associated with this 
peak approaches the density of the main peak at times. This behavior is 
observed in all three data sets. All three data sets also contain very leu*ge 
density gradients at the bottomslde of the layer. In some cases, the sodlm 
density Increases from a very small to a very large value In 0.9 km. The 
spatial profile In Figure 6.16 contains such a feature. Sodxum densities have 
been calculated from temporally filtered ohotocount data and two-dlmenslonally 
filtered data* In each case, the sodium density Is observed to Increase from 
0.07 of Its peak value at 82*35 km to 0.63 of Its peak value at 83*25 km. As 
a result, the vertical sodium density gradient in this region Is approximately 
6.5 X 10^ Na atoms/m**. 
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Flgijre 6*13 Altitude profllea^of eodlia density observed on June 17<-18, 1961 
St point 1 (o| s90 ). The spatial and teaporal filter outoffs 
Mere 0.4 loa**' and 0.033 aln* t respectively. Interpolated data 
are Indicated by dashed lines. 
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Plot of the spatial profile of density observed at 23^8 EDST at 
position 1. The circles are range corrected photocount data and 
the smooth curve is two-dlmensionally filtered photooount data. 
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Figure 6.17 shows the temporal variations in sodium density assoolated 
with data set 1 (the zenith measurement). These variations are plotted In the 
same fashion as the temporal profiles presented In Chapter 5. It Is difficult 
to dlsoern any evidenoe of a systematlo movement of features. No phase 
progression is clearly evident anu therefore It Is difficult to conclude that 
wave activity is causing the density variations observed In these data sets. 
By plotting the locations of peaks and valleys in the layer , It should be 
possible to determine the direction of propagation and horizontal wavelength 
assoolated with waves Inducing density perturbations in the layer. Figure 
6.18 contains such a plot. Unfortunately • it is not clear that a single trave 
Is responsible for the observed motions of the layer's peaks and valleys. If 
a single wave were responsible for these density variations (assuming there 
were no horizontal structure to the steady state sodium layer)) the curves In 
this figure would differ from each other only by a phase delay. However) this 
Is not the case. The observed activity could be the result of the 
Interactions of several waves or It could Indicate that a sodium layer with 
horizontal structure was being swept overhead by the mesospheric nfiun winds. 
The motion of the layer peak near 88 km is consistent with the motion of a 
patch of enhanced sodium density moving at approximately 90 m/s towards the 
southwest (at a heading of 220°). The extent of the patch In the direction of 
motion Is approximately 300 km. It Is not possible to accurately determine 
the extent of the patch In the direction perpendicular to the direction of 
motion. These deductions were made by examining the upward movement of the 
layer peak near 0230 EDST. The motion of the layer peak which) at 2330 SDST) 
Is located at 98 km does not Indicate the presence of a region of enhanced 
sodium density sharing the characteristics listed above. In fact) the peak 
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Figure 6.18 





observed at position 1 la consistently above the peaks observed at positions 2 
and 3> The conjectured region of enhanced sodlun density responsible for the 
motion of the layer peak at 88 kn does not appear to extend to 98 Ion. A more 
oomplete explanation of the observed density variations would Include the 
effects of wave activity and the effects of a mean wind on a horizontally 
structured layer. 

Figure 6.19 shows the time variations In column abundance for each 
pointing angle. The spatial profiles of sodium density used wo c:.lculate the 
column abundance have been spatially , but not temporally > filtered. The same 
spatial filter cutoffs used to obtain the profiles In Figures 6*13) 6.14 and 
6.15 have been used. The laser oialf unction which interrupted data collection 
at 0120 also rendered the laser tuning suspect. It Is possible that the 
spectral output llnewldth of the laser was Increased as a result of optical 
misalignment In the laser cavity. This would have the effect of decreasing 
the effective sodivn backscatter cross section. As a result) the estimates of 
column abundance for data collected after 0200 may be too low. Relative 
In these estimates after the interruption In data) however, will remain 
accurate. An examination of Figure 6.19 seems to confirm this suspicion. The 
calculated column abundances Indicate a sudden drop between 0120 and 0200. 
After 0200, however, a steady rise In column abundance at all three pointing 
angles Is apparent. This Is consistent with data collected at Urbana, 
Illinois which Is presented In Chapter 5 Indicating a morning enhancement In 
sodium density. As Indicated in these data sets, this enhancement begins near 
0300 (about 3 hours before ground sunrise). The differences observed In the 
column abundance measurements after 0200 between each pointing angle may 
Indicate that there was significant horizontal structure to the sodium layer* 
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A general eaattfard Increase In ooluan abundance la observed* The horizontal 
gradient In oolwn abundance calculated at O^IOO EDST Indicates the BaxlBin 
Increase occurs at a heading of 110^ (Just south of due east}* This Is 
Interesting since the region of hipest meteor deposition occurs In this 
general direction at 0400* The region of maximum meteor deposition rate la 
Inferred from meteor radar observations* These observations typically 
indicate a maximum meteor radar echo rate occurs near 0600 local time 
[McKl-iley, 1961]* The differences In column abundance at each pointing angle 
become smaller as the sodiin density enhanoement prepresses and sunrise 
approaches* 

Before 0130, the column abundances at each pointing angle were roughly 
equal. The sharp decrease In abundance near 0000, observed at all three 
points, is probably an artifact that arose as a oonsequence of a drift in 
laser tuning. Laser tuning was checked after each zenith photocount profile 
was collected to oilnialze problems of this type* After 0200, the abundanoe 
variations at points 2 and 3 closely track each other. A sudden increase lu 
abundance at 0300 EDST Is seen at these two points but not at point 1. No 
periodic varlatioiis in colimin abundanoe which could olearly be associated with 
atmospheric waves are evident. 


7. LIDAR OBSERVATIONS OP STRATOSPHERIC AEROSOLS OVER ORBANA, ILLINOIS 


7.1 INTRODUCTION 

A Mjor «ruption destroyed the suamlt of Mt. St. Helens in southern 
Hsshington (46.20 Nt 122.16 tf) and projeoted massive amounts of voloanlo ash 
and dust into the stratosphere on May 18, I960. Prevallinf winds oarried the 
dust clouds eastward over muoh of the United States. The vertioal plume of 
dust and ash initially rose to more than 19 km above sea level (as measured by 
Portland airport's weather radar). The subsequent advanoe of olouds of dust 
at 3kmt 9kat I6km and 16 km altitudes are shown in Figure 7.1 CSEAN Bulletint 
I960]* Some of these olouds reached central Illinois by noon on May 20. On 
May 25f the dust olouds reached Europe [Meixner et al., 1981} Hauoheoorne et 
al.f I960; Aokeraan et al.| I960). Initiallyi layers of dust were located at 
low altitudes (near 13 km). Near July lOf however, stratospheric layers near 
21 km were observed over France [Lefrere et al., 19&1] and Italy [D'Altorio et 
al., 1981]. The initial dust olouds passed Japan on May 27 [Iwasaka and 
Hayashida, I960] and were observed over Hawaii on July 7 [Deluiai et al., 

1980]. 

In order to take advantage of this geophysical event, the Orbana lidar was 
reconfigured to observe scattering from stratospheric heights* Itoasurements 
of the height distribution of the volcanic dust began on Nay 20^^. Receiver 
range gate and delay factors of 2 {is and 4{is, respectively, were ooomonly 
utilized in collecting the data presented here. The laser output wavelength 
was 589*0 nm* The combined effects of receiver ran^ gate and laser pulse 
shape limited the altitude resolution to approximately 750 meters. Receiver 
aperture and laser mergy were reduced from the values listed in Tables 2.1 
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Figure 7*1 Paths traveled over North Amerloa by clouds of volcanic ash at 3 
kBy 9 km» 16 km, and 18 km altitude. Tick marks along each line 
show positions of ash cloud fronts at that altitude every 24 
hours (at noon GMT). The date at each tick mark Is indicated. 
[Sean Bulletin, 1980]. 
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and 2*2 to prevent overloading of the PUT froa the atrong low altitude 
Rayleigh returns. Because the laser and reoeivlng telescope were separated by 
10 deters ( their fields-of-vlew did not completely overlap for altitudes below 
6 ka. Individual profiles were obtained by firing the laser and aoouaulatlng 
the returns froa each laser pulse over periods ranging froa 2 to 10 minutes. 

7.2 LIDAR MEASUREMENTS OF THE BACKSCATTER RATIO 

Lidar aerosol measurements are usually made by utilizing elastic 
backscattering In which the received signal is at the same wavelength as the 
transmitted laser pulse. The lidar equation (equation 2.1) may be used to 
express the relationship between the expected photooount and the voluae 
backscatter coefficient. 


<ar(3)> * 





z 

Ta 

2 

exp 

-2 I (&a(s) 

J 


J 

0 


(7.1) 


tdiere 


<x> = expected photons detected in the Interval (t-2*, t) 
n = overall system efficiency 

2 

s receiver aperture area (a ) 
hv s energy of one photon at frequency v (J) 

s photocount rate due to background and dark noise (s'b 
P X laser pulse energy (J) 

8 X voluae backscatter coefficient (a“^) 
a X ataospherlc extinction coefficient (a**^) 

Q X speed of lli^t (a/s) 
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T 8 reoelver gate tlM (a) 


Because elastic backscatter Is considered rather than resonant scattering, the 
llfetimj of the excited state Included In equation (2.1) Is no longer of 
concern. In addition, we have Ignored the effects of laser pulse shape on the 
received photocounts. 

Both the volune backscatter and atmospheric extinction coefficients are 
composed of molecular (Rayleigh) and aerosol (particulate) components. 




8(a) » 6^ (a) + 3^ (a) 

(7.2) 



o(a) - o^(s) + o^(s) 

(7.3) 

It Is often 
backscatter 

convenient to 
ratio, i?(s) . 

analyze Ildar aerosol data In 
This Is defined as the ratio 

terms of the 
of the total 


backscatter coefficient to the Raylelj^i backscatter coefficient. 


i?(2) 


3^(2) + 3^ (a) 6^ (a) 

3^(a) “ ^ ^ 6^(2) 


(7.4) 


A common means of determining /;( 2 ) exists which depends upon scattering from 

an altitude at which 6 and the losses due to atmospheric attenuation are known 

[McCormick, 1975; Russel et al., 19791* This altitude Is known as the 

* 

normalization hel^t, a . The losses due to atmospheric attenuation were 
accounted for In equation (7.1) by the factor 

- cxp[-2x 1 (7.5) 

A. 

Where Q(a)ls actually the atmosi^erlc transmissivity and t Is the optical 

z 

* 

thlokneas of the atmosphere between 0 and z • Quite often t z is In excess of 
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20 km or 30 km* At these high altitudes most of the baoksoatter oan be 
attributed to Rayleigh baoksoatter since few aerosols are present. R(a) oan 
be derived from the lidar measurement as follows; 


Hz) 






(7.6) 


where <x^(s)»><a;(s)>-X^T. Q(z) is usually provided by an atmospherio model 

which may be updated using Ildar data.. 6 (a) is determined from radiosonde or 

R 

satellite data or fr<»i an atmosiriierlo model. A scattering ratio of 1 
indicates that few aerosols are present and large scattering ratios generally 
indicate large aerosol densities. The aerosol backscatter coefficient) B^(z) , 
can be determined from R(z) and Bjf(z ) . If the scattering properties of the 
aerosol are also known, the aerosol density oan be computed once ^J^(z) is 
determined . 

If Rayleigh scattering is present from an altitude where there are few 

aerosols, computing i?(s) is stralghtfoniard. Since the atmospheric aerosol 

* 


content is not known beforehand and niust be known at ^ , an initial 
choice of the normalization altitude above regions where aerosols often reside 
is made. For this reason, 30 km is typically used. Using an atmospherio 
model, ^(z) can be estimated. For these measurements, values for Q(z) were 
assumed based upon direct measurements of aerosol oonoentrations and 
properties frcm 5 to 28 km by Rosen [1966] in a series of six balloon flights. 

6 can be estimated if the baoksoatter cross section and the atmospherio 

R 

density are known. Radiosonde measurements often supply atmospherio ;»‘essure 
and temperature as a function of altitude for heights less than 30-35 km. 
Using the ideal gas law and these data, the atmospheric density oan be 
approximated. 
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P(«) 

P(a) - 


(7.7) 


iAior« p s number density of the ateosj^ere (aoIe*a*3) 
p a ataospherlo presaure (N*a*^) 

T a ataoai^erlo teaperature (°K) 

P* a gas oonstant (8*3143 J*aole~^*K**b 

If radiosonde data are not available up to the desired height) an 
exponentially decreasing density Is assumed v/hidti aatohes the density computed 
fron the highest available radiosonde measurements* The atmospheric scale 
helj^t Is used to determine the rate of decrease of the assumed density 
profile* The volume backsoatter coefficient) 8p(x) ) Is given by 


e^(a) - p(a) (7.8) 

where Is the backsoatter cross section. Traperature effects on the cross 
section are Ignored with little resulting error* Combining Q(z) t $i,(x) 

A 

the photooount datS) P(a)can be computed according to equation (7*6)* 

* 

Hie Initial choice of z must sometimes be altered since the aerosol 

content of the atmosphere Is not known a priori* 6^ and 8^ are non-negative 

so It Is clear from equation (7*4) thatP cannot be less than 1* At a*) si 

because It was assumed that no aerosols were present at this altitude* As a 

* 

result) an Iterative procedure may be required to select a ) compute i?(a) , and 
Insure that the computed i^lues of i?(a) are not consistently less than one* 
In effect) a search for an aerosol-free altitude takes place* In most oases ) 
scattering from 30 km la due almost entirely to molecular scattering and no 
Iterations are necessary* 
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7*3 EXPERIMBMTAL DATA 

The first clouds of voloanio dust assooisted %rltb the eruption of Nt« St* 
Helena passed over Urbansy Illinois at approxiaately noon on May 20 y I960* A 
prelininary report on lidar observations of these olouds is contained in 
Ctardner et al.y [1980]* Figure 7*2(a) shows the scattering ratio computed 
from the average of 5 photocount profiles collected between 0004 and 0034 COST 
on May 21. Four distinct layers at 7*5 kmy 9*9 ksy 12 km and 13 km are 
apparent* The scattering ratios associated with these peaks are 2*5y 5*5y 33 
and 36y respectively. The layers at 7.5 km and 9*9 km nay be associated with 
particulate matter from a large forest fire which was burning in northwestern 
Ontau?lo at this time* The two peaks near 12*5 km actually appear to be one 
bifurcated layer* Radiosonde measurements obtained over Sal me y Illinois y 
approximately 185 km south of Urbansy at 1800 COST on Nay 20 were used to 
compute atmospheric density andy thusy * Since measurements such as 
these are routinely taken at Salem and are readily available y radiosonde data 
is used to deduce for each data set presented in this Chapter* Because the 
density associated with these dust clouds was much greater than the aerosol 
concentration included in the model provided by Rosen y Q^(z) was 
underestimated* This generally results in estimates ot R(z) trtiich are too 
low* In addltiony the dust clouds on the evening of May 20-21 were so thick 
and dense that the transmitted laser pulse failed to penetrate than* As a 
consequence y scattering ratios above 14 km cannot be estimated since no 
backsoattered radiation was received from above this altitude* The scattering 
ratios may be seriously underestimated above the regions where losses due to 
scattering and absorbtion frcm the dust clouds was great* However y other 
lidar measurements of i?(s) made over Illinois on the same evening tend to 
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MAY 21, I960 
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MAY 21, 1980 
0105-0117 COST 



SCATTERING RATIO 


Scattering ratio profiles derived from Urbane lidar data 
oolleoted on May 20-21, 1980. 
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oonfira the eoattering ratios plotted here. During the late evening of May 21 
and early morning of May 22 an airborne aerosol Ildar operated by NASA*Langley 
made measurements approximately 75 miles south of Urbans [MoCormlok et al., 
1980]. At 2325 COST) aerosol layers at 13*81 and 12*61 km assoolated with 
soatterlng ratios of 71 and 3*5 were observed (Figure 7*3)* These data, 
oollected with a ruby laser operating at a wavelength of 694.3 nm, uompare 
favorably with the Urbana Ildar measureiBents. 

The soatterlng ratios observed later during the same evening, between 0105 
and 0117 (3)ST, are plotted In Figure 7.2(b). The peaks below 11 km have not 
changed In the span of 30 minutes separating the plots In Figures 7.2(a) and 
(b). The peaks associated with the bifurcated layer near 12.5 km have changed 
In magnitude, however. The maxlmuD soatterlng ratio, observed at an altitude 
of 13 km, is now 40. It appears that this large layer may have some 
horizontal structure assoolated with it. 

Balloon measurements of temperature and wind speed and direction show an 
Interesting correlation with i?(z) . Figure 7*4 shows the wind speed and 
direction measured over Salem at 1800 COST on May 20, I98O* Note the marked 
change In wind direction above 10*5 km, shifting from 270 ^ at 10.8 km to 238^ 
at 13 km. It appears that the main bifurcated layer between 11 and 14 km was 
associated with a southwesterly wind component over southern Illinois. The 
bifurcated layer also occurs at the same altitude as a wind speed minimum (20 
knots at 13 km). Figures 7.5 and 7.6 show the temperature profiles measured 
over Salem on the evenings of May 20 and May 21, I960, respectively. While a 
temperature minimum and a maximum are present near 11 km and 12.5 km, 
respectively, on May 20th, there Is only an Inflection present on May 21st. 
The bifurcated layer near 12.5 km Is no longer present on May 21st. This Is 
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Figure 7.3 



SCATTERING RATIO 


The profile of scattering ratio observed fron an airborne Ildar 
75 miles south of Urbans. These data were collected at 2325 
CDST on May 21, 1980 CMcCormlck et al.. 1980l . 
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illustrated in Pifurs 7.7(a)* 3avantaan profilaa wart oollaotad during tha 
avaning and in aaob oaaa tha aaroaol layara wara panatratad bp tha laaar pulaa 
allowing tha noraalisation haight to ba ohoaan abova tha lapara of voloanio 
dust, k aingla distinot layar at 13.5 la ia praaant aaaooiatad with a •axiaia 
aoattaring ratio of 27. At tiaaa during tha obaarvation pariod thia proainant 
layar did ahow avidanoa of a ainor bifurcation. Howavar, at no tlaa waa thia 
bifurcation as pronounoad as that obsarvad on Nay 20th. by ooaparing tha 
raoaivad signal strangtha at 12 ka and 15 IoB| tha optical thicknaas of tha 
layer at 13.5 ka waa calculatad to ba 0.82. Tha layara avidant balow 11 la on 
tha pravioua avaning ara no longer diatinot. 

Poor waathar pravantad additional aeaauraaanta until tha avming of May 
25>26| i960. The aoattaring ratloa obaarvad on thia avaning ara plotted in 
Figure 7.7(b). Onoa againi a bifuroatad layar oantarad at 12.5 iai ia evident. 
Howavar y tha aa x i a ua aoatterlng ratio has daoraasad by an order of aagnituda 
over a period of 4 days to 2.9« The layar at 12.5 km saaaad to persist for 
aavaral waaks. In Figure 7.6(a) the profile of aoatterlng ratio obsarvad on 
June 17t i960 la plotted* Onoe again« a layar at 13 ka is present* The 
•maximum scattering ratio haa daoraasad to 2.0. It is not ulaar whether this 
layar resulted from a aubxaquant eruption on June 13t I960 or was a residual 
from earlier voloanio eruptions. 

A stratoaphario layar of voloanio dust waa obsarvad at 21.5 km on July 25t 
1980 (Figure 7.8(b)). The aoattaring ratio at 21.5 ioi waa 3.0* Tha 
appaaranoa of thia high altitude layar waa oonfirmad by other lidar groups at 
NOAA-Bouldar [Daluiai at al., I960], HASA-Langl«y [HoCwviok at al., I960], 
and tha Onl varsity of Viaoonain [Bloranta, I960]. This high altitude layar 
waa also obsarvad on tha following avaning (Figure 7.6(o))« Tha peak 








scattering ratio haa decreased to 2.1 over the preceding 24 hours. This layer 
is probably associated with another eruption at Nt. St. Helens on July 22 y 
1980. By August 31, 1980, no distinct layers of aerosols are apparent. 
Figure 7.8(d) shows the scattering ratio profile computed tram photooount data 
collected on this evening. The peaks near 10 km are apparently associated 
with sub- visible cirrus clouds. Subsequent observations have also failed to 


reveal evidence of distinct layers. 
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8. CONCLUSIONS AND RECOttlENDATIONS 


8*1 Conclusion 

Laser radar studies have revealed the dynamic nature of the mesospheric 
sodivmt layer. These data consist of photooount profiles whioh are used to 
estimate spatial and temporal profiles of sodium density. Because fAiotooount 
data is subject to the degrading effects of shot noise, appropriate filtering 
techniques are required to increase the spatial and temporal resolution of the 
data. Rowlett and Gardner [1979] developed a spatial processing procedure 
that deconvolved the effects of laser pulse shape and receiver range gate as 
well dis removing shot noise effects and improving spatial resolution. This 
technique has been extended to Include a two-dimensional filtering technique 
which utilizes spatial and temporal filtering to further reduce shot noise and 
Increase profile resolution. The most critical aspect of the two-dimensional 
filtering process is the determination of the appropriate spatial and temporal 
cutoff frequencies. This is accomplished by examining the spatial 2 ud 
temporal perlodograms of the photocount data. Filter cutoffs are selected at 
the frequencies where the signal spectrum no longer dominates the periodogram 
and shot noise contributions are important. Expressions for the expected 
value and variance of the temporal and spatial perlodograms before and after 
spatial and temporal filtering, respeotively, are developed in Chapter 3* In 
addition to their use in determining filter cutoffs, the perlodograms are 
useful in interpreting the photooo'Jit profiles. 

An expression for the layer response to gravity waves has been developed 
that gives insist into the nature of the response and is required to 
correctly interpret the spatial periodogram. In addition to describing 
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spatial and tesporal variations In sodlua density f variations In oolusn 
abundance Induced by gravity wave perturbations are explained. If the layer 
has no horizontal structure, the linear layer response is oharaoterized by 
density oscillations with the sane horizontal wavelength and temporal period 
as the gravity wave. The vertical variations In density, however, are very 
depmdent on the vertioal sodium density gradient. A phase reversal In the 
linear layer response occurs near the layer peak. Below this point 
oscillations in sodium density are 180 degrees out of phase with atmospheric 
density; above this point sodium and atmospheric density perturbations are In 
phase. Because of the vertioal structure of the linear response, a local 
minimum may be observed In the spatial perlodogram at the vertioal wavenumber 
of the gravity wave. Consequently, the layer response must be understood 
before the spatial perlodogram can be Interpreted. 

At points where the linear response Is small or In regions where large 
sodium density gradients exist, nonllnearltles In the layer response become 
important. It Is convenient to use the perturbation series solution presented 
In Chapter 4 to determine the relative magnitudes of the nonlinear 
oscillations in density. The nonlinear terms are associated with sodium 
density oscillations with temporal frequencies which are harmonics of the 
temporal frequency of the gravity wave. Density perturbations consistent with 
nonllnearltles in the layer response have been observed experimentally. 

Lldar observations over Urbana, Illinois are presented which contain 
wavelike features that are apparently a layer response to gravity waves. In 
addition to nonllnearltles in the layer response, which are most readily 
observed In temporal profiles of sodium density, many other layer features 
associated with a layer response to gravity waves are observed. These 
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features Include phase raversals in density perturbations which are apparent 
In both spatial and temporal plots of sodium density and larger density 
perturbations below than above the layer peak. Simulations of the layer 
response based upon theoretical and experimentally observed sodium profiles 
are consistent with observations. Steerable observations have been made at 
Goddard Space Flight Center In Maryland which may be used to deduce the 
horizontal structure of the layer. In addition to Its Influence on the layer 
response to atmospheric waves » the horizontal structure of the layer can be 
used to Infer background winds at mesospheric heights. 

Column abundance variations Induced by gravity wave perturbations In 
sodium density have been observed. In addition) a systematic Increase In 
sodium density above 90 km has been observed during the pre>sunrlse hours. 
The resulting Increase In column abundance Is quite dramatic. In one case the 
column abundance almost tripled In a span of four hours prior to sunrise. The 
Increase Is apparently related to the Influx of atomic sodlui rather than the 
response of the sodium layer to atmospheric wave motions. Meteoric ablation 
Is probably the major source of mesospheric sodium [Hunten, 1981]. It Is 
Interesting that the mean meteor rate and the column abundance both Increase 
In the early morning hours. The meteor rate reaches a maximum at 
approximately 0600 local time [McKinley, 1961]. This Is Illustrated In 
Figure 8.1. Because the Ildar measurements presented here did not extend 
beyond sunrise, the daytime behavior of the sodium column abundance and the 
time of maximum abundance are not known. 

Early dayglow measurements Indicated a marked diurnal variation with the 
oolimm abundance Increasing from sunrise and reaching a maximum at noon. The 
evening abundance declined to the same value as the morning abundance. 


HOURLY METEOR RATES 


Figure 8*1 



LOCAL TIME 


Diurnal variations in nonshower aeteor rates obtained by 
(a) visual observations » (b and o) baoksoatter radar 
observations ( and (d) forward scatter radar observations* 
[McKinley , 1961]* 






DaytlM absorption aeasuremsnta indloatad auoh saallar diurnal variations than 
the dayglow observations. These observations are susaarized in Brown [1973]* 

V 

Daytine Ildar observations tOibson and Sandfordt 1972] failed to reveal any 
substantial diurnal variation* In making dayglow measurwBents it is 

generally assumed that all the observed emission is due to a resonance line 
(l.e.} sodium D1 or D2). This assumption may not be Justified [Grainger and 
Ring, 1962; Noxon, 1968]. If the observed emission is oontaminated by souroes 
other than resonant emlssloni the abundance will be overestimated, nils may 
explain some of the disorepanoles in different experimental methods. Lidar 
observations at Urbana and Goddard Space Plight Center have also revealed an 
early morning increase in column abundance* It is olear that further 
observations are required to resolve the discrepancies* 

8.2 SUGGESTIONS FOR FUTURE WORK 

8*2.1 Experimental Efforts 

In order to resolve questions concerning diurnal variations in sodium 
density) daytime measurements should be made* Hhile the observations reported 
here indicate a pre<-sunrlse enhancement in sodiimi density, daytime Ildar 
observations have failed to reveal any systematic variation [Gibson and 
Sandford, 1972]* Because of the narrow fleld«of-vlew possible with the 
telescope at Goddard Space Flli^t Center, the colleotl(m of high quality 
daytime data is feasible. By narrowing the telescope field-of>vlew and the 
bandpass of the Interference filter in the receiver, the daytime background 
noise level may be reduced sufficiently to permit data collection. 

i 

« Observations spanning the period of morning increase In abundance and 

I extending into the daylight hours should be made. 

i 

r 

i- 
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Additional ataorabla obsarvatloas should bo node* B]r soannlng tho 
tolosoopo in aziauth uhllo oolleotinf data, a sorlos of pi^riloa is obtainod 
whioh desoribes the vertioal and horizontal variations in the sodiw Isfsr as 
tiell as tenporal variations in density. If a suffioient nuaber of 
horizontally si>aoed points are observed, data prooessing aay be extended to 
Include spatial filtering in the horizontal dlreotion. This would yield 
three-dlaensionally filtered data and would further reduoe shot noise and 
increase profile resolution. Profiles of .^odlua density detailing the 
vertioal amd horizontal structure of the sodlia layer as well as teaporal 
variations in density at each observed point could be estlaated. These 
estiaates would be very useful in determining the horizontal structure of the 
layer thus penal ting layer dynamics to be modeled more accurately. A 
horizontal gradient resulting from the morning enhancement in sodium density 
would be revealed in this type of observation. Because a large number of 
{Mints would be scanned, these measurements would require rapid data 
oolleotlon and would have to take advantage of the 10 Hz pulse re|>etltlon rate 
of the Candela laser. 

Correlative observations have been made in whioh the University of 
Illinois meteor radar and the sodium Ildar made near siaultaneou& observations 
[Tetenbaua, 1981]. Mesospheric winds revealed by the meteor radar can be 
ooapared to Ildar observations of sodlua density to study the the effects of 
ataospherio motions on the sodium layer. In effect, two independant 
techniques oan be utilized to estiaate mesospheric winds. Radar aethods, used 
to track the ionized trails left by meteors, and Ildar observations of the 
sodium layer density response to ataospherio waves aay both be utilized to 
deduoe ataos{>herio aotlons at aesospherio heights. Beosuse the lidar is 
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adjacent to the neteor radar « Intarferenoe i>reventa the Ildar froa operating 
Hhen the aeteor radar la tranaalttlng* Thwt, extended » uninterrupted Ildar 
observations oannot be aade when the aettrar radar operates. Consequently) It 
is not praotloal to teaporally filter the i^otooount data and) In general) the 
resulting resolution and short observation periods aalM It dlffloult to 
observe any perlodlo variations In sodliai density. This aakes ooaparlsona of 
aeteor radar winds and sodlua density profiles dlffl<*ult. The oorrelatlve 
observations aade to date are promising but not oonoluslve. Comparisons of 
the two aeasur«Mnts oan be aade aore readily If the Ildar systea oan be 
shielded from the aeteor radar sufficiently to permit slaultaneous operation. 

Sodium temperature measurmnents oan be made. Beoauae omny of the ohemloal 
reactions responsible for the nature of the sodium layer are temperature 
dependent) further Insights Into sodium layer ohemlstry could be obtained from 
these measurements. Another means of monitoring atmospheric motions would be 
provided If the oscillations In temperature associated with the propagation of 
atmospheric waves could be observed. Sodium tmaperature may be Inferred by 
obtaining estimates of the shape of the D2 emission line at 5890 %• 'me 
temperature may be deduced from the doppler broadening of the line. The line 
shape Is be observed by scanning the laser output frequency through the D2 
line. In order to obtain the required spectral resolution) this method 
requires the laser to have a very narrow spectral output (much less than the 
width of the D2 line). Since the D2 line la approximately 1 pm wide) output 
llnewldths near .1 or .2 pm are necessary. In addition) aoourate measurements 
of the center frequency and output llnewldth must be made. Because of laser 
Instabilities) It may prove dlffloult to maintain laser tuning for the tine 
necessary to make extended observations. An alternative approaoh la to soan 
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tiM F«oelv«r tuning aorosn the D2 line while the laser la aligned for a broad 
speotral output. In thla oase, the epeotral width of the laaer output should 
be at least as broad as the D2 line. Consequently, laser tuning la aore 
easily aalntalned. The reoelver efflolenoy aay be reduced taoause of the very 
narrow bandpass and tunabllity requlreaents of this systea. The reoelver 
bandpass oust be aueh leas than the width cf the sodlua 1>2 line, so the 
required bandpass Is roughly *1 or .2 pa. 

8.2.2 Theoretloal Development 

Because the sodlua layer response to gravity waves Is In part nonlinear, 
the layer response to several gravity waves will inolude Interactions between 
the layer responses to Individual waves. The nonlinear Interaction aay aalce 
Interpretation of experiaental data difficult. For exaaple, the data 
collected on March 11-12, 1981 (presented In Chapter 5) contains aany wavellke 
features which cannot be attributed to the layer response to a single gravity 
wave. It Is probable that a spectrua of gravity waves rather than a single 
wave Is usually present. As a result, the layer response to several gravity 
waves should be developed. It Is anticipated that a "beating" effect occurs 
In which density oscillations appear whose teaporal frequencies are the sub 
and difference frequencies of the gravity waves. The effects of the wave 
Interactions on spatial and teaporal perlodograas require Investlgstlon. 

In Chapter 4, the effects of a background wind have been Ignored. If the 
sodlua layer has horizontal structure, the background wind will cause ohanges 
In sodlua density as features are swept overhead. These density ohanges aay 
be used to deduce the aean wind If the horizontal structure of the layer Is 
knowi. This Is one of the aotlvatlons for asking steerable Ildar 
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observations* Thors will bo Intsraotlons bstwson density perturbations 
Induced by ataospherlo waves and density ohanfss assoolated with the 
Interaotlon of a baokground wind and a horizontally struotured layer. Correot 
Interpretation of Ildar observations requires that the layer response to both 
baokground winds and ataospherlo waves be known* 

The Bomlng enhanoeaent observed In data presented In Chapters 5 and 6 la 
not understood* It appears to be assoolated with the Influx of atoalo sodlua 
Into the sodliai layer* The aeohanlsa for this Influx aust bo explained* In 
addition I the Influx af foots the sodlua layer response to ataoaidierlo waves by 
changing the steady state layer profile* In developing the layer response * 
source and loss terns were Ignored In order to amplify the development of 
the solution* It alght be Instructive to Include a source tern which ref loots 
the observed Inoreases In sodlua density during the early aomlng hours* 

The layer density response has been used to Investigate sodlua layer 
dynaalos In terns of a layer density response to gravity waves* Ataospherlo 
tides will also Induce density perturbations* The notion of the layer peak on 
March 16>17t 1981 (with a period of 12 hours) was suggestive of tidal 
aotlvlty* The response of the layer to tidal osolllatlons should be 
developed* In addition f tldal-gravlty wave Interactions in the layer response 
should be Investigated* 

A oonplete study of sodlua layer dynaalos aust Inolude sodlua layer 
obMlstry as well as the oorksorew aeohanlsa. The oorksorew aeohanlas aay be 
mportant only In the upper half of the layer* but ohealoal effeota oould be 
quite mportant near the sink of atoalo sodliai at the layer bottonslde* The 
layer ohealstry and descriptions of the sink for atoalo sodlua are still not 
oonpletely understood* Bven though these effects are negleoted* aany of Uie 
djmaaloal features of the sodliai layer are explained by the analyam m 
Chapter A* 
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APPENDIX I 

OSSCRIPnON OP PHI BLAHKINQ AND TIMINO CONTROL CIRC0IT3 

Beoaua* th« PNT la subject to ovarloadins as a oonsaquanoa of strong 
bsoksoattsp froa bolow 30 Iob» s blanking circuit has boon added to the 
reoelvlott syatea. The PMT gain la reduced by decreasing the first dynode 
voltage to the oathode voltage. The circuit which perfoms this task la given 
In Figure 1.1. Nhen a low TTL level (0 volts) Is present at the blanking 
Input) the output of the opto-laolator (TIL 11?) will be high (5 volts). As a 
result) transistors Q2 and Q3 will be turned off and the first dynode voltage 
will not be Influenced by the gating circuitry* The current through the chain 
of bias resistors In the PMT housing is auoh greater than the leakage current 
froa oolleotor to ealtter In Q3> When the blanking Input la *5 volts ) the 
opto-lsolator output Is near 0 volts and Q2 and Q3 are turned on* The first 
dynode voltage Is pulled do*^n to the oathode voltage and the PMT Is blanked. 

For safety reasons the PMT anode la aalntalned at approxlaately 0 volts 
and the cathode voltage la adjusted to -2000 volts* The opto-lsolator 
protects the tlalng controller and operators froa the hl^ oathode voltage. 
The 5 volt supply which provides the bias voltage for the opto-lsolator and Q2 
Is referenced to the oathode voltage (l.e. the bias voltages are 4>5 volts 
with respect to the PMT cathode). 

The tlalng controller synohronizee laser firing and PMT blanking. Its 
soheaatlo Is given In Figure 1.2. The LM555 operates as a free running 
oscillator and establishes the laser pulse repetition rate. The oscillator 
frequency and) therefore) the laser pulse repetition rate aay be varied froa 
roughly 0*3 to 30 Hz by adjusting RA. An external Input has bean suppllad 
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which enables the Ildar systea ooaputer to ^te the osolllator output. When 
SI la In the "Coaputer Gate On” position and 0 volts Is placed on the ooaputer 
gate Input} the osolllator output Is Interupted and the laser tdll not be 
triggered. A •t-5 volt Input will enable the osolllator and the laser will 
pulse at the rate established by the setting of R4. The osolllator la 
continuously enabled when SI Is In the "Computer Gate Off" position. 

The PMT gating control output Is normally at a 0 volt level. When an 
osolllator pulse triggers the monostable multivibrator (C) and S2 Is In the 
"Gating On" position} a ••>5 volt pulse Is sent to the blanking clroult. This 
blanks the PHI and prevents overloading. The duration of this pulse Is 
nominally 205 yS} but nay be adjusted with R6. If S2 Is In the "Gating Off" 
position} the PMT Is not blanked. However} the laser trigger continues to be 
generated. 

The oscillator pulse which triggers the blanking pulse also triggers 
another monostable (A). The output of this monostable Is a 5 us positive 
pulse. The falling edge of this pulse triggers the remaining monostable (B). 
The output of (B) Is delayed by 5 ys as a result. The duration of the pulse 
generated by (B) IS nominally 6-10 ys. Q1 is used to Increase the pulse 
voltage to approximately 12 volts. The resulting 12 volt} 10 ys pulse 
triggers the laser. The PMT gating control remains at *5 volts for another 
200 ys (as determined Iqr R6) insuring that the PMT Is protected from 
overloading below 30 km. 
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